JLAB-THY-05-371 



Few— Nucleon Forces and Systems in 
Chiral Effective Field Theory 

Evgeny Epelbaum *^ 

Jefferson Laboratory, Theory Division, Newport News, VA 23606, USA 

Abstract 

We outline the structure of the nuclear force in the framework of chiral effective field theory of QCD 
and review recent applications to processes involving few nucleons. 

Keywords: effective field theory, chiral perturbation theory, nuclear forces, few-nucleon systems, 
isospin-violating effects, chiral extrapolations. 



Contents 

1 Introductionl 



2 Effective field theories in nuclear phvsicsl 



2.1 Chh'al perturbation theory 



2.2 EFT for nucleons at very low energy 



3 

3 
9 



2.3 



Chiral EFT for few~nucleon svstemsl 13 



3 Nuclear forces in chiral effective field theory 



18 



3.1 Nuclear potentials from field theory! 18 



3.2 



Two-nucleon forcel 22 



3.2.1 Regularization of the pion-exchange contribution j 22 



3.2.2 



Pion-exchange contributions! 24 



3.2.3 Contact termd 27 



3.2.4 Relatiyistic correctiond 28 

30 



3.2.5 Isospin-breaking corrections! 



3.2.6 



Towards a relativistic NN potential! 35 



3.3 



Few-nucleon forced 37 



3.3.1 Three- and four-nucleon forces in the isospin limit! 37 

3.3.2 Isospin-breaking corrections! 40 



3.4 



Role of the A-excitatioiil 43 



*^email: epelbaum@jlab.org 



4 Applications to few nucleon svstema 46 

4.1 The two-micleon svstemi 46 

4.1.1 Peripheral micleon-micleon scattering! 46 

4.1.2 Regularization and renormahzation of the Schrodinger eauationi 51 

4.1.3 Two nucleons up to N^LOl 53 

4.1.4 Resonance saturation for NN contact interaction j 62 

4.1.5 Quark mass dependence of the nuclear forcel 63 

4.2 The three- and four-nucleon systems up to N^LOl 65 

4.2.1 The formalismi 66 

4.2.2 Elastic Nd scatterind 67 

4.2.3 Nd break-upl 71 

4.2.4 Bound stated 73 

4.3 More micleonsl 74 



5 Miscellaneous omission d 76 



6 OutlooJ 77 



A 2PE potential from the third order TriV amplitudel 77 



1 Introduction 

One of the basic problems in nuclear physics is determining the nature of the interactions between the 
nucleons which is crucial for understanding the properties of nuclei. The standard way to describe the 
nuclear force is based on the meson-exchange picture, which goes back to the seminal work by Yukawa 
01 . His idea as well as experimental discovery of vr- and heavier mesons (p, w, . . .) stimulated the 
development of boson-exchange models of the nuclear force, which still provide a basis for many modern, 
highly sophisticated phenomenological nucleon-nucleon (NN) potentials. 

According to our present understanding, the nuclear force is due to residual strong interactions between 
the color-charge neutral hadrons. A direct derivation of the nuclear force from QCD, the underlying 
theory of strong interactions, is not yet possible due to its nonperturbative nature at low-energy. In 
order to provide reliable input for few- and many-body calculations, research on the NN interaction 
proceeded along (semi-)phenomenological lines with the aim of achieving the best possible description of 
the low-energy NN data. In the case of two nucleons, the potential can be decomposed in a few different 
spin-space structures, and the corresponding radial functions can be parameterized using an extensive 
set of data. Although the resulting models provide an excellent description of experimental data in 
many cases, there are certain major conceptual deficiencies that cannot be overcome. In particular, one 
important concern is related to the problem of the construction of consistent many-body forces. These 
can only be meaningfully defined in a consistent scheme with a given two-nucleon (2N) interaction 
Notice that because of the large variety of different possible structures in the three-nucleon force (3NF), 
following the same phenomenological path as in the 2N system and parametrizing the most general 
structure of the 3NF seems not to be feasible without additional theoretical guidance. Clearly, the same 
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problem of consistency arises in the context of reactions with electroweak probes, whose description 
requires the knowledge of the corresponding consistent nuclear current operator. Further, one lacks 
within phenomenological treatments a means of systematically improving the theory of the nuclear force 
in terms of the dominant dynamical contributions. Finally, and most important, the phenomenological 
scheme provides only loose connection to QCD. 

Chiral Effective Field Theory (EFT) has become a standard tool for analyzing the properties of hadronic 
systems at low energies in a systematic and model independent way. It is based upon the approximate 
and spontaneously broken chiral symmetry of QCD, which governs low-energy hadron structure and 
dynamics. In addition, it provides a straightforward way to improve the results by going to higher 
orders in a perturbative expansion. In the past two decades, this framework was successfully applied to 
a variety of low-energy reactions in the meson and single-baryon sectors. Fifteen years ago Weinberg 
0, 0| proposed a generalization of this approach to the few-nucleon sector, where one has to deal with a 
nonperturbative problem. He demonstrated that the strong enhancement of the few-nucleon scattering 
amplitude arises from purely nucleonic intermediate states and suggested to apply EFT to the kernel 
of the corresponding scattering equation, which can be viewed as an effective nuclear potential. This 
idea has been explored in the last decade by many authors. In this work, we will review the current 
status of research along these lines, focusing on the description of few-nucleon systems. The manuscript 
is organized as follows. In section El we discuss effective field theories which are of relevance for the 
topics considered in this work. In particular, we give a brief account of chiral EFT for the pion and 
single-nucleon sectors and then discuss how it can be generalized to the few-nucleon sector. In section |31 
we consider the structure of the nuclear force based on chiral EFT. Applications to systems with 2 ... 7 
nucleons are presented in section |1J Section [S] lists some related further topics which are not covered in 
this work. A brief outlook is presented in sectional Finally, appendix^ contains explicit expressions for 
certain contributions to the 2N force. 



2 Effective field theories in nuclear physics 
2.1 Chiral perturbation theory 

Chiral Perturbation Theory (CHPT) is the effective theory of QCD and, more generally, of the Standard 
Model, which was formulated by Weinberg and developed in to a systematic tool for analyzing low- 
energy Quantum Chromodynamics (QCD) by Gasser and Leutwyler |a, l3|. In this section, we give a brief 
overview of the foundations of this approach. 

Consider the QCD Lagrangian in the two-flavor case of the light up and down quarks 

>Cqcd = q {il^^D^ - M)q - \g%G-^'' , (2.1) 

where = d^ — igsG'^T'^ with T", (with a = 1 ... 8) are the SU(3)coior Gell-Mann matrices and q the 
quark fields. Further, G^^, are the gluon field strength tensors, and the quark mass matrix is given by 
M = diag(mu, m^). We do not show in Eq. (|2.H) the 6- and gauge fixing terms which are not relevant 
for our consideration. The left- and right-handed quark fields are defined by gL,R = 1/2(1 ± 75)?- The 
chiral group G is a group of independent SU(2)fiavor transformations of the left- and right-handed quark 
fields, G = SU(2)l x SU(2)r. Expressing the quark part in the QCD Lagrangian (|2.1() in terms of qL,R-, it 
is easy to see that the covariant derivative term is invariant with respect to global chiral rotations, while 
the quark mass term is not. The running quark masses in the MS scheme at the renormalization scale 
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H = 1 GeV are rriu ~ 5 MeV and rrid 9 MeV 8]. Given the fact that the masses of the up and down 
quarks are much smaller than the typical hadron scale of the order of 1 GeV, chiral SU(2)l x SU(2)r 
symmetry can be considered as a rather accurate symmetry of QCD. 

There is a strong evidence on both experimental and theoretical sides that chiral symmetry of QCD is 
spontaneously broken down to its vector subgroup (isospin group in the two-flavor case). Perhaps, the 
most striking evidence of the spontaneous breaking of the axial generators is provided by the nonexistence 
of degenerate parity doublets in the hadron spectrum and the presence of the triplet of unnaturally 
light pseudoscalar mesons (pions). The latter are natural candidates for the corresponding Nambu- 
Goldstone bosons which acquire a small nonzero mass due to the explicit chiral symmetry breaking by 
the nonvanishing quark masses. Further, on the theoretical side, recent (quenched) QCD determinations 
of the vacuum expectation value of the scalar quark condensate S = (0|gg|0), a natural order parameter 
of the spontaneous chiral symmetry breaking, yields the nonvanishing value 

S = -(262 ±12MeV)^. (2.2) 

This value is based on the MS scheme at the renormalization scale /x = 2 GeV. Also based on rather general 
arguments, it has been shown that the vector subgroup of the chiral group cannot be spontaneously broken 
jlOt |. Further, in the three-flavor sector, spontaneous chiral symmetry breaking is consistent with the 
so-called anomaly matching condition 11]. For arguments based on the large Nc limit the reader is 
referred to Ref. |12i] . 

The low-energy dynamics of QCD can be studied using the method of external sources 0, [^. The idea 
is to couple quarks to external classical fields which formally allows to compute Green functions of the 
corresponding quark currents in a straightforward way. The extended QCD Lagrangian takes the form 

Cqcb = >C^cD + 97m i^'' + 750^) q-qis- «75P) q , (2.3) 

where the external fields u^, a^, s and p are hermitian, color neutral, traceless 2x2 matrices in flavor 
space and Cqqy) refers to the Lagrangian in Eq. (|2.1j) in the absence of the quark mass term. Notice 
that one can also include an external vector singlet field which becomes particularly useful for studying 
electromagnetic processes. The transformation properties of external sources follow from the require- 
ment that the extended QCD Lagrangian is invariant under local chiral rotations. The ordinary QCD 
Lagrangian is recovered by setting = = p = 0, s = diag(mu, m^). The QCD Green functions 
built from the associated quark currents can be derived by taking functional derivatives of the generating 
functional T{v,a, s,p) defined as 

with respect to the sources. It is not presently possible to evaluate the Green functions in a closed 
form. At low energy, however, one can calculate the generating functional within effective field theory 
formulated in terms of the observed asymptotic states. As proven by Leutwyler j^sf, the gauge-invariant 
generating functional can be represented by a path integral constructed with a gauge-invariant effective 
Lagrangian for the Goldstone bosons Ce{f{U; v,a,s,p): 

^r[v,a,s,p] _ h^jj-^^J id^xCcsiU;v,a,s,p) ^ ^2.5) 



Here, the 2x2 unitary matrices U, which satisfy det U = 1, collect the triplet of pseudoscalar Goldstone 
bosons. The Green functions can be evaluated from the effective Lag rangian >Ceff i^ ^ systematic way 
by expanding in powers of the external momenta q and the quark mass matrix A4 and keeping the ratio 
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M./cf' fixed. This procedure to evaluate S'-matrix elements is called chiral perturbation theory and can, 
in principle, be carried out to arbitrarily high orders in the low-energy expansion. Since, at present, 
one cannot derive £efr from QCD directly, one writes down its most general form including all terms 
consistent with the required symmetry principles. This can be done following the lines of 3 3l- '^^^ 
effective Lagrangian takes the form 

£eff = 4'^+4')+4') + ..., (2.6) 

where the superscripts d = dg -\- 2dM refer to the number of derivatives dg and/or quark mass matrices 
dM- To be specific, let us parametrize the matrix U as 



U(x) = exp 



r • 7V(x } 
i- 



(2.7) 



F 

where Tj are the Pauli isospin matrices, tTj are pion fields and F is a constant. Notice that one also often 
uses the so-called sigma-representation U{x) = \/l — 7v'^{x) / F"^ + iT ■7v{x)/F. The matrix U is required 
to transform under local chiral rotations as 



U RUL^ . {21 



Here, the 2x2 matrices L and R represent local SU(2)l and SU(2)r rotations: L{x) = exp[— iT-0L(x)/2], 
R{x) = exp[—iT-6ji(x}/2]. One can show from Eqs. 1)2. 8|) . 1)2. 7|) that pions belong to a nonlinear realization 
of the chiral group [l4l] and transform linearly under its vector (isospin) subgroup given by a subset of 
rotations with Ol = Or. Different realizations of the chiral group turn out to be equivalent to each other 
by means of nonlinear field redefinitions [3] • It is convenient to define the quantity u 



u 



2-C/, u^Ruh\ (2.9) 



where the 2x2 matrix h{x) depends on L, R and U . More precisely, it is given hyh = VLWR^RVU. 
The effective Lagrangian is constructed out of the following building blocks, see e.g. |lfi] : 



X- = 



: u^D^'xu^ -uD^'x^u, 
= u^Fj!^''u±uF^''uK (2.10) 

Here, x = 2B{s + p) with B being a constant, F/"" = d^'F^ - d^Fj^ - i[F^, Ff] with I = L,R is the 
field strength tensor associated with F^ = + a^^, F^ = — a^, and the covariant derivative D^X is 
defined via 

D^X = d^X - i{v^ + a^)X + iX{v^ - a^) . (2.11) 

All quantities in Eq. 1)2. lUj) transform covariantly under G, i. e. I hlh^ . Chiral invariant terms in the 
Lagrangian can therefore be easily constructed via building the traces (denoted in the follo wing by (...)) 
of the products of these objects. The leading and subleading Lagrangians take the form Il6|| : 

4'^ = ^{u,u^ + X+), 

+ ^|(/^K,n.])-i^(X-)^ (2.12) 



5 



where h,...,7 are low-energy constants (LECs). Further, the constant F can be identified with the pion 
decay constant in the chiral hmit while the constant B is related to the scalar quark condensate via 
(0|un|0) = {0\dd\0) = —BF"^ + 0[M). Notice that we are followirig here the standard CHPT scenario 
with 2B{mu + md)/M^ ~ 1. The generalized CHPT scenario jjlJl, in which 25(m„ + md)/M^ <C 1, 
is ruled out by the recent determination of the S-wave, isospin-zero vrvr scattering length Cq from the 
kaon decay K — > Trneu [l^ . [l9| . We stress that there are further terms in which do not contain 
Goldstone boson fields and are therefore not directly measurable. In addition, one has also to account for 
the chiral anomaly which can be done along the lines of Ref. 0, see also Q. The effective Lagrangian 
in Eqs. H2.12() can be used to describe the interaction between pions among themselves and between pions 
and external fields in the low-energy regime. For a reaction involving Nj^ external pions, the transition 
amplitude M is related to the S'-matrix via S = S'^{pi +P2 + ■ ■ • +P7V^)-^- The low-momentum dimension 
of M, i.e. the power of a soft scale Q associated with the pion mass or external momenta, is given by |^ 

i^ = 2 + 2L + ^V^{d,-2), (2.13) 

i 

where L {V^) is the number of loops (vertices of type i). The chiral dimension di is given by the number 
of derivatives and/or quark mass insertions. Diagrams with loops and/or vertices with more derivatives 
and pion mass insertions are therefore suppressed by powers of Q/A^ with being a hard scale, which 
is sometimes referred to as the chiral symmetry breaking scale. This scale sets the (maximal) range of 
convergence of the chiral expansion. The appearance in the spectrum of the p, the first meson of the 
non-Goldstone type, suggests A-,^, ~ Mp ~ 770 MeV. Another estimate based on consistency arguments 
|23 | leads to A^ ~ AttFj^ with F^^ = 92.4 MeV being the pion decay constant. The leading term in the 

(2) 

low-energy expansion of the scattering amplitude results by evaluating tree diagrams with£r. The first 
corrections arise from tree graphs with exactly one insertion from and one-loop diagrams with all 

(2) 

vertices from Cn . They are suppressed by two powers of momenta or one power of the quark masses 
compared to the leading terms. Notice that all ultraviolet divergences in loop diagrams are cancelled by 
counterterms from . The divergent parts of the LECs /j have been worked out in 6] using the heat- 
kernel method. The finite parts of the /j's are not fixed by chiral symmetry and have to be determined 
from measured data. This then allows one to make predictions for other observables. At next-to- 
next-to-leading order, one must include tree diagrams with one insertion from (and all remaining 
vertices from C^^^) or two insertions from as well as one-loop graphs with a single insertion from 
>ci^^ and two-loop graphs with all vertices from The Lagrangian contains 53 independent 

terms plus 4 terms depending only on external sources and 5 terms (in the absence of a singlet external 
vector current) of odd intrinsic parity [2^ . [3, l2^ . [2^. The renormalization at this order is carried 
out in I24II . At present, several two-loop calculations (i.e. at order p^) have already been performed, see 
e.g. j27i.l28i| for some examples. One of the most impressive theoretical predictions is given by the precision 
calculation of the isoscalar S-wave vrvr scattering length Og , a fundamental quantity that measures explicit 



breaking of chiral symmetry. The results of the two-loop analysis 2^, sO, ^ combined with dispersion 



relations in the form of the Roy equations \32\ allowed for an accurate prediction: Oq = 0.220 it 0.005 

^0 



33l |. To compare, the leading-order calculation by Weinberg yielded Og = 0.16 while the next-to- 



leading value obtained by Gasser and Leutwyler is ag = 0.20 6]. The results of the E865 experiment at 
Brookhaven beautifully confirmed the prediction of the two-loop analysis of Ref. [s^l yielding the value 
ag = 0.216 ± 0.013 (stat) ± 0.002 (syst) ± 0.002 (theor) |li,E3- 

It is clear that CHPT can, in principle, be carried out to an arbitrarily high order in the low-energy 
expansion. The predictive power is, however, limited due to the rapid increase of the number of new 
LECs.'^^ It is, therefore, particularly important to be able to estimate the values of the LECs. One 

*^ Clearly, not all LECs contribute to a particular process/observable, so that there is usually no need to determine all 
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possible way is to assume that the LECs are saturated by low-lying resonances such as the triplet of 
jO-mesons. The form of their couphng to Goldsone bosons is dictated by chiral symmetry and can be 
parametrized in terms of a few parameters. At low energy, the resonance fields can be integrated out 
which gives rise to a series of terms in the effective Lagrangian whose strength is given in terms of 
resonance couplings and masses and can be used as an estimation for the correspondin g L ECs. For more 



details on resonance saturation the reader is referred to |35|, |3g, |33] , see also |38|, |39|, |40|, |41( for recent 
works on meson resonances in the chiral EFT framework. We also emphasize that the LECs in the 
effective Lagrangian are, in principle, calculable in QCD, see 0, 0, |4^, ,42|] for some recent attempts 
using the framework of lattice gauge theory. 

So far we have only discussed interactions between Goldstone bosons and external fields. We will now 
consider the extension of CHPT to the single-nucleon sector. We enlarge the effective Lagrangian 

C-eS = + C,-kN , (2-14) 

to include terms which couple mesons to nucleons. It it is convenient to introduce the nucleon field N in 
the isospin-doublet representation which transforms under the chiral SU(2)l x SU(2)r group as 

N ^ h{L,R,U)N , (2.15) 

with the matrix h being defined according to Eq. 1)2. 9() . The above equation together with Eq. 1)2. 8|) 
specifies the nonlinear realization of the chiral group in terms of pions and nucleons. Notice that for 
vector-like transformations with Oy = Ol = ^i?; the matrix h does not depend on U and reduces to 
the usual isospin transformation matrix h = exp[— zr • 0y/2]. Eq. (|2.15|) is, therefore, consistent with 
the transformation properties of the nucleon field under isospin group. We stress that there is no loss of 
generality in the requirement for the nucleon field to transform under G according to Eq. 1)2. 15|) . Different 
realizations of the chiral group can be reduced to the one specified in Eqs. ¥2^^ by means of field 

redefinitions , . The covariant derivative of the nucleon field is given by 

D^N = d^N + T^N, T^ = ]^[u\ d^,u]-\\v^ + a^)u-\{v^-a^)v) . (2.16) 

G 

The so-called chiral connection F^ ensures that L)^ transforms covariantly under G, i.e.: hD^h!^ . 

To construct the effective Lagrangian ^CttAt one simply combines and the building blocks in Eq. (|2.1U|) , 
which transform covariantly under G, with the appropriate nucleon bilinears. To first order in the 
derivatives, the most general pion-nucleon Lagrangian takes the form ji^l 

6^1 = N (^i^D^ -m + ^^^j,u^) N , (2.17) 

where m and are the bare nucleon mass and the axial-vector coupling constant. Further, the super- 
script of CttN denotes the power of the soft scale Q related to a generic nucleon tree-momentum, pion 
four-momentum or pion mass. Notice that contrary to the pion mass, m does not vanish in the chiral limit 
and introduces an additional large scale. Consequently, terms proportional to D n an d m in Eq. (|2.17|1 
are individually large. It can, however, be shown that {ij^^D^ — m)N ~ 0{Q) |4q|. The presence of 
the additional hard scale associated with the nucleon mass makes the power counting significantly more 
complicated since the contributions from loops are not automatically suppressed. To see this consider 



the correction to the nucleon mass tjin due to the pion loop which in the chiral limit takes the form 42 1 
(ruN - m)L% ^= ° + ^ 1" ^) + - 4) ' (2-18) 



LECs at a given order. 
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where /U is the mass scale introduced by dimensional regularization (DR), d is the number of dimensions 
and the quantity L is given by 



^ = f6^{d^"^^^^^^^^+^'^^) + ^)} ' r'(l) = -0.577215... . (2.19) 

The result in Eq. (|2.18j) shows that mjv receives an (infinite) contribution which is formally of the order 
~ ni (m/'iirF)'^ and is not suppressed compared to m. In addition, the parameter m in the lowest- 
order Lagrangian C^j^ needs to be renormalized. These features of the relativistic ttN EFT should be 
contrasted with the purely mesonic sector where loop contributions are always suppressed by powers 

(2) 

of the soft scale and the parameters F and B in the lowest-order Lagrangian ' remain unchanged 
by higher-order corrections (if dimensional regularization is applied). This problem with the power 
counting in the baryonic sector can be dealt with using the heavy-baryon formalism which 
is closely related to the nonrelativistic expansion due to Foldy and Wouthuysen j5]|. The idea is to 
decompose the nucleon four-momentum according to 

Pf, = mv^ + , (2.20) 

with Vfj^ the four-velocity of the nucleon satisfying = 1 and /c^ its small residual momentum, v-k <^ m. 
One can now decompose the nucleon field in to the velocity eigenstates 

Ny = e^^'-^P+A^ , K = e*'"''-^p- , (2.21) 

where = (1 ±jfj,v^)/2 denote the corresponding projection operators. Notice that for the particular 
choice Vf^ = (1,0,0,0), the quantities and coincide with the usual large and small components 
of the free positive-energy fields (modulo the modified time dependence), see e.g. jH3]- One, therefore, 
usually refers to Ny and hy as to the large and small components of A^. The relativistic Lagrangian C^^lf 



in Eq. (|2.17|) Ccin be 6xprGSSGci in terms of Ny and. hy hs'. 

4iv = NyAN^ + KBNy + KloB^oK - KCK , (2.22) 

where 

A = i{v-D)+gA{S-u), B = --1^ 2i{S ■ D) + ^-{v u) , C = 2m + i{v D) + gA{S ■ u) . (2.23) 

Here = i^f^a^yv'^ refers to the nucleon spin operator. One can now use the equations of motion for 
the large and small component fields to completely eliminate hy from the Lagrangian. Utilizing the more 
elegant path integral formulation [s^l, the heavy degrees of freedom can be integrated out performing 
the Gaussian integration over the (appropriately shifted) variables hy^ hy. This leads to the effective 
Lagrangian of the form 

reS _ AT 
'-■kN — 



A + {joBho)C-^B] Ny = Ny [i{v ■ D) + gA{S ■ u)] Ny + O (-\ . (2.24) 
J \m J 



Notice that the (large) nucleon mass term disappeared from the Lagrangian, and the dependence on m 
in resides entirely in new vertices which can be classified according to their powers of 1/m. Clearly, 
the formalism outlined above can be extended to the relativistic pion-nucleon Lagrangian beyond the 
leading order in derivatives/quark masses. The resulting heavy-baryon Lagrangian can be expressed as 

C^N = + 4'i + 4'i + • • • > (2-25) 



^■^This statement applies for dimensionally regularized loop integrals. 
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where /^^^jy given by the terms in the right-hand side of Eq. H2.24() and the superscripts refer to the 
power of the soft scale Q. Higher-order terms in the Lagrangian will be discussed in sectional We stress 
that in the single-nucleon sector, relativistic corrections are usually treated on the same footing as the 
corresponding chiral corrections, i.e. one counts 1/m ~ Notice further that some of the 1/m-terms 

in the heavy-baryon Lagrangian are protected from extra counterterm contributions as a consequence 
of the so-called reparametrization invariance associated with the freedom in parametrizing the nucleon 
momentum p^. It relies on the fact that the same physics should be described using = mv'^ + k'^, 
v'^ = 1, instead of Eq. (|2:2n|l ^,^,^]. 

The advantage of the heavy-baryon formulation of CHPT (HBCHPT) compared to the relativistic one 
can be illustrated using the example of the leading one-loop correction to the nucleon mass 

(m^ - m)HB = _4ciM^^ - , (2.26) 

(2) 

where the counterterm contribution oc ci stems from C^^f Contrary to the relativistic CHPT result 
in Eq. (|2.18j) . the loop correction in HBCHPT is finite (in DR) and vanishes in the chiral limit. The 
parameters in the lowest-order Lagrangian do not get modified due to higher-order corrections which are 
suppressed by powers of Q/^.^. Notice further that the second term in Eq. (|2.26j) represents the leading 
contributions nonanalytic in quark masses and agrees with the relativistic CHPT result 47 1, see js^ for 
an earlier determination of this correction. In general, the power of a soft scale Q for the scattering 
amplitude in the single-nucleon sector HBCHPT is given by 

u = l + 2L + Y. ^^^^^ - 2) + E ^^"'^('^^ - 1) ' (2.27) 

i i 

where V^^ is the number of vertices from Ct^n with the chiral dimension di. Notice that no closed 
fermion loops appear in the heavy-baryon approach, so that exactly one nucleon line connecting the 
initial and final states runs through all diagrams in the single-baryon sector. 

While most of the calculations in the single-nucleon sector have so far been performed in HBCHPT, it was 
realized a few years ago that its range of convergence is rather limited in some kinematical regions. The 
problem can be traced back to the fact that certain analytical properties of the relativistic amplitude are 
destroyed in the heavy-baryon approach, see e.g. [5^]. This can be avoided using a manifestly Lorentz 
invariant formulation. Various methods like e.g. the infrared regularized CHPT have been developed 
which allow one to stay covariant and, at the same time, to preserve a consistent power counting j58l . 

In the formulation of 58], this is achieved by keeping the infrared singular contributions 
of the loop integrals and simultaneously discarding the polynomial terms that are responsible for the 
breakdown of the power counting and which can be absorbed by local counter terms. More details on 
the foundations and th e appli cations of CHPT in the meson and single-nucleon sectors can be found in 
the review articles lecture notes and recent conference proceedings 0,0. 

A pedagogical introduction is given in 0]. 



2.2 EFT for nucleons at very low energy 

So far we have only dealt with the low-energy processes in the mesonic and single-baryon sectors. Per- 
turbation theory works well in these cases due to the fact that Goldstone bosons do not interact at 
vanishingly low energies in the chiral limit. In the few-nucleons sector one has to deal with a nonper- 
turbative problem. Indeed, given the fact that there are shallow few-nucleon bound states, perturbation 
theory is expected to fail already at low energy. To understand how this difficulty can be handled in 
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the EFT framework it is instructive to look at the two-nucleon system in the kinematical regime where 

Q < 0, El, Ei, Ei, E3- 

Then, no pions need to be taken into account explicitly, and the only 
relevant degrees of freedom are the nucleons themselves. The corresponding EFT is usually referred to as 
pionless EFT. The most general effective Lagrangian consistent with Galilean invariance, baryon number 
conservation and the isospin symmetry takes in the absence of external sources the following form: 



/ 1 1 

C = NUido + ] N - -Cs (iV"^iV)(Af"^iV) - -Ct {N'<aN){N^aN) + ... , 



(2.28) 



where Cs,t are LECs and ellipses denote operators with derivatives. Isospin-breaking and relativistic 
corrections to Eq. (|2.28|) can be included perturbatively [t^. Notice further that in certain cases it turns 
out to be convenient to introduce, in addition to the nucleon field, the auxiliary "dimeron" fields with the 
quantum numbers of the two-nucleon system [t^.*^ Let us consider NN scattering in the ^Sq channel. 
The S'-matrix can be written as 

^2i5 



S 



1 



km 



(2.29) 



where k is the magnitude of the nucleon momentum in the center-of-mass system (CMS) and b (T) is 
the phase shift (T-matrix). Utilizing the effective range expansion (ERE) for (A; cot 6), the T-matrix can 
be expressed as 



^ 47r 
T = 



1 



m k cot 6 — ik 



47r 



m 



1 



+ |roA;2 + V2k'^ + vsk^ + ...)- ik 



(2.30) 



where a, rg and Vi are the scattering length, effective range and shape parameters, respectively. While 
the effective range is bounded from above by the range R of the nuclear potential, the scattering length 
can take any value. In particular, it diverges in the presence of a bound state at threshold. It is then 
useful to distinguish between a natural case with \a\ ~ R and an unnatural case with \a\ ^> R, where the 
range of the nuclear potential is of the order M^^. In the natural case, the T-matrix can be expanded 
in powers of k as: 



T = To + Ti + To + . . . 



m 



1 



iak + ( — - 
V 2 



(2.31) 



A natural value of the scattering length implies that there are no bound states close to threshold. The 
T-matrix can then be evaluated perturbatively in the EFT provided one uses a regularization and sub- 
traction scheme that does not introduce an additional large scale. A convenient choice is given by DR with 
the minimal or the power divergence subtraction (PDS) 7^, 7^ or momentum subtraction at k"^ = — /i^ 
(t^ . In the PDS scheme, the power law divergences, which are normally discarded in DR, are explicitly 
accounted for by subtracting from dimensionally regulated loop integrals not only l/{d — 4)-poles but 
also e.g. l/{d — 3)-poles. The typical loop integral takes then the form 7^, 7^: 



4-d 



m 



d^4 



{2n) 



d-l 



cp' + ie 



m 
47r 



(2.32) 



where p = \p\, q = \q\. The choice /i = leads to the result of the minimal subtraction scheme (MS). 
Taking fi ^ k <^ M^, the leading and subleading terms Tq and Ti are given by the tree- and one-loop 
graphs constructed with the lowest-order vertices from Eq. 1)2. 28(1 . T2 receives a contribution from both 
the two-loop graph with the lowest-order vertices and from the tree graph with a subleading vertex 
"tsI . 7^ . Higher-order corrections can be evaluated straightforwardly. Matching the resulting T-matrix 



^''The auxiliary dimeron fields can be integrated out which leads to a completely equivalent form of the EFT with only 
nucleonic degrees of freedom. 
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Figure 1: The leading and subleading contributions to the S-wave T-matrix in the case of unnatu- 
rally large scattering length. Solid dots (filled rectangles) refer to contact vertices without (with two) 
derivatives. Lines represent the nucleon propagators. 



to the ERE in Eq. 1)2. 31() order by order in the low-momentum expansion allows to fix the LECs Cj. At 
next-to-next -to-leading order (N'^LO), for example, one finds: 



Co 



At: a 



m 



1 + 0{an) 



C2 



27ra^ 



m 



(2.33) 



where the LECs Co and C2 are defined Co + C2A;2 + .... The LEG 

Co is related to Cs,t in Eq. as C = C5 - 3Ct. We stress that the ERE in Eq. can also be 

reproduced in the cut-off EFT framework choosing A ~ and resumming loop diagrams to all orders 
(i.e. solving the Lippmann-Schwinger equation with the nuclear potential given by contact interactions). 

For the physically interesting case of np scattering, the two S-wave scattering lengths take unnaturally 
large values: 



ai 



So 



-23.714 fm 



-16.6M. 



-1 



5.42 fm ~ 3.8 M 



-1 



(2.34) 



Instead of using the low-momentum representation in Eq. (I2.3ip which is valid only for k < 1/a, one can 
expand the T-matrix in powers of k keeping a/c ~ 1 ^75, 7(i]: 



4lT 



T0 + T1 + ... 



(2.35) 



ro 



+ 



V2 



k' + 



m {a-^+ik) ' 2{a-'^ + ik) \'i{a-^ + ik^ ' {a'^ + ik) ^ 

The EFT expansion of the T-matrix in the unnatural case is illustrated in Fig. ^ The leading term 
r_i results from summing an infinite chain of bubble diagrams with the lowest-order vertices. The 
corrections are given by perturbative insertions of higher-order interactions dressed to all orders by the 
leading vertices. Matching the resulting T-matrix with the one in Eq. 1)2. 35() one finds at NLO: 

„ 47r 1 „ Att \ rn , , 

Co = — — , C2 = — 7—^ —-^ . (2.36) 



m a 



m (a ^ — n)^ 2 



Notice that since fi ^ k and <C the LECs Co, 2 scale as Co ~ 1/A; and C2 ~ 1/k'^. More generally, a 
LEG C2n accompanying a vertex with 2n derivatives can be shown to scale as C2n ~ l/k""^^ [l^jEHl- This 
has to be contrasted with the scaling C2n ~ in the case of a natural scattering length, cf. Eq. ()2.33() . 
Notice further that the LECs Co 9 ta ke very large values in the MS scheme (i.e. for = 0) which destroys 
the manifest power counting 



OO P la 

mm 



The three-nucleon problem within pionless EFT has attracted a lot of scientific interest during the past 



few years, see e.g. [80|, |8l|, |8j, |83, |8J, |85|, |8a, |83|- The ultimate question is to what extent the low- 



energy behavior of the 2N system constrains the properties of the three-nucleon (3N) system. Here, it is 
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Figure 2: The first line is the graphical representation of the integral equation describing nucleon-dimeron 
scattering. The second line shows the dressed dimeron propagator. Double (bold) lines correspond to the 
dressed (bare) propagator of the dimeron. Shaded rectangles refer to the 3N T-matrix. For remaining 
notation see Fig. ^ 



particularly interesting that one can identify universal properties of systems, where the scattering length 
in the two-body system is large. This situation is not only realized in the NN system, but also for '^He 
atoms and atomic systems close to a Feshbach resonance, see (s^ for more details. The integral equation 
for the T-matrix describing nucleon-dimeron scattering and including the leading three-nucleon force 
is schematically depicted in Fig. [2 For simplicity, we here discuss the case of three interacting bosons, 
which gathers the main aspects of the problem. For the state with total orbital angular momentum 
L = 0, the integral equation takes the following form in the three-body CMS: 

16 4 /"^ 1 

T{k,p;E) = —M{k,p;E) + - dqq^T{k,q; E) =M{q,p;E), (2.37) 

3a vr Jq _a-i + ^/3q^/A - niE - le 

where the inhomogeneous term reads 

1 , f k"^ + kp + - mE\ H(A) 
M(k,p; E) = _ in (^^-^^-^^ ^ . (2.38) 

Here a (H) is the two-body scattering length (the strength of the three-body force), k = \k\ (p = \p\) 
is the magnitude of the dimeron incoming (outgoing) momenta and E = ^ — 7 is the total energy in 
the incoming state with 7 ~ (ma^)~^ being the two-body binding energy. The incoming and outgoing 
bosons are taken on the energy shell. The on-shell point corresponds to k = p and the phase shift can 
be obtained via 

T{k,k;E). (2.39) 



k cot 6 — ik 

For H = 0, A ^ CO, Eq. ()2.37|) has been first derived by Skorniakov and Ter-Martirosian [s^. It is 
well known that Eq. (|2.37|) has no unique solution in this limit [9o[-'^^ The regularized equation has a 
unique solution for any given (finite) value of the ultraviolet cut-off A but the amplitude in the absence 
of the three-body force shows an oscillatory behavior on A. Cut-off independence of the arnplitude 
is restored by an appropriate "running" of H{A) which turns out to be of a limit cycle type [sol . Isif . 
Adjusting H(A) to a single three-body observable for large enough A (or even in the limit A 00) 
allows to determine all other low-energy properties of the three-body system. Alternatively, this can 
also be achieved by choosing H = 0, tuning A to reproduce a three-body data point and using the same 
cut-off to calculate other observables jo^l- It has also been conjectured that the behavior of the physical 
amplitude at asymptotically large momenta has to satisfy certain constraints which might be used to 
extract a unique solution of Eq. 1)2. 37() in the case H = and A — > 00 8^, 9^ . 



*^Whether Eq. l|07|l with // = and A 00 possesses a unique solution or not depends on the value of the factor which 
multiplies the last term in Eq. 112.371 . 
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The 3N problem can be considered as generalization of the bosonic case. For S-wave Nd scattering in 
the spin-3/2 channel, the corresponding equation for the T-matrix has a unique solution for H = 
and A — > oo so that there is no need to include a 3N force. In the spin-1/2 channel two nucleons can 
form both spin~l and spin-0 dimeron fields which leads to a pair of coupled integral equations for the 
Nd r~matrix. Including the leading contact 3N force allows one to solve these equations in a cut-off 
independent way [s^]- Thus, one needs a new parameter which is not determined in the 2N system in 
order to fix the (leading) low-energy behavior of the 3N system in this channel. Higher-order corrections 
to the amplitude including the ones due to 2N effective range terms can be included perturbatively*^ 
Extension to 3N channels with different quantum numbers is straightforward [s^ . For 
the current status of these applications see [o^, [ol, |9^ and references therein. Universal low-energy 
properties of few-body systems with short-range interactions and large two-body scattering length are 
reviewed in [s^, see also 97 1 for an early work on this subject. First results in the four-body sector 
within pionless EFT are presented in [Sa, Recently, this approach has also been applied to halo 
nuclei, see for an overview. For more details on these and further topics including applications to 
a variety of electroweak processes in the 2N sector see the recent review articles 9^, 9^ and references 
therein. 



2.3 Chiral EFT for few— nucleon systems 

So far we have considered few-nucleon processes at very low momenta k <^ which can be well treated 
within pionless EFT. We now wish to go to higher momenta k ~ M-,^ where the inclusion of explicit pions 
is mandatory. The interaction between pions and nucleons is governed by the spontaneously broken 
approximate chiral symmetry of QCD as explained in section ITTl One would, therefore, like to have an 
approach which utilizes both resummation of certain classes of Feynman diagrams in order to describe 
the nonperturbative features of few-nucleon systems as well as chiral expansion familiar from the single- 
nucleon sector. While the leading NN contact interaction has to be resummed to all orders at least in 
the case of an unnaturally large scattering length, it is not clear a priori whether the interaction resulting 
from the exchange of pions between the nucleons is weak enough to be treated perturbatively. We will 
now outline two basic EFT appro aches with explicit pions to few-nucleon systems: the one due to Kaplan, 
Savage and Wise (KSW) [75|, |7a| which treats pion exchange in perturbation theory, and t he ot her one 
due to Weinberg 0, Q| based on its nonperturbative treatment. For yet another scheme see 

The KSW formalism represents a straightforward generalization of the pionless EFT approach for the 
case of large scattering length discussed in section to perturbatively include diagrams with exchange 
of one or more pions. The scaling of the contact interactions is assumed to be the same as in pionless EFT 
(provided one uses DR with PDS or an equivalent scheme to regularize divergent loop integrals). For 
k ~ M-n- ~ a~^, the leading-order S-wave amplitude T_i is still given by the diagrams shown in the first 
line of Fig.n The first correction Tq arises from perturbative insertions of subleading contact interactions 
(i.e. the ones with two derivatives and tx M^) and one-pion exchange (IPE) dressed to all orders by the 
leading contact interactions, see Fig. 01 The undressed static IPE contribution corresponding to the 
second graph in Fig. |21 and based on the Lagrangian in Eq. i)2.24() has the form 

^iTT, undressed 1 ( 9a\^ (<7l • g ) ((72 " 9 ) ,„^ 

where q = p' — pis the nucleon momentum transfer and cjj (r^) are spin (isospin) matrices of the nucleon 
i. The overall normalization of the T-matrix is consistent with Eq. I|2.29jl . It is clear from Eq. (|2.4U|) 

*^Resumming the effective range corrections to all orders in the dimeron propagator leads to an unphysical pole which 
might cause problems in the solution of the 3N scattering equation 0]. 
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Figure 3: The first correction to the NN scattering amplitude in the KSW approach. Dashed hues refer 
to pions, solid rectangles denote insertions of vertices with two derivatives or proportional to M^. For 
remaining notation see Fig. ^ 



that this IPE contribution as well as the contributions from the last two graphs in Fig. are of the 
order 0{k^). Notice that the coefficients of contact interactions with 2m derivatives and oc M^^~'^"^ are 
assumed to scale as 1/A;"^^. Two-pion exchange (2PE) is suppressed compared to IFF and starts to 
contribute at N^LO. At each order in the perturbative expansion, the amplitude is made independent on 
the renormalization scale by an appropriate running of the LECs Ci, Di. As a nice feature, the KSW 
approach allows to derive analytic expressions for the scattering amplitude. In order to conclude on the 
usefulness of the KSW expansion with perturbative pions, it is crucial to understand the scale at which 
it fails. While in the single-nucleon sector, this scale is associated with the chiral symmetry breaking 
scale, ~ Mp ~ 4:TtF-^ ~ 1 GeV, the chiral expansion in the NN sector was found to entail the new 
scale A^r^v associated with the iterated IFF contributions. Estimations based on dimensional analysis 
yield Ann = 167rF^/{g\m) ~ 300 MeV EES- 

In IQRH . an even more conservative result was obtained: 
A-NN = 47rF^/((7^m) ~ 70 MeV. The small estimated values of Ann already indicate that the expansion 
based on perturbative pions might converge poorly. Clearly, dimensional analysis only provides a fairly 
rough estimation for the scale Ann- The convergence of the KSW expansion can ultimately be only 
tested in concrete calculations. The 2N system has been analyzed at N^LO in |lOll |. While the results 
for the ^Sq and some other partial waves including spin-singlet channels were found to be in reasonable 
agreement with the Nijmegen partial-wave analysis (PWA), large corrections show up in spin-triplet 
channels already at momenta ~ 100 MeV and lead to strong disagreements with the data. This is 
exemplified in Fig. 0] The perturbative inclusion of the pion-exchange contributions does not allow to 
increase the region of validity of the EFT compared t o the pionless theory. The failure of the KSW 
approach in the spin-triplet channels was associated in | lQl | with the iteration of the tensor part of the 
IPE potential. Further evidence o f the poo r convergence of the KSW expansion with perturbative pions 
was given by Cohen and Hansen |l02l . Il03| | who obtained predictions for the effective range and shape 
coefficients in the effective range expansion at lowest nontrivial order. These coefficients are sensitive to 
pion dynamics and were found to be poorly described, which indicates that the chiral expansion is not 
converging. More details on the KSW approach with explicit pions and its applications in the two- and 
three-nucleon sectors can be found in 9 4l. 9 f) | and referen ces therein. For further discussion on the role 
of the pion-exchange contributions see Il04ni0,ill0^ . ll07| . 



A suitable way of including the pion-exchange contributions nonperturbatively was proposed in the 
seminal work of Weinberg [j, Q| which preceded the development of the KSW approach and caused a 
flurry of activities to apply EFT in the few-nucleon sector. Weinberg's original arguments are formulated 
in terms of "old-fashioned" time-ordered perturbation theory, see e.g. which is an appropriate tool 

since we are dealing with nonrelativistic nucleons. Consider the S-matrix for few-nucleon scattering 

Saf3 = S{a -13)- 27ri6iE^ - Ep)T^p , (2.41) 

where a and (3 denote the final and initial few-nucleon states and E'q,, En are the corresponding energies. 
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Figure 4: The NN phase shift versus CMS momentum p. The solid Hne is the Nijmegen multi- 
energy fit lOSl . 109( 1 ■ the long dashed line is the LO EFT result, the short dashed line is the NLO result, 
and the dotted line is the N^LO result. The dash-dotted line shows the result of including a higher-order 
contact interaction. Figure courtesy of Ian W. Stewart. 



The T-matrix can be evaluated in "old-fashioned" time-ordered perturbation theory via 

^ rzT \ , Sr^ iHj)aa{Hl)al3 , {Hl)aaiHi)ab{Hl)bl3 , . 

Iai3 = [^i)ai3 + ^ — + / . TE' f i ■ TP > ■ \ + • • • ' (-•^•4^) 



a 



where Hj is the interaction Hamiltonian corresponding to the effective Lagrangian for pions and nucle- 
ons.*^ Here, we use Latin letters for intermediate states, which, in general, may contain any number of 
pions, in order to distinguish them from purely nucleonic states denoted by Greek letters. We remind the 
reader that no nucleon-antinucleon pairs can be created or destroyed due to the nonrelativistic treatment 
of the nucleons. Consequently, all states contain the same number of nucleons. It is useful to represent 
various contributions to the scattering amplitude in terms of time-ordered diagrams. For example, the 
Feynman box diagram for NN scattering via 27r-exchange can be expressed as a sum of six time-ordered 
graphs, see Fig. [51 which correspond to the following term in Eq. ()2.42() : 

" ^ {Ep -Ea + ie){Ep -Eb + ie){Ep - E, + ie) ' ^ ' > 



abc 



where H-j^jsji^ denotes the ttNN vertex. It is easy to see that the contributions of diagrams (d-g) are 
enhanced due to the presence of the small (of the order Q'^/m) energy denominator associated with the 

*^Notice that in contrast to the purely quantum mechanical consideration in section lT^ one has now to account for nucleon 
self-energies. This can be achieved by a proper separation between the unperturb ed Hamiltonia n and the interaction and 
using the formulation in terms of the corresponding "in" and "out" states, see e.g. |llfllllll . lll2| . 
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(a) 



(b) (c) (d) (e) (f) (g) 



Figure 5: Two-pion exchange: Feynman diagram (a) and the corresponding time-ordered graphs (b-g). 
Solid (dashed) hnes correspond to nucleons (pions). 



purely nucleonic intermediate state |6) which in the CMS takes the form:*^ 

1 _ 1 

E/3-Eb + ie p^/m- p^/m + ie' 



(2.44) 



Notice that the energy denominators corresponding to the ttNN states \a) and |c) are of the order M,r ~ 
Q. According to Weinberg, the failure of perturbation theory in the few-nucleon sector is caused by the 
enhanced contribution of reducible diagrams, i.e. those ones which contain purely nucleonic intermediate 
states. To see how this difficulty can be dealt with, it is useful to rearrange the expansion in Eq. (|2.42)1 
and to write it in the form of the Lippmann-Schwinger equation 

T„, = (Kfr).,+E#%T7;' ('-^'^ 

^ hfj-E^ + ie 

with the effective potential {Ves)ai3 defined as a sum of all possible irreducible diagrams (i.e. the ones 
which do not contain purely nucleonic intermediate states): 

^ Ep - Ea + te ^ [E^ - Ea + ie){Ep - E-^ + le) 

ab 

Here, the states |a), |6) contain at least one pion. The effective potential in Eq. 1)2. 46(1 does not contain 
small energy denominators and can be obtained within the low-momentum expansion following the usual 
procedure of CHPT. The contribution of a given irreducible time-ordered diagram can be shown to be of 
the order {Q/A)'^ 0)Q| with A being the scale which enters the values of the renormalized LECs, where 

z/ = -2 + 2Ar + 2(L-C) + ^yiAi, where Aj = + - 2 . (2.47) 

i 

Here N , L, C and Vi are the numbers of nucleons, loops, separately connected pieces and vertices of type 
i, respectively. The quantity Aj gives the chiral dimension of a vertex of type i. Further, di denotes 
the number of derivatives or M^r insertions and n,- is the number of nucleon lines at the vertex i. Notice 
that eq. JUB is modified compared to the one given in Refs. % 0, ^ in order to account for the 
proper normalization of the A^-nucleon states. Chiral symmetry guarantees that Aj > 0. Consequently, 
the chiral order u is bounded from below and for any given v only a finite number of diagrams needs 
to be taken into account. Notice that Eq. (|2.47j) supports a rather natural view of nuclear dynamics, in 
which nucleons interact mainly via 2N forces while many-body forces provide small corrections. After the 

**Equivalently, evaluation of the Feynman graph (a) in Fig. lousing the standard heavy-nucleon propagator of the form 
i/{p'^ + ie) leads to infrared divergences resulting from a pinch singularity associated with the poles — ±ie. These infrared 
divergences are avoided (but still leading to the enhancement in the amplitude) by the inclusion of the kinetic energy term 
in the heavy-nucleon propagators. 
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potential is obtained at a given order in the chiral expansion, few-nucleon observables can be calculated 
by solving the Lippmann-Schwinger equation (|2.45j) . which leads to a nonperturbative resummation of 
the contributions resulting from reducible diagrams. It is easy to see from Eq. ()2.47|) that the leading- 
order (u = 0) potential results from contact interactions without derivatives and the Ivr-exchange. This 
has to be contrasted with the KSW approach, where the exchange of pions is suppressed compared to the 
lowest-order contact terms. It should be understood that the power counting rules in Eq. 1)2. 47() apply 
to renormalized matrix elements. After removing the ultraviolet divergences by a redefinition of the 
LECs in the effective Lagrangian, the remaining integrals are effectively cut off at momenta of the order 
of the soft scale Q. The power counting described above is based on an assumption, sometimes referred 
to as the naturalness assumption, that a renormalized coupling constant C of dimension [mass]"" can be 
written in terms of a dimensionless coefficient c ~ 0{1) as C = cA"". Clearly, higher-dimensional 
terms in the amplitude are only suppressed if the hard scale A that enters the values of the LECs is 
sufficiently large, i.e. if Q ~ <^ A. The validity of the naturalness assumption can, at present, only 
be verified upon performing actual calculations. 

The presence of shallow bound states in few-nucleon systems suggests that the perturbative (iterative) 
solution of Eq. (|2.45|) does not converge. As pointed out by Weinberg |^, Q], this requires for the nucleon 
mass to be counted as a much larger scale compared to the hard scale A. To see that consider the iteration 
of the leading order potential V^^ in the Lippmann-Schwinger equation (|2.45|) which, in operator form, 
can be written symbolically as 

T = + Go V}^ + VS' Go V}^ Go V^^ + . . . , (2.48) 

where Go is the free 2N resolvent operator. One can estimate the size of the leading-order potential by 
the size of the static Ivr-exchange potential leading to V^^^ ~ Since each momentum integration 

in Eq. ()2.48|) gives an additional factor (5^/(47r)^ and Go ~ m/Q'^, one finds that the (n + l)-th term in 
the above equation is suppressed compared to the first term by (Qm/A^)", where we used the estimation 
A ~ AttFtj. The requirement that all terms in the right-hand side of Eq. ()2.48l) are of the same order in 
order, which justifies the necessity of the nonperturbative treatment and enables to describe the physics 
associated with the low-lying bound states, therefore leads to the following counting rule for the nucleon 
mass |3,i,Eil: 

-~T^- (2.49) 
m A2 

This counting rule will be adopted in the present work. Clearly, this estimation based on the naive 
dimensional analysis is fairly crude. A somewhat different estimation can be found in j^H]- Notice that it 
is hardly possible in such an estimation to keep track of various numerical factors, even of the large factors 
such as 47r. For example, the leading and subleading 27r-exchange potentials, both arising from 1-loop 
diagrams, differ by a factor 4tt, see section Hj. 2. 21 Fortunately, the particular way of counting the nucleon 
mass is not crucial from the practical point of view since it only determines the relative importance of the 
relativistic corrections to the nuclear force but does not affect the lowest-order potential and, therefore, 
also not the dominant contribution to the scattering amplitude. Finally, we stress that Weinberg's 
power counting does not explain the unnaturally large values of the NN S-wave scattering lengths or, 
equivalently, the unnaturally small binding energies of the deuteron and the virtual bound state in the ^5*0 
channel. This has to be achieved via an appropriate fine tuning of the lowest-order contact interactions. 



*^We stress that while perturbative renormalization of the scattering amphtude in the pion and single nucleon sectors 
is a straightforward task, both from the conceptual and practical points of view, nonpertu rbative renormalization in the 
few-nucleon sector still attracts the interest of many researchers, see e.g. |ll4llll5Llll6llir^ for some recent work. We will 

address this issue in some detail in section Pi. 1.21 

*^°For LECs accompanying NN contact interaction the expected scaling is C = cF^^A"""*"^, see e.g. |ll8ll . 
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To summarize, the "Weinberg program" for describing the low-energy dynamics of the few-nucleon 
systems proceeds in two basic steps which will be discussed in detail in the next sections of this review. 
First, the few-nucleon potential has to be derived from the effective Lagrangian for pions and nucleons 
using the framework of chiral perturbation theory. Secondly, the corresponding dynamical equations with 
the resulting potential as an input have to be solved. 



3 Nuclear forces in chiral effective field theory 



In the previous section we have introduced the basic concept of the Weinberg approach to few-nucleon 
systems. We will now discuss the structure of the nuclear force in the lowest orders in the chiral expansion 
based on the effective Lagrangian 



^(1) 



1 



1 



4F2 



(tt X tt) 



N 



^CsiN^N){N^N) - ^CTiN^^N){N^aN) + 



C4 

2F2 



(^ijkeabcCTiTaiVj Vr;,)(Vfc VTc) 



^{N'<N){N'<aTN) ■V7r--E {N^N){N'<tN) ■ (A^ViV) + . . . 
AF 2 



N 



(3.1) 



where the superscripts denote the vertex dimension Aj, see Eq. ()2.47l) . q, Cs,t, D and E are LECs and 
ellipses refer to terms with more pion fields. Notice that the nucleon kinetic energy contribute, according 
to Eq. (|2.49j) . to C'^'^\ The above terms determine the nuclear potential up to N^LO (with the exception 
of the NN contact terms at NLO) in the limit of exact isospin symme try. More complete expressions for 

the Lagrangian including higher-order terms can be found e.g. in 0, E3, EH, Eil, Eii, iim, dsa] . 



3.1 Nuclear potentials from field theory 

The derivation of a potential from field theory is an intensively studied problem in nuclear physics. 
Historically, the important conceptual achievements in this field have been done in the fifties of the last 
century in the context of the so called meson field theory. The problem can be formulated in the following 
way: given a field theoretical Lagrangian for interacting mesons and nucleons, how can one reduce the 
(infinite dimensional) equation of motion for mesons and nucleons to an effective Schrodinger equation for 
nucleonic degrees of freedom, which can be solved by standard methods? It goes beyond the scope of this 
work to address the w hole variety of different techniques which have been developed to construct effective 
interactions, see Ref. |l26j | for a comprehensive review. We will now briefly outline a few methods which 
have been used in the context of chiral EFT. 



We begin with the approach developed by Tamm |l27l | and Dancoff |l28l | which in the following will be 



referred to as the Tamm-Dancoff method. Consider the time-independent Schrodinger equation 

{H^ + Hi)\^)=E\^), (3.2) 

where |^') denotes an eigenstate of the Hamiltonian H with the eigenvalue E. One can divide the full 
Fock space in to the nucleonic subspace \(j)) and the complementary one and rewrite the Schrodinger 
equation (|3.2|) as 

■qHr] rjHX 



XHtj XHX 



E , (3.3) 
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where we introduced the projection operators rj and A such that 
state 1^) from the second hne of the matrix equation (|,3.3|) as 



1 



E- XHX 



-H\ 



Al^'). Expressing the 



(3.4) 



and substituting this in to the first hne one obtains the Schrodinger-hke equation for the projected state 

{Ho + V™iE))\^) = m, (3.5) 
with an effective potential Vcs{E) given by 



1 



E-XHX 



XHiri . 



(3.6) 



It is easy to see that the above definition of the effective potential is identical with the one given in 
Eq. (|2.46|) in the context of "old-fashioned" time-ordered perturbation theory. We stress that in order 
to evaluate V^{E) one usually has to rely on perturbation theory. For example, for the Yukawa theory 
with Hi = gHi, the effective potential V^{E) up to the fourth order in the coupling constant g is given 

by 

Ai „ Ai „ A2 „ Ai 



Vl'iE) 



-T] 



-Hi + 0{g' 



(3.7) 



Hq — E Hq — E Hq — E Hq — E 

where the superscripts of A refer to the number of mesons in the corresponding state. It is important 
to realize that the effective potential Vcs{E) in this scheme depends explicitly on the energy, which 
makes it inconvenient for practical applications. In addition, the projected nucleon states \(p) have a 
normalization different from the states |^) we have started from, which are assumed to span a complete 
and orthonormal set in the whole Fock space: 



\HiX 



1 



E-XHX 

Note that the components ipi in this equation do, in general, not vanish. 



XHM^ 



(3.8) 



The abo ve rn entioned deficiencies are naturally avoided in the method of unitary transformation jl29l |. 
see also |l30l |. In this approach, the decoupling of the rj- and A-subspaces of the Fock space is achieved 
via a unitary transformation U 

r]Hr] 



Following Okubo [129|, the unitary operator U can be parametrized as 



U 



7?(l+AtA)-l/2 _^t(l + ^^t)-l/2 
A(l+AtA)-l/2 A(l + ^^t)-l/2 



with the operator A = XAt]. The operator A has to satisfy the decoupling equation 

X{H- [A, H] - AHA) r] = 



(3.9) 



(3.10) 



(3.11) 



in order for the transformed Hamiltonian H to be of block-diagonal form. The effective ry-space potential 
can be expressed in terms of the operator A as: 



f>UT 



rj{H - Ho 



[1 + A^A)-^/^{H + A^H + HA + A^HA){1 + A^A)' 



1/2 



Hn 



(3.12) 
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For the previously considered case of the Yukawa theory, the operator A and the effective potential 
can be obtained within the expansion in powers of the coupling constant g, which leads to: 



UT 



eff 



2 / 

-9 V 



-Hi 



\H,-^—-H, + h. c. 



2 {Hq — Erfi ) 



4 / 

V- 9 V 



A2 



A^ 



Ai 



2 {Hq — Ej^) {Hq — E^) {Hq — En 
Hi 



-H, 



H — Hi 7- Hi fj Hi — Hi 

8 [Ho — Erf') [Hq — Efi){HQ — Eri) 



1 



H 



Ai 



8 {Ho — Erf'){Ho — Efj 



-HifiHi 



Ai 



{Ho — Efi 



-Hi + h. c. 



V + 0{g'). 



(3.13) 



In contrast to given in Eq. H3.7|) . does not depend on the energy E. Another difference to the 
Tamm-Dancoff method is given by the presence of terms with the projection operator fj which give rise 
to purely nucleonic intermediate states. These terms are needed to ensure the proper normalization of 
the few-nucleon states. It should be understood that, in spite of the presence of the purely nucleonic 
intermediate states, such terms are not generated through the iteration of the dynamical equation and are 
thus not reducible in the language of section EHl Since all energy denominators in Eq. (|3.13j) correspond 
to intermediate states with at least one pion, there is no enhancement by large factors of m/Q, which is 
typical for reducible contributions. 

The two methods of deriving effective nuclear potentials are quite general and can, in principle, be ap- 
plied to any field theoretical meson-nucleon Lagrangian. In the weak coupling case, the potential can 
be obtained straightforwardly via the expansion in powers of the corresponding coupling constant(s). 
Generalization to the effective chiral Lagrangian requires the expansion in powers of the coupling con- 
stants to be replaced by the chiral expansion in powers of Q/A. For practical applications, it is helpful 
to use time-ordered diagrams to visualize the contributions to the potential. In "old-fashioned" per- 
turbation theory or, equivalently, the Tamm-Dancoff approach, only irreducible diagrams are allowed. 
Their importance is determined by the power counting in Eq. (|2.47|) and the explicit contributions can 
be found using Eq. 1)3. 6(1 . In the method of unitary transformation one can draw both irreducible and 
reducible graphs, whose importance is still given by Eq. (|2.47|) . Notice that these graphs have a different 
meaning from the time-ordered ones arising in the context of "old-fashioned" perturbation theory and 
will only be used to visualize the topology associated with a given sequence of vertices. The structure 
of the operators contributing to the potential can, in general, not be guessed by looking at a given dia- 
gram and has to be determined by solving the decoupling equation (|3.11|) for the operator A and using 
Eq. H3.12() . This is discussed in detail in Refs. 



where it is also demonstrated how to derive 
the effective potential from the effective chiral Lagrangian at an y giv en order u in the low-momentum 
expansion using the method of unitary transformation, see also |l33| | for a different but closely related 
scheme. The explicit expressions for the operators contributing to the potential in the few lowest orders 
can be found in these r efere nces. For issues related to renormalization within the method of unitary 
transformation see Ref. |l34l |. Another Hamiltonian approach, which is similar to the method of unitary 
transf ormation and isusually referred to as the dressed particle approach, is extensively discussed in 
Refs. [i3i . Bli37l . li3i . 



To illustrate how the above ideas work in practice, let us consider the contribution to the leading 27r- 
exchange potential at order u = 2 arising from diagram (a) in Fig. El The Hamilton operator H^^^ 
describing the ttNN vertex of the lowest possible dimension, Aj = 0, corresponds to the last term in 
Eq. (|2.17j) . In "old-fashioned" perturbation theory, the potential arises from diagrams (b) and (c) in 
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Fig. [5] and can be obtained evaluating the appropriate matrix elements of the operator 



UJ UJl + UJ2 UJ 

where uji = \Jk'^ + denotes the energy of a pion with the momentum k. Notice that at the order 
considered, it is sufficient to treat nucleons as static sources. Explicit evaluation of Eq. ()3.14() yields the 
following result in the CMS [n^: 

^2™ = --Axi I + (3.15) 




4(2F^)4 j (27r)3 ujIuj^ 

+ 4 { ^ + ^^^^ ) {d2 ■ [g X / ])(ai • [qx I]) ) , 

\UJ+ + UJ 

with uj± = \J {q^ lY + 4M^ and q being the nucleon momentum transfer. In the method of unitary 
transformation, the potential is due to diagrams (b-g) in Fig. El and is given by: 

VZ" = + \H(^^^H^^^r^H(''^-H(^^rj+-r]H^'^^-H(^^r^H^''^^H'-'^^rj. (3.16) 
2 UJ 2 UJ uj^ 

The first term on the right-hand side of the above equation, V^™, gives the contribution of the irreducible 
graphs (b) and (c) which, for static nucleons, is the same as in the previously considered case. The 
contribution of reducible diagrams (d-g) is given by the last two terms in the above equation. The 
resulting potential has the form |l39l |: 

= -Ji^. I ^ i^"""--^": (r. . If - iA' + 6(S. ■ [fx tW, • |f x 1])] , (3.17) 



2(2F^)4 7 {27t)^ ujIuj^_{uj+ + UJ-) 

Notice that the isoscalar central and isovector tensor components in Eq. ()3.15|) cancel in the method 
of unitary transformation against the contributions from reducible diagrams. We will discuss the chiral 
27r-exchange potential in more detail in section 113.2.21 

Before closing this section, let us mention two other methods to derive energy-independent potentials 
used in the context of chiral EFT. Historically, energy-independ ent e xpressions for the chiral 27r-exchange 



potential at order u = 2 were first obtained by Friar and Coon Il4 [fl using the method described in jl41|. 



Yet another approach was applied e.g. in Refs. |l43 . ll43 . ll4i . ll45| to study chiral 2tt~ and 37r-exchange 



forces, in which the potential is determined through matching to the 5-matrix. 

Last but not least, one should always keep in mind that, in contrast to the on-shell scattering amplitude, 
nuclear potentials them selve s are not experimentally observable and can always be modified via a unitary 
transformation, see e.g. jl46l | for some explicit examples. This non-uniqueness of the nuclear forces should, 
of course, not be considered as a conceptual problem. A similar sort of non-uniqueness at the level of 
the Lagrangian is well known in quantum field theory, where one has the freedom to perform nonlinear 
field redefinitions. Notice that unitary transformations will, in general, affect not only few-nucleon forces 
but also the corresponding nuclear current operators. It is, therefore, important to have a consistent (in 
the above mentioned sense) formulation for 2N, 3N, . . ., forces and current operators. Such a consistent 
formulation is provided by the chiral EFT framework. 
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3.2 Two— nucleon force 



The chiral NN force has the general form 

^2N = K + ^cont, (3.18) 

where Kont denotes the short-range terms represented by A^A^ contact interactions and K- corresponds to 
the long-range part associated with the pion-exchange contributions Both and Kont are determined 
within the low-momentum expansion as will be discussed in the next sections. 



3.2.1 Regular izat ion of the pion exchange contributions 

Let us now take a closer look at the pion-exchange contributions. The explicit form of the correspond- 
ing non-polynomial functions of momenta'^^^ depends, to some extent, on the way one regularizes the 
corresponding loop integrals. Consider, for example, the isoscalar central part of the 2PE potential at 
order v = 3 which results from the triangle diagrams and is given by 

where g = |^|, / = |/ 1 and Cj are the corresponding LECs. The integral is cubically divergent and needs 
to be regularized. Applying dimensional regularization one finds: 

Vc{q) = (2M2(2ci - cs) - esq') {2M^ + q') A{q) , A{q) = 1 arctan ^ . (3.20) 

Here, we do not show polynomial terms of the kind a + /3g^ which contribute to V^ont- It is instructive 
to express the potential using the spectral function representation: 



Fc(,) = ^r*4^. (3.21) 

71" J2AU + ^ 



where the spectral function p{fi) can be obtained from Vc{q) in Eq. 1)3.20(1 via 

pifi) = 9 [yc(0+ - ip)] = {2Ml{2ci - cs) + cs/i') {2Ml - p')U{p - 2M^) . (3.22) 

In Eq. (|3.21() . the twice subtracted dispersion integral is given which is needed in order to account for 
the large-// behavior of p{p). Eq. (|3.21|) shows that the 2PE potential resulting from Eq. ()3.19|) upon 
applying DR contains explicitly the short-range contributions associated with the integration over large 
values of p. These short-range contributions to the potential are an artifact of the chosen regularization 
procedure (i.e. DR). They are model dependent and cannot be predicted in the chiral EFT framework 
since the chiral expansion for the spectral function p(ij,) is invalid for large values of p. Instead of keeping 
the spurious short-range physics in the 2PE potential, one can perform the spectral function integral 
only over the low-p region, where chiral EFT is applicable. This can be achieved using the regularized 
spectral function 

p{p)^p''{p)=p{p)e{A-p), (3.23) 

Polynomial contributions to the potential are represented by a series of contact interactions V^cont. 
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with the reasonably chosen finite ultraviolet cut-ofF A which prevents that the regularized 2PE potential 
Vq given by 

= ^-^r di.\ 4t^2 = (2M.^(2ci - C3) - c.q') {2Ml + A%) + ..., (3.24) 

with the ellipses referring to polynomial (in q^) terms has components with the range r < A^^. In the 
above equation, the regularized loop function A^{q) turns out to be 

A^{q) = e{A - 2M^)— arctan ~ ^^"^ . (3.25) 

In wha t follows, we will refer to the above described regularization scheme, which has been introduced 
in 



143], as to spectral function regularization (SFR). 



What is the relation between the DR and SFR potentials in Eqs. (|3.2fl|) and (|3.24|) ? To see that one can 
rewrite the spectral function integral in Eq. H3.21() as follows: 

2q^ r , 1 p(/i) M 1 p(/i) , 2g4 1 p{,x) 



dfi K + aq^ + pq" + ... , (3.26) 



71" J2M^ ' /^^ /^^ + q"^ 
where 

It is, therefore, obvious, that the DR and SFR potentials differ from each other by an infinite series of 
higher-order contact interactions. In fact, they might be viewed as two different conventions to define the 
non-polynomial part of the two- and more-pion exchange potential. DR corresponds to the convention, 
according to which the non-polynomial part includes components of arbitrarily short range which are 
strongly model-dependent. In contrast, the SFR approach uses the convention, according to which only 
the components with the range r > A~^ are explicitly kept in the non-polynomial part of the potential 
while all shorter-range contributions are represented by the contact interactions. In general, for quickly 
converging expansions, both the DR and SFR methods are completely equivalent provided the ultraviolet 
cut-off A is chosen to be large enough. For example, we will see in section l4. 1 . ![ that both schemes lead 
to similar results for peripheral NN scattering at order u = 2. If, however, the convergence for some 
well understood physical reason is slow and (some) observables become sensitive to higher-order counter 
terms, it is safer to avoid the spurious short-distance contributions kept in DR. In such a case, SFR is 
a preferable choice. An example of such a situation will be considered in sectio n 14.1.11 Notice that a 
similar approach based on a finite momentum c ut-of f was used in Ji^, H EE^ to deal with the slow 



convergence in the SU(3) bar yo n CHPT, see also [I5ll | for a recent application to the chiral extrapolation 



of the lattice QCD results and 152l . 1531 1 for a discussion on cut-off schemes in CHPT. 



It is also instructive to compare the DR and SFR potentials in configuration space. For r > 0, the inverse 
Fourier-transform can be expressed in terms of the spectral function p(/x) via 

Vc{r) = 1^ r d/i/^e-^X/x). (3.28) 
27r^r J2M^ 

Substituting p{ii) from Eq. (|3.22l) in to Eq. 1)3. 28() . one obtains the following expression for the potential 
corresponding to DR: 



Vc{r) 
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-2x 



327r2Fi 



2ci x^{l + xf + C3(6 + 12x + lOx^ + Ax^ + 



(3.29) 
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Figure 6: The potential Vc in r-space. The sohd hne (band) shows the DR (SFR, A = 500 . . . 800 MeV) 
result. The short-dashed (long-dashed) line refers to the phenomenological a (a + oj + p) co ntributions 
based on the isospin triplet configuration space version (OBEPR) of the Bonn potential 



where we have introduced x = M-j^r. Using the regularized expression for the spectral function in 
Eq. H3.23() . one obtains for the SFR potential: 



Voir) = Voir) 



^9a 



1287r2Fi 



4cix^ 2 + y(2 + y) -2x^ 



+ C3 (24 + y(24 + I2y + ^y"^ + y^) - 4x^(2 + 2y + y^) + Ax 



(3.30) 



where y = Ar. 



In Fig. El we compare the isoscalar central part of 2PE obtained using DR and SFR for the central 
values of the LECs ci,3, ci = -0.81 GeV"^ and C3 = -4.70 GeV~\ from Ref. 154]. Clearly, the large- 
distance asymptotics of the potential, which is constrained in a nontrivial way by chiral symmetry of 
QCD, is unaffected by the cut-off procedure (provided A » M^). The strongest effects of the cut-off 
are observed at intermediate and shorter distances, where 2PE becomes unphysically attractive if DR is 
used. In contrast, removing the large components in the mass spectrum of the 2PE with the reasonably 
chosen cut-off A = 500 . . . 800 MeV greatly reduces this attraction and yields the potential of the same 
order in magnitude as the one obtained in phenomenological boson-exchange models. We will see in 
section |4. 1 . II how a choice of regularization affects the results for peripheral NN scattering at N^LO. 

We further stress that other regularization schemes may be applied as well. For example, one can 
regularize divergent loop integrals using an ordinary momentum-space cut-off. The prominent feature of 
the SFR scheme is given by the fact that it only affects the two-nucleon contact interactions. One can, 
therefore, directly adopt the values for various LECs resulting from the single-nucleon sector analyses, 
where dimensional regularization has been used. This, in general, is not the case for a finite momentum 
cut-off regularization. 



3.2.2 Pion exchange contributions 



Consider now pion-exchange contributions to the potential 



(3.31) 
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where one-, two- and three-pion exchange (3PE) contributions Vit^, V2-k and Vst^ can be written in the 
low-momentum expansion as 

y^^ = yS+yS+y^ + -. 

1/3. = ^3? + .... (3.32) 

Here, the superscripts denote the corresponding chiral order and the elhpses refer to (Q/A)^- and higher 
order terms. Contributions due to the exchange of four- and more pions are further suppressed: n- 
pion exchange diagrams start to contribute at the order (Q/A)^""^. Notice further that in this section 
we restrict ourselves to isopin-invariant contributions. Isospin-breaking corrections will be discussed in 
section 13.2.51 The corresponding relativistic corrections will be considered in section 13.2.41 

The static IPE potential at N^LO has the form 

..(0) , ..(2) , ..(3) , ..(4) / 9a\^ , ..2 Oi-qa2-q , . 

^iV + ^iV + ^iV + ^iV = - j (1 + ^) ^1-^2 ^2^M^ • ^ ^ 

Here 6 denotes an isospin-conserving Goldberger-Treiman discrepancy 

2d 



-Ml + kM^ , (3.34) 



■'tt ' 

gA 

where dig is a LEG from the dimension three vrA^ Lagrangian and the constant k determines the size of 
the first correction to the Goldberger-Treiman discrepancy. Here and in what follows, the expressions for 
the nuclear force should be understood as operators with respect to spin and isospin quantum numbers 
and matrix elements with respect to momentum variables. We further stress that all one- and two-loop 
Ivr-exchange diagrams at this order lead to renormalization of various LEGs without introducing any 
form-factor-like behavior. The derivation of th e IPE potential to one loop in the method of unitary 
transformation is discussed in detail in Ref. 11." 



We now turn to the 2PE contributions. It is convenient to express V2tt in the GMS in the form: 

V27r = Vc + Ti . T2 Wc + [Vs + Ti ■ T2 Ws] <?! • (72 + [Vt + T"1 • T2 Wt] <71 • q(J2 ' Q (3.35) 

+ [Vls + ri-T2 Wls] «(<?! + ^2) ■ iqxk) + [V^l + ti ■ T2 W^l] <?! • (g x k)a2 -{qxk), 

where the superscripts C, S, T, LS and aL of the scalar functions Vc, ■ ■ ■, Wo-l refer to central, spin- 
spin, tensor, spin-orbit an d quadratic spin-or bit components, respectively. The chiral 2PE potential 
vj^^ + Vg*!^^ is discusse d in Il3ll . Il39l . Il4(]| . Il42l | and in |ll9l | using an energy-dependent formalism. For 
a related work see also The NLO 2PE potential is given by the contributions of the box, crossed- 

box, triangle and football diagrams shown in the first line of Fig. [7| which in the energy-independent 
formulation read 

W!,'\q) = -,,^l^LHg)\4Mli5g\-Agl-l) + qH23g\-10gl-l)+'^^ 



3847r2F4 { / . V . . 4M2 + g2 

y?\i) = -^2yPi^) = -^.LH,). (3.36) 

Here, the loop function L^{q) is given by 

[2, ,2 I „2„2 



L%) = e{A - 2M.) f In ''^,^1-^^'''" > - = V^n^^ , . = - ■ (3-37) 
2q 4M2(A2 + g2) v 
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Figure 7: Leading, subleading and sub-subleading contributions to the chiral 27r-exchange potential. 
Solid (dashed) lines correspond to nucleons (pions). Solid dots, filled rectangles and filled diamonds 
represent vertices with Aj = 0, 1 and 2, respectively. Shaded blob denotes the next-to-next-to-leading 
order contribution to the pion-nucleon scattering amplitude. 



These expressions are based on the SFR approach. The corresponding DR expressions can be obtained 
taking the limit A — > oo. We further emphasize that a significant part of the NLO 2PE contributions was 
considered much earlier in the context of meson theory of nuclear forces, see e.g. 

N^LO one has to take into account the contributions of the triangle and football graphs in the second 
line of Fig. [7|with a single insertion of the subleading ttttNN vertex. One finds: 



V!,'\q) = -^{2M,^(2ci-C3)-C3g2}(2M,2 + ,V^(g), 
W^'\q) = -^wf\q) = -^^^c,{AMl+q^)AHq). (3.38) 
where the N^LO loop function A^[q) has been defined in Eq. ()3.25|1 . 

The N^LO corrections to the 2PE potential Vg^^^ have been recently calculated by Kaiser |l45l | and are 
schematically depicted in the third line of Fig. [3 They arise from two groups of diagrams, the one- 
loop football graphs with both dimension two ttitNN vertices of the ci^...^4-type and the diagrams which 
contain the third order pion-nucleon amplitude and lead to one-loop and two-loop graphs. We begin 
with the first group of corrections, for which one finds: 



V^,'\q) = Y^,LHq)\[^u.' + Cs{2M'^+q')-4c,M'^ 



2 a 
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The expressions for the second group of corrections were obtained by Kaiser [I45l | and are given in 
appendix in terms of the corresponding spectral functions. 

Three-pion exchange starts to contribute at N^LO and is given by diagrams shown in Fig. |S1 The 
corresponding expressions for the spectr al function s and the pot ential (obtained using dimensional regu- 
larization) have been given by Kaiser in see also [16li ] for a related work. It has been pointed 
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Figure 8: Leading contributions to the 37r-exchange potential. For notation see Fig. d 



out in these references that the resulting 3PE potential is much weaker than the N^LO 2PE contribu- 
tions at physically interesting distances r > 1 fm. Having the explicit expressions for the 3PE spectral 
functions, it is easy to calculate the potential in the SFR scheme. It is obvious, even without performing 
the explicit calculations, that the finite-range part of the 3PE potential in the SFR scheme is strongly 
suppressed at intermediate and short distances compared to the result obtained using DR. This is because 
the short range components which dominate the 3PE spectrum are explicitly excluded in this approach. 
This is exemplified in Fig. IHl for the ca se of the N^LO isoscalar spin-spin contribution proportional to 
g\, which has been found in |l43l . Il44l | to yield the strongest 3PE potential for 0.6 fm < r < 1.4 fm. 
For r > 0.5 fm, it reaches at most 2% — 8% (depending on the choice of the spectral function cut-off) 
of the correspondi ng N^LQ 2PE contribution [l62]. Similar results have been found for other leading 
3PE contributions 
N^LO analyses ^ 

3PE diag rams proportional to LECs q and obtained using DR was found to be sizable at intermediate 
distances Il64l l . The strength of the 2PE and 3PE contribution s considered above and t he corresponding 



143l. 144l | . For that reason the 3PE contributions have been neglected in the present 
" l63l |. Notice, however, that the potential resulting from subleading (i.e. N^LO) 



expressions in coordinate space are discussed in detail in Refs. HHHEii, 



3.2.3 Contact terms 



The short-range part of the potential is represented by a series of contact interactions 



I'cont — "cont "cont "cont ) 



.(2) 
' cont 



.(4) 
' cont 



(3.40) 
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Figure 9: The ratio of the isoscalar spin-spin 3PE and 2PE N^LO contributions using dimensional and 
spectral function regularization. The cut-off in the spectral function varies in the range A = 500 . . . 700 
MeV. 



where the superscripts denote the corresponding chiral order as defined in Eq. (|2.47|) . These terms feed 
into the matrix-elements of the two S-, P- and D-waves and the two lowest transition potentials in the 
following way: 

{P\Vcont\P) = Cppp' + Dppp'{p^+p'^) , 

{D\V,ont\D) = DdP^p'^ , 
{^Si\V,ont\'Di) = C3di-3SiP^+dIoi^ssiP^p'^ + dIdi-3Sip\ 

{^P2\VcontfF2) = D3F2-3P2/P, 

(^i^2|^cont|^P2) = D3F2-3P2PP''^ , (3.41) 

where p = \p\, p' = \p'\ and the subscripts 5 = {150, 3S1}, P = {IPl, 3P0, 3P1, 3P2}, D = 
{1-D2, SDl, 3D2, 3L>3} refer to the corresponding channel. The relations between the spe ctroscopic 
LECs in the above equations and the ones that occur in the Lagrangian can be found in jl62| . Isospin- 
breaking short-range corrections will be specified in section Hi. 2. 51 



3.2.4 Relativistic corrections 



The first relativistic corrections to the nuclear force appear at order v = A provided that the nucleon mass 
is counted according to Eq. (|2.49|) . They result from both 1/m^-corrections to the static IPE potential 
and 1/m-corrections to the order u = 2 static 2PE potential. In addition, at this order one also needs 
to correct the nonrelativistic expression for the nucleon kinetic energy 



i^Hn = V^?^-m~|^-^, (3.42) 

which enters the corresponding dynamical equation. Equivalently, one can use the full, not expanded 
expression for the nucleon kinetic energy which leads to the following Schrodinger equation for two 
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nucleons in the CMS: 



2m] +V 



(3.43) 



where m = nip for the pp, m = ran for the nn and m = 2mpmn/ {rap + rUn) for the np case, see jl62[ | for 
more details on the kinematics. Notice that Eq. (|3.43|) can be cast into equivalent no nrelativis tic forms, 
i.e. into the Schrodinger equations with the nucleon kinetic energy 

advantage of Eq. 1)3. 43() versus the corresponding nonrelativistic forms is that it can easily be generalized 
to the case of three and more nucleons. Changing the form of the Schrodinger equation causes changes 
in the relativistic corrections to the nuclear potential. Relativistic corrections to the interaction can, 
therefore, only be defined within a particular framework. For the Schrodinger equation 1)3. 43() . the 
corrections to the leading IPE potential V^'^ take in the NN CMS the form: 



(4) 



Itt 



p2 j^^l'i 

2rr? 



(0) 



Itt 



(3.44) 



where and p' are the NN initial and f inal CMS momenta. This choice of corrections is sometimes called 
the "minimal nonlocality" choice, see |l46| | and references therein. The corresponding 1/m-corrections 
to the 2PE potential read: 



3^1 r2M5 



5127rmF4\ u;2 



9\ 



1287rmF4 1 uj'^ 



AMI + - 9a ( + 



q 

-^Wf[q) 



-^i^(AMl + -q^\A^{q), 



9l 



9\{^-9\) 



2567rmF4 
{2Ml+qM\<l) 
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{AMl+q^)A%). 



{2Mi + q')A^{q) 



(3.45) 



Notice that the above expressions differ from the ones given in |142l | due to the different form of the 
dynamical equation and relativistic corrections to the IPE potential employed in the present work. For 
an extensive discussion of this issue the reader is referred to Ref. |l4d| where the dependence of relativistic 
corrections on certain kinds of unitary transformations is studied and the general expressions for 
corrections to the IPE potential and 1/m-corrections to the leading 2PE potential are obtained. We 
further stress that if the nucleon mass is counted, contrary to Eq. (|2.49|) . as m ~ A, 1/m^-corrections to 
the IPE potential, 1/m^-corrections to the order v = 2 2PE potential and 1/m-corrections to the order 
u = 'i 2PE potential have to be taken into account at NfLO in addition to terms shown in Eqs. (|3.44)1 . 
()3.45|) . The corresponding expressions can be found in |l65l |.^^^ 

differ from each other by the 



^^■^The two forms of the resulting nonrelativistic Schrodinger equation discussed in [l£ 
definition of the nonrelativistic CMS momentum. 

*^^Notice that the expressions given in |l65(| probably need to be adjusted in order to be made consistent with Eqs. 13. 43 II . 
llOIl and lin^ . 
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3.2.5 Isospin breaking corrections 



Within the Standard Model, isospin violation has its origin in the different masses of the up and down 
quarks and the electromagnetic interactions. Consider first isospin-breaking in the strong interactions. 
The QCD quark mass term can be expressed in the two-flavor case as: 

-^Sass = -^^("iu + ?n-d)(l + eT3)g , where e = , ^'^ ^ ■ (3.46) 

The above numerical estimation corresponds to the light quark mass values based on a modified MS 
subtraction scheme at a renormalization scale of 1 GeV The isovector term (oc t^) in Eq. ()3.46|) breaks 
isospin symmetry and generates a series of isospin-breaking effective i ntera ctions oc [eM^y^ with n > 1. It 
is, therefore, natural to count strong isospin violation in terms of eM^ |l2ll ]. Electrom agnet ic terms in the 
effective Lagrangian can be generated using the method of external sources, see e.g. jl67l . llGS . 169, 17o| 
for more details. All such terms are proportional to the nucleon charge matrix Qch = e (l + r3)/2, where e 
denotes the electric charge. More precisely, the vertices which contain (do not contain) the photon fields 
are proportional to Q^-^ (Qch)' where n = 1, 2, . . .. For processes in the absence of external fields, in which 
no photon can leave a Feynman diagram, it is convenient to introduce the small parameter ~ 1/10 
for isospin-violating effects caused by the electroma gneti c interact i ons. Isospin-violat ing t erms in the 
effective Lagrangian at lowest orders can be found in |i2il . B E Il24| |. see also Ref. |l25l |. 



Isospin-breaking nuclear forces can, in principle, be derived in the EFT framework performing indepen- 
dent expansions in Q/A, e and e. It is, however, convenient to relat e these small parameters with each 
other in order to have a single expansion parameter. In Ref. |l25l ]. the following counting rules were 
adopted: 

Q , , 

The power counting expression in Eq. (|2.47|) can be easily extended to include the contributions due to 
n-y virtual photons: 

1/ = -2 + 2ny + 2iV + 2(L - C) + ^ V^Ai . (3.48) 

i 

Here, the Q-power Aj of the vertex i defined in Eq. ()2.47|) has to be adjusted according to the rules given 
in Eq. (|3.47j) . It should be understood that the counting rules on Eq. (|3.47|) are by no means unique 
and represent an attempt to relate the sizes of the isospin-breaking and isospin-conserving nuclear forces 
with each other in a reali s tic w ay. Different rules are usually adopted in the meson and single-nucleon 
sectors, see e.g. |l22l . Il23l . Il24i |. Counting rules very similar to t he ones in Eg. (13.471) fbut no t exactly 
the same) have been used in applications in the 2N sector in Refs. |l2ll . llTll . Il72 . Il73l . Il74l . llT^ . Clearly, 



changing the counting rules shifts various nuclear force contributions to different orders but does not 
affect their explicit form. The most realistic set of counting rules, i.e. the one that finally leads to the 
most natural values of the LECs, can only be figured out in practical calculations. 

The 2N forces fall in to four classes with respect to their isospin structure ^ 

Class I: Vi = ai + (3i ri ■ T2 , 

Class II: Vu = au rf t| , 

Class III: ^111 = am (rf + t|) , 

Class IV: Viy = aiv {rf - t|) + Av [ti x , 

where aj. Pi are space and spin operators. The operator Piy has to be odd under a time reversal 
transformation. While class (I) forces are isospin-invariant, all other classes (II), (III) and (IV) forces are 
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Figure 10: One-loop tt^ diagrams contributing to the isospin-breaking NN force of the IPE range. Wavy 
hnes represent photons. Sohd dots refer to the leading strong and electromagnetic vertices. Only one 
representative topology is depicted for each graph. For the remaining notation see Fig. EJ 



isospin-breaking. Class (II) forces, Vu, maintain charge symmetry but break charge independence. They 
are usually referred to as charge independence breaking (GIB) forces. Class (III) forces break charge 
symmetry but do not lead to isospin mixing in the 2N system. Finally, class (IV) forces break charge 
symmetry and cause isospin mixing in the 2N system. 

We will now discuss various contributions to the isospin-violating 2N force which have been extensively 
studied in the chiral EFT framework and worked out up to order v = b. It can be expressed as: 

^2N = ^EM + + VL + Vl + Vi,, , (3.50) 

where the terms in the right-hand side refer to the long-range electromagnetic force, pion-photon ex- 
change, isospin-breaking one- and two-pion exchange potentials and contact terms, respectively. The 
superscript I is used in order to distinguish the isospin-breaking from the corresponding isospin-invariant 
contributions considered in sections 13.2.21 13.2.31 Let us first comment on the long-range electromagnetic 
force whose dominant contribution is given by the static Coulomb interaction at order v = 2. The first 
long-range corrections are suppressed by (relativistic corrections to the static one-photon exchange). 
At this order, the long-range electromagnetic NN interaction is given by 

Vem{pp) = Vbi + Vb2 + Vvp + 14dM(pp) , 
VEuinp) = VuMinp) , 

VEMinn) = VuMinn) , (3.51) 

where Vci and Vc2 are usually referred to as "improved Coulomb potentia l". T hev include the relativistic 
1/m^-corrections to the static Coulomb potential worked out in Ref. 177 1. The expressio ns for the 



vacu um p olarization potential Vyp and magnetic moment interaction Vmm can be found in Refs. |178l . ll7£ 
and [iSd], respectively. Notice that ^em contains classes (II), (III) and (IV) forces. The class (IV) force 
is given exclusively by the magnetic moment interaction. We also stress that the effects of Vem are 
enhanced at low energy due to the long-range nature of this force. Even the effects due to Vmm, which 
is suppressed by factor ~ (Q/m)'^ compared to the static Coulomb interaction and thus contributes at 
orde r u = 6, might be large for certain scattering observables under specific kinematical conditions, see 
e.g. US- 

The one-loop diagrams contributing to the isos pin- breaking NN force of the IPE range are shown in 
Fig. Eland were considered by van Kolck et al. jl8ll |. Some of the graphs depicted in this figure lead to 
renormalization of the isospin-breaking IPE potential (i.e. change its strength). Notice further that due 
to isospin, only charged pion exchange can contribute to the nj potential 14-7 ^^'^ thus it only affects 
the np system. The resulting ttj potential is CIB and has the form 

VTr^iq ) = - ^pl^2 ("^1 • '^2 - rf t|) <Ti ■qa2-q V^^{(3) , 
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(3.52) 



Here, /? = \q\/M^+ and 7 is a regularization scheme dependent constant. The analytical form of is 
similar to the one of the IPE potential but differs in strength by the factor a/vr ~ 1/400. 

Isospin-breaking co rrections to the IPE potential were extensively studied within the EFT framework 
ml . B \vrk |l75|. It is convenient to express the static IPE potential in Eq. ()3.33|) in a more general 
form, which already incorporates some (but not all) of the isospin-breaking corrections: 

Vi. + Vl{pp) = {l + 5pfV{M^o), 



-(1 + 5p){l + 5„) y(M,o) + (-1)^+1 2(1 + 5,f V{M^i 



where we utilize the notation of Ref. jl25| . In Eq. (|3.53() . T 
two-nucleon system and V{M^) is defined as: 

in x 2 

V{Mi 



(3.53) 



0, 1 denotes the total isospin of the 



9A 



2F^ 



((Ti • g)((T2 • q, 



(3.54) 



The constants 5p, 5n and 6c in Eq. (|3.53j) specify the isospin dependence of the pion-nucleon coupling 
constant. For the isospin-symmetric IPE potential, the quantity Sp = 6n = 5c gives an isospin-conserving 
Goldberger-Treiman discrepancy.*^^ Charge-symmetry conservation implies dp = 6n Sc- The leading 
isospin-breaking contribution to the pion-nucleon coupling constants is of the strong origin and leads 
to the class (III) force at order 1/ = 3. Up to order v = 5, isospin- violating contributions to 6i^s result 
from various tree- and one-loop diagrams and lead to both classes (II) and (III) forces, see |l25l . Il72l | 
for more details. Apart from the corrections due to Sp ^ 5n Sc, the expressions for the IPE potential 
in Eq. 1)3. 53() also incorporate the order v = 2 class (II) contribution due to the mass difference of the 
exchanged pions, 6M^ = — M^o = (36 MeV)^. This is, in fact, the dominant contribution to the 
isospin-breaking nuclear force. Clearly, the corrections oc (SM^)'^ at order u = 4 are also included in 
Eq. (|3.53|) . Further contributions to the isospin-breaking IPE potential not included in Eq. (|3.53jl arise 
at order = 4 du e to the proton-to-neutron mass difference, 5m = nip — m-n = —1.29 MeV, and are 
given by 



6m ( QA 



2m \ 2K 



1 



Ti X T2 



(fi' + M2) 



((71 • qi )((72 • qi 



P2' 



- ((<?! • Ql){^2 ■ ip2 +P2')) + (o^l • {pi +Pl'))(<?2 • qi\ 



(3.55) 



where pi {pi ') denotes the incoming (outgoing) momentum of the nucleon i and qi = pi —pi = —{P2—P2)- 
Notice that the first term in the square bracket vanishes in the CMS. The potential in Eq. ()3.55|) represents 
the dominant clas s (IV) force. The dependence of Vj^^^^ on the total momentum of the 2N system is 
explained in |l75l |. In addition to the corrections linear in 6m/ m in Eq. ()3.55|) . one has to take into 
account the CIB contribution oc {6m)'^ which reads 121 . 125l . 175l . ll8ol |: 



V/J^^ = -{6mf 



9A_ 
2F^ 



7^1 ^2] 



((Tl • q){a2 ■ q] 



(3.56) 



(g2+M2)2 • 

Notice that the isospin-violating piece has the same structure as the correction due to the pion mass 
difference at order u = 2 but is 6M^/(6m)'^ ~ 660 times weaker. No new structures in the IPE potential 
appear at order = 5. 

Notice that in addition to terms (X (m„ + rrid) in Eq. II3.34L one has now to include corrections cx |m„ — md\ and oc a. 
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Figure 11: Leading isospin-breaking corrections to the 2PE potential. A light-shaded circle inserted at a 
pion (nucleon) line refers to a single insertion of the pion (nucleon) mass difference. Crossed circles denote 
the leading isospin-breaking -mrNN vertices of dimension Aj = 2. For graphs with pion/nucleon mass 
difference insertions, only one representative topology is depicted. For remaining notation see Fig. d 



Let us now switch to the 2PE potential. The dominant contributions arise at order i/ = 4 from 1-loop 
diagrams constructed from the leading ttNN and ttttNN vertices and a single insertion of the pion mass 
difference, the proton-to-neutron mass shift or the order Aj = 2 isospin-breaking ttttNN vertex. The 
corresponding diagrams are shown in Fig. 111! Let us first disc uss t he class (II) contributions arising 
from diagrams with a single insertion of SM^. As shown in Ref. |l73| |. the potential can be expressed in 
terms of the corresponding isospin-invariant contributions given in Eq. 1)3. 36() without performing any 
additional calculations. To that aim, one can first decompose the isospin-invariant 2PE potential in to 
the isoscalar and isovector pieces 

V2n = Vl + VlTi-T2 . (3.57) 

The leading isospin-breaking effects due to M^± ^ M^o are incorporated properly if one uses M^r, defined 
as 

= ^M,± + Im^o , (3.58) 
in the scalar part V^l adopts for the isovector part: 

V2^^(Mj^±) for pp and nn , 

VliM^o) for np, T = l, (3.59) 
Vi^{Mn) for np, T = . 

These results are valid modulo (JM^/M^) ^-corrections which first contribute to the 2PE potential at 
order = 6. Clearly, one can also express the corresponding potential in a more t radit ional form in terms 
of isospin matrices and without referring to particular isospin channels, see Ref. jl25l |. 



All remaining 2PE contributions at order u = A lead to c l ass (III) f orces and were derived independently 
and using 
One finds 



and using different methods by several groups, see |l25l . 1 1741 . Il82l | and |l83l | for a related earlier work 



to c l ass (III) f orces and v 
EH, Em Em and |l8^ 
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2g\5mMl 



Ag\ 5m - {5mr' (2M^ + q') A^{q) 



4:M^ + 

+ (^1 + ^1) [(^1 • • q) - [ai • as) f] A^q) . (3.60) 
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Figure 12: Subleading isospin-breaking corrections to the 2PE potential. Crossed rectangles denote 
isospin-breaking vertices of dimension Aj = 3. Graphs resulting from the interchange of the vertex 
ordering are not shown. For remaining notation see Figs. 171 ITTl 



Here the LEG accompanying the order Aj = 2 isospin-breaking tttiNN vertex is expressed in terms of 
to the strong contribution [SmY^^ to the nucleon mass shift dm = {5mY^'^ + (^m)*^™, where 

{5mf^' = {nip - ninT' = -2.05 ± 0.3 MeV , 

{Smf^ = {nip - mnf^ = 0.76 ± 0.3 MeV . (3.61) 

These values are taken from Ref. jl84l | and based on an evaluation of the Gottingham sum rule. Notice 
that the vertices corresponding to {SmY^^ and {Smy^ have, according to the counting rules in Eq. (|3.47jl . 
dimensions Aj = 2 and 3, respectively. Notice further that some of the diagrams shown in Fig. 1111 in 
particular the planar box and the football diagrams with a single insertion of dm, lead to vanishing 
contributions. 

Gonsider now the subleading isospin-breaking 2PE potential which is due to 1-loop diagrams (some of 
which lead to vanishing results) shown in Fig. 1121 The class (II) contributions oc can be obtained 
using Eq. (|3.59|) from the corresponding isospin-invariant 2PE potential in E q. (|3.38() . The remaining 
terms were calculated recently using the method of unitary transformation |l2fj | . The resulting class (III) 
potential reads: 

Kl(^) ^ --^{rl^rMi-ol^ni'^m^ 



2- - 96^^'i ^'2^^ ^'^^l ^^"'^ 4M2 + g2 

g\ bm (18ci + 2c2 - 3c3) + (ibm - {bmf^^ ) (6ci - C2 - 3c3) 



+ 4M2 



g\ bni (5c2 - I8C3) - (25m - {bmf^-"^ (ca + 6C3)] | 



2 

%4 (^1 + ^1) [(<?! • 9 )(<?2 • q) - (<?i • S2) f'] L\q) ibm C4 + <7A /?) , (3.62) 



167r2F4 

where /3 is a combination of the LECs accompanying the leading isospin-breaking ttNN vertex. At order 
= 5, it can be expressed in terms of the constants bi defined in Eq. 1)3. 53(1 as /5 = (1/A)( bp — br,)qA ■ The 
complete expressions for the isospin-breaking 2PE potential in r-space can be found in [l25|. 
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Finally, isospin-breaking contact terms up to order v 
P-waves in the following way: 



5 feed in to the matrix-elements of the S- and 



(^5*0, PPlKontl^'S'o, VV) 
(^5o, nnlVcontl^S'o, nn) 
ePo, pp\ViJ'Po, pp) 
fPi, pp\ViJ^Pi, pp] 

(3P2, pp\Vi^tfP2, pp] 
{^Pi,np\Vij^Pi,np] 



Prso + Piso{p'+p") 



J2\ 



Piso-f^isoip'+p'- 
-ePo, nn\ViJ^Po, nn] 
-(M, nn\ViJ'Pi, nn] 
-(3P2, nn\Vi^,\^P2, nn] 
/5iPi-3Pi pp' , 



Pspo pp' 
P3P1 pp' 

P3P2 PP 



(3.63) 



where f3i, Pi are the corresponding LECs. Here, we use the convention according to which the np matrix 
elements (with exception of the last term in Eq. H3.63() ') do not change by switching off isospin-violating 
contact terms. Notice further that all terms quadratic in momenta are of the order u = 

Last but not least, we would like to emphasize that chiral EFT supports the following hierarchy of the 



nuclear forces 



12l|: class (I) > class (II) > class (III) > class (IV). 



3.2.6 Towards a relativistic NN potential 

In this section, we will discuss the problem with the formal inconsistency of the heavy-baryon (HB) 
expansion mentioned in section 12.11 and its impact on the nuclear potential and NN scattering. For 
that, consider the contribution to the NN scattering amplitude arising from the triangle diagrams in the 
second line of Fig. [Tj We ass ume exact isospin symmetry throughout this section. Further, we follow 
closely the procedure of Ref. |l85l | and restrict ourselves to the central part of this contribution, which 
is proportional to the LECs ci^s and can be expressed in terms of the vrA^ scattering amplitude and the 
nuclear scalar form factor. For a discussion of the nuclear scalar form factor in relativistic CHPT, the 
reader is referred to Ref. [13], see also |l58l] for the results in the infrared regularized version of C HPT. 
Using the relativistic expression for the nucleon propagator and utilizing the notation of Ref. |l85| , this 
contribution to the scattering amplitude can be written as: 

TF'it) = (2M2(2ci - C3) + c,t) [j,,,(t) - J,%(t)] . (3.64) 

Here, t = , Q = p' — P and the integrals Jc,c{t) and J^^gi^it) read: 

(4vr)2 f dH 1 



i J (27r)4[(^-g/2)2-M2+ie][(/ + g/2)2-M2+ie] 
(1) _ (47r)2 f dH 2mV-l 



iy2 J (2TT)mi-q/2y - M^+ie][{l + q/2)'^ - M^+ie][P + 2mV -l-q^/A + ie]' 

where V = {p' +p)/2m. In the above equation, we have replaced the bare nucleon and pion masses and 
the p ion decay constant by their physical values. The corresponding NN potential V^'^\q ) is defined as 



|18,^ : 



VS'\q)=T^c'\t) 



(3.66) 

Neglectin g ter ms proportional to q"^ /m?, the spectral function /9^'(^) corresponding to Eq. (|3.64|1 takes 
the form [l85|: 

Qfl2 n 11^2 _ ..2 

P''S\^^) = (2M2(2ci - C3) + C3/z2) ^ arctan:^ , (3.67) 
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where we have introduced the abbreviation 



For the sake of simphcity, we restrict ourselves to the SFR cutoff A = cxd. The choice of the cut-off is of 
no relevance for the present discussion. Since x ~ m/Q, the standard HB approach corresponds to the 
expansion of the arctanx in the above equation for large x jH^: 

TT 1 1 , , 

arctanx = — h + . . . . (3.69) 

2 X 6X'^ 

Keeping only the first term in this series, one reproduces the HB expression for the spectral function 
given in Eq. 1)3. 22|) . The series in Eq. (|3.69|) . however, only converges for |x| > 1. This condition is not 
met for ^ in the threshold region 2M^ < /i < 2M^(1 + A), where A = Ml/iSm^) + 0{M^/m^).*'^^ The 
HB expression for the spectral function is, therefore, not correct in this region, as pointed out in Ref. [j^ ]. 

How does this failure of the HB expansion affect the NN potential? In order to answer this question, it is 
useful to switch to coordinate space. For r > 0, the inverse Fourier-transform can be expressed in terms 
of the spectral function: 



1 



VF'ir) = ^ / d^Me-'^>^^'(/i) (3.70) 



1 /■2M^(1+A) 



2vrV ./2M^ 27r2r j2Af^(i+A) 



The first integral in the second line of the above equation goes over the range of /i, where the spectral 
function is incorrectly described in the HB approach. The integration in the second term goes over the 
interval of ^u, where the HB expansion is valid (for not too large values of /x). In general, the contribution 
of the first integral is expected to be suppressed relative to the one from the second integral by a factor 
M'^/rm? and thus can be safely neglected at the considered order. This suppression, however, does not 
occur at asymptotically large distances, i.e. for r > m? /M'^, where the potential is determined by the 
threshold behavior of the spectral function. On the other hand, the 2PE potential falls off exponentially 
and becomes very weak at large distances. Consequently, the problem with the formal inconsis tency of 
the HB approach is expected to be of little relevance for practical applications. Recent work Il86l | on 
peripheral NN scattering based on the relativistic approach confirms this expectation. 

To see how the above qualitative arguments and estimations work in practice, we again follow Ref. |l85l ] 
and calculate the leading and next-to-leading order HB approximations for the potential in coordinate 
space which correspond to the first and the first + second terms in the series in Eq. ()3.69() . In Fig. 
we show the ratios of these leading and next-to-leading order HB potentials to the relativistic result 
corresponding to the spectral function in Eq. 1)3. 67(1 . Here, we adopt the same values for the LECs ci^s as 
in section TS. 2. 11 As expected, the deviations from the correct result increase for large values of r. Notice 
that large deviations at short distances r < 1 fm are due to the fact that the HB expansion diverges at 
large ^. The leading HB approximation deviates by ~ 20% from the relativistic result in t he re gion of 
physical interest. For further discussion of these and related issues, the reader is referred to |l8f)l |. 

To conclude, the formal inconsistency of the HB approach strongly affects the behavior of the potential at 
asymptotically large distances. In this region, however, the potential is very weak and its contribution to 

*^^This condition is also violated for large /i of the order ^ ~ 2m, which is, however, outside of the validity region of the 
HB approach and therefore of no relevance. For the role of the large /x-components, see also the discussion in section 1^.2. II 
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Figure 13: Ratios of the leading (dashed line) and next-to-leading order (solid line) HB potentials to 
the relativistic result versus the distance r as explained in the text. 



observables is negligible. At physically interesting distances, the effects due to the formal inconsistency 
of the HB expansion are small and expected to be irrelevant at the consid ered order. The 2PE potential 



1871 . Il88l | using the formalism 



has been derived at order u = 4 in t he re lativistic version of chiral EFT 
based on the S-matrix, see also Ref. |18(t| | for the application to peripheral NN scattering 



3.3 Few— nucleon forces 



Two-nucleon forces discussed above provide the most important contribution to the Hamilton operator. 
Three- and more-nucleon forces are suppressed compared to the two-nucleon ones (2NFs) by powers of 
Q/A and thus appear as corrections. It is important to take these corrections into account in order to 
understand the properties of few-nucleon systems at the quantitative level. The strength of the chiral 
EFT approach is that it provides a framework to derive few-nucleon forces in a systematic way and fully 
consistent with the 2NF. 



3.3.1 Three— and four— nucleon forces in the isospin limit 

According to the power counting in Eq. (|2.47|) , three- and four-nucleon forces (4NFs) start to contribute 
at order u = 2, i.e. they are suppressed by a factor (Q/A)^ compared to the leading 2NF. The leading 3NF 
contribution arises from tree diagrams constructed with the order Aj = vertices, see graphs (a), (b) and 
(c) in Fig. ll4[ Th e lea ding 4NF is represented by the disconnected tree diagrams (d) and (e). As pointed 
out by Weinberg |ll3l |. graph (a) does not contribute at order u = 2 due to an additional suppression 
factor Q/m. The origin of this suppression is easily understood in terms of "old-fashioned" time-ordered 
perturbation theory. Since graph (a) does not give rise to reducible topologies, its contribution to the 
3NF is given by the sum of all possible time-ordered graphs, which build up the corresponding Feynman 
diagram. Since energy is conserved at each vertex of a Feynman graph, the time derivative, which enters 
the Weinberg-Tomozawa ttitNN vertex, yields a difference of nucleon kinetic energies which scales as 
Q'^/m instead of Q. 
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(a) (b) (c) (d) (e) 

Figure 14: Leading contributions to 3NF (graphs (a)-(c)) and 4NF (graphs (d) and (e)) at order v = 2 
that vanish. For notation see Fig. [7| 




Figure 15: Time-ordered 2PE diagrams at order v = 2. The first and second rows show all possible 
irreducible topologies, while the third row subsumes reducible topologies. Graphs resulting from the 
interchange of the vertex ordering at the middle nucleon line are not shown. For notation see Fig. [Tj 



The remaining graphs (b)-(e) in Fig. lead to vanishing contributions to the nuclear force when the 
latter is defined w i thiii an energy-independent formulation such as the method of unitary transformation, 
see e.g. |l32l . Il8^ . ll9Cl| |. Consider, for example, the 2PE 3NF resulting from diagram (b) which can be 
obtained evaluating the corresponding matrix elements of the operators in Eq. (|3.1(i|) . The first term 
in this expression gives rise to irreducible time-ordered graphs (i.e. the ones without purely nucleonic 
intermediate states) shown in Fig. Eland defines the 3NF potential in the Tamm-DancofF method, which 
can be expressed schematically as 

4 _^ 2 _^ 2 

_UJl{uJi + UJ2)UJ2 UJi{uJl+U)2) (Wl + W2)W2 

where u)i = Hq = \Jk'^ + denotes the pion energy and we pulled out the common factor M representing 
the spin, isospin and momentum structure, which is obviously the same for all graphs in Fig. 1151 The 
remaining two terms in Eq. (|3.16|) are specific for the method of unitary transformation and correspond 
to reducible topologies shown in Fig. 1151 We remind the reader that although these diagrams contain 
purely nucleonic intermediate states, their contributions are not enhanced in the limit m — > oo and 
cannot be identified with the iteration of the potential in the dynamical equation. These diagrams are 
thus not truly reducible in the sense of "old-fashioned" time-ordered perturbation theory. The resulting 



M 



-2 ' 



(3.71) 
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(a) (b) (c) 

Figure 16: Three-nucleon force at order v = 3. For notation see Fig. [7| 



contribution to the potential reads: 



LO1UJ2 



UJ1UJ2 



M 



(3.72) 



The cancehation between irreducible and reducible diagrams in the method of unitary transformation is 
now evident. ^^"^ The same sort of cancellation occures for the remaining graphs (c)-(e) in Fig. 1141 In 
other words, the entire contribution to the scattering amplitude represented by the Feynman diagrams 
(b)-(e) in Fig. El is reproduced by iteration of the leading 2NF in the corresponding dynamical equation 
with no need to introduce additional few-nucleon forces. We therefore conclude that the 3NF and 
4NF at order v = 2 vanish completely. 

The first non-vanishing 3NF contribution appears at order = 3, i.e. at N^LO, and arises from diag rams 
shown in Fig. 1161 The contribution from graph (a) in this figure is given (in the 3N CMS) by jl92l |: 



(3) 



2lT 



E 



9A 
2F^ 



where qi = pi 



(g-;2+M2)(g,2 + M2) 

Pi] Pi [pi ') is the initial (final) momentum of the nucleon i and 
4ciM2 



(3.73) 



F 



0/3 



ijk 



i?2 



+ 



2C3 



9i 



p2 ^ 



[qi X Qj] ■ 



The subscripts of the Pauli spin and isospin matrices refer to nucleon labels. The form 1)3. 73(1 was 
shown to match with the low-momentum expansion of various existin g phe nomenological 3NFs provided 
they respect chiral symmetry. This issue is extensively discussed in |l93l |. The contributions from the 
remaining graphs (b) and (c) in Fig. 1161 take the form (192, .194] 



y(3) 
Itt, cont 



t^j^k 



D 



+ M2 



(3) 



cont 



Tfej 



(3.74) 



where D and E are the corresponding LECs from the Lagrangian of order A = 1. We stress that the 
proper incorporation of the Pauli principle allows to substantially reduce the number of independent 
operators yielding only two terms in Eq. 1)3. 74(1 194]. Notice further that no 4NFs appear at order u = 3. 

The first corrections to the 3NF as well as the first contribution to the 4NF arise at order 1/ = 4. Some 
examples of diagrams which contribute at this order are shown in Fig. I17L For the 3NF, one has to take 

*^®Notice that although the above canceUation at the level of the nuclear potential does not take place in the Tamm-Dancoff 
method, the resulting non-vanishin g 3NF wa s shown to cancel against the recoil corrections to the 2N potential upon the 



iteration in the dynamical equation |191l . ll92 
*^''This statement does not apply to nuclear forces defined in the Tamm-Dancoff method. In that case, the IPE 2N 
potential receives 1/m-corrections, which are absent in the method of unitary transformation and require additional few- 
nucleon forces in order to cancel the corresponding additional contributions to the scattering amplitude. 
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Figure 17: Examples of the three- and four-nucleon force contributions at order = 4. For notation 
see Figs. [3 



Chiral order 


2N force 


3N force 


4N force 


1^ = 


^Itt + Kont 






u = 1 








v = 2 


Vitt + V2tt + Vcont 






= 3 




^27r + ^Itt, cont + Kont 




1/ = 4 


Vl^ + V2n + + Fcont 


work in progress 


work in progress 



Table 1: Isospin-symmetric nuclear forces up to order v = A. 



into account all possible one-loop graphs constructed with the lowest-order vertices and tree diagrams 
with one insertion of the order Aj = 2 vertices. In addition, one has to include the leading 1/m-corrections 
to diagrams (a)-(c) in Fig. 1141 The leading contributions to the 4NF arise at order = 4 from connected 
tree diagrams with the lowest-order vertices. One should also consider various disconnected diagrams 
which might, in principle, also contribute to the 4NF at this order. Work along these lines is unde rway . 
The leading 1/m-corrections to the static 2PE 3NF were already evaluated and can be found in Ref. jl4fl |. 
The relative importance of various isospin-invariant contributions to few-nucleon forces discussed in this 
section and in sections 13.2.21 1X2. 3l and r3. 2. 41 are summarized in Table^based on the chiral power counting 
in Eq. (|2.47jl . We stress that these results rely on the energy-independent formulation. 



3.3.2 Isospin breaking corrections 

Let us now discuss isospin-breaking corrections to the 3NF which were recently worked out up to the 
same order = 5 as the corresponding 2NFs. No isospin-breaking four- and more nucleon forces appear 
up to this order. Similar to the 2NF, it is useful to classify the 3NF with respect to its isospin structure. 
Following the lines of Ref. jl95l |. one distinguishes between the following three classes: 

Class I: Vi = Ei^j^k + f^'i'' '^i ' '^j + ['^i x "^jl " '^k) , 

Class II: Vu = Y.^+,+k («/^ + [r. x r,f r|) , (3.75) 

Class III: Vni = E.^.y. [^ti rf + PYu [r^ x r,f + 7^, rf r, • r, + rf r| r^) , 

where a, P, 7 and k are space and spin operators and the indices i, j, k refer to the nucleon labels. 
The class (I) forces are isospin invariant while the class (II) 3NFs break isospin but respect charge 
symmetry. Finally, the class (III) forces are charge-symmetry-breaking. Contrary to the commonly used 
classification scheme in the 2N sector, see Eq. 1)3. 49|) . conservation of the operator with T being the 
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Figure 18: Isospin-breaking 3NF at orders = 4 (graphs (a)-(d)) and z/ = 5 (graphs (e)-(g)). Graphs 
resulting by the interchange of the vertex ordering are not shown. For remaining notation see Figs. [7| [TTl 
andini 



total isospin operator, which ensures that there is no isospin mixing, is not used in Eq. 1)3. 75() . This 
is because this property depends on the number of particles in the system under consideration. For 
example, all isospin-breaking two-nucleon forces, which do not cause isospin mixing in the two-nucleon 
system, lead to isospin mixing in the three-nucleon system. 

The leadin g arid subleadi ng isospin-violating 3NFs arise at orders u = 4 and z/ = 5 from diagrams shown 
in Fig. [THl jlQSl . Il96| . Il97l |. Graphs (a)-(c) result from a single insertion of the proton-to-neutron mass 
difference into the order v = 2 3NF diagrams shown in Fig. 1141 Contrary to the corresponding isospin- 
invariant contributions discussed in section 13.3.11 diagrai ns fa )-(c) in Fig. yield nonvanishing 3NFs. 
The origin of this difference is explained in detail in Ref. jl2Rl |. One finds the following class (III) 3NFs 
resulting from diagrams (a) and (c): 



(4) 



2n 



25m 



QA 



2FJ {e+M^y{q,^+M^) 



[qi X qj] ■ dk [Ti X Tj 



(4) 



Itt, cont 



25mCT 



9A 



2FJ {e+MiY 



[Tfc X Ti\ [Gj X 0-fcJ • qi . 



(3.76) 



The contributions from the 27r-exchange diagram (b) can naturally be combined with the contributions 
from graphs (d) and (g), which have the same structure. The resulting 3NF reads: 



v.. 



(4,5) 



27r 



9A 
2K- 



[aj ■ qi){aj • qj) 

(%2 + M2)(g,2 + M2) 



(Sm 



istr 



4F2 



2{Ti • Tk)Tf 



(3.77) 



While the term oc {6mY^^ breaks charge-symmetry, the contribution proportional to the LEG /i and 
arising from diagram (g) is charge-symmetry conserving, i.e. class (II). Isospin-violating t tttNN vertex 
oc /i plays an important role in the analysis of isospin violation in pion-nucleon scattering |l2Cl| | and the 
evaluation of the ground state characteristics of pionic hydrogen [l2J]. In the 2N sector, it only leads 
to an isospin-invariant contribution to the 2PE potential, which has the same form as the ci-term in 
Eq. 1)3. 38() . The charge-symmetry-breaking part of V2^'^^ in Eq. (|3.77|1 was also obtained by Friar et 
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(a) (b) (c) (d) (e) 




(f) (g) (h) (i) 

Figure 19: Leading (a-c) and subleading (d-i) isospin-violating contribution to the 3NF which vanish, 
as discussed in the text. Graphs resulting from the interchange of the vertex ordering are not shown. For 
notation see Figs. [3 and El 



al. jl97j | using a completely different approach. In their method, the neutron-to-proton mass difference, 
which is inconvenient to handle in practical calculations, is replaced by a seri es of new isospin-breaking 
vertices in the Lagrangian via an appropriate field redefinition, see Ref. jl75l | for more details. Finally, 
diagrams (e) and (f) due to the pion mass difference lead to the following class (II) 3NF at order u = 5: 



Vr 



(5) 



2n 



QA 



(5i ■ qi){aj ■ q.j) 



2FJ ie+M^y{q,^+M^) 




p2 ^ p2 



(5) 



Ivr, cont 



2 9A 



D 



8F2 (gl2 + M2)2 



(3.78) 



These results are consistent with the isospin-invariant 3NFs V2J and V^iVcont ™ ^<is- (|3.73|) . (|3.74|) being 
expressed in t erms of the charged pion mass. The corresponding expressions in coordinate space can be 
found in Ref. Iioi. 



In addition to graphs shown in Fig. 1181 diagrams (a)-(c) and (d)-(i) in Fig. formally contribute 
to the leading = 4) and subleading (u = 5) isospin-breaking 3NF, respectively. Their pertinent 
contributions, however, vanish. In particular, for graphs (a), (b), (d), (f) and (i) one observes the same 
sort of cancellation between various time orderings as in the case of the corresponding isospin-invariant 
3NFs, see section TJ. 3. II for more details. The contributions of diagrams (c) and (e) are suppressed by a 
factor of Q/m due to the time derivative entering the Weinberg-Tomozawa vertex. Explicit evaluation 
of the remaining contributions of graphs (g) and (h) can be performed along the lines of ref. ^95] and 
yields vanishing result. 

Finally, we point out that there are many 1/m-corrections to the obtained results which start to con- 
tribute at order z/ = 6. Notice, however, that if one would adopt the counting rule m ~ A, various 
1/m-corrections (including the ones due to virtual photons) would contribute already at order {Q/Kf'. 
Some 3NF diagrams due to virtual photon exchange were considered by Yang and found to provid e rel - 
atively small contributions of the order of ~ 7 keV to the ^He-^H binding-energy difference jl98i . il99l |. 
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2N fnrrp 




1^ = 






v = l 


- 


- 


u = 2 






1^ = 3 


Vln + Vcont 




1/ = 4 


K-7 + Vitt + V2tt + V^ont 


^277 + ^Itt, cont 


1^ = 5 


Vitt + V27r + K;ont 


^277 + ^Itt, cont 



Table 2: Isospin-breaking two- and three-nucleon forces up to order v = 5. 



Furthermore, we remind the reader that the long-range electromagnetic 3NFs might, potentially, cause 
large effects in scattering at low energy. The relative sizes of various isospin-breaking contributions to 
the two- and three-nucleon forces discussed in this section and in section ri.2.51 are summarized in Table 
[21 We stress that these results rely on the energy-independent formulation and on the counting rules 
specified in Eqs. (ITlIlll . (ITTtI) and (10^1) . 



3.4 Role of the A— excitation 



The A-isobar is well known to play an important role in hadronic and nuclear physics. This is due to its 
low excitation energy. Am = ttia — m = 293 MeV, and strong coupling to the vr system. Because of the 
small value of At w, it is not clear a priori whether in EFT it should be included explicitly, treating Am as 
a small quantity |200|, or integrated out. On the one hand, inclusion of A yields a scheme, which differs 
from the chiral expansion since Am does not vanish in the chiral limit. On the other hand, it might 
provide a useful phenomenological extension and a systematic power counting has already been worked 



out assuming Am ~ |2flll ]. see also |2n2l | for an alternative scheme. In the chiral EFT discussed so far, 
the eff ects of the A's are only taken into account implicitly, i.e. through the values of the corresponding 
LECs j203j| . We will now discuss the implications of treating the A as a dynamical degree of freedom in 
the EFT. 

The leading contributions to the nuclear force due to the A-excitation can be obtained in the heavy— 
baryon approach from the following terms in the Lagrangian [lift . i2nii] 

£ = A^ {ido - Am) A + ^ (^iV^STA + h. c.) • Vtt - DrN^arN ■ (^N^STA + h. c.) , (3.79) 

where Ha and Dt are LECs and A is a four-component spinor in both spin and isospin spaces that 
represents the delta. Further, Si [Ta) are the 2x4 spin (isospin) transition matrices which satisfy the 



relations SjSj = {26ij — ieijkO'k)/3 {TaT^ = {26ab — KabcTc)/^) |204| |. Throughout this section, we will 
assume exact isospin symmetry. Notice that the large mass scale disappears from the Lagrangian 
after performing the heavy-baryon expansion. Only the small scale Am enters the resulting Lagrangian. 

The leading contributions to the 2NF due to intermediate A excitations arise at NLO, u = 2 , and are 
shown in Fig. I2UI In the context of chiral EFT, they were first discussed by Ordonez et al. |ll9| | based on 
"old-fashioned" time-ordered perturbation theory.*^^ These contributions were reconsidered by Kaiser 

#18 



Certain contributions to the potential were considered much earher, see e.g. jl& 
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Figure 20: Leading contributions to the 27r-exchange potential with single and double A-excitations. 
Double lines represent the A~isobar. For remaining notation see Fig. |7| 



et al. |205| using the Feynman graph technique. For completenes s, we collect here the expressions for the 
corresponding non-polynomial pieces of the 2PE potential from \2Q^ (generalized to SFR with arbitrary 
A) for the three distinct groups of terms: 

• A-excitation in the triangle graphs: 

^irtriangie = "^^^^ (^1 " ^2) { (6i? " u;')L^{q) + l2{Amf E (q)'^ , (3.80) 

with 



" <''' = ASlL,.^^""""'"^'^2ASr^- E = 2MiW-2(^''^)'. (3.81) 
• Single A-excitation in the box graphs: 

(n • T2) I (12(Am)2 - 20M2 - lW)L^{q) + QE^D^{q) 



8647r2F4 

' (<?i • q){a2 ■ q) - q\Si ■ ^2)] \ -2L^{q) + (t^' - 4{Amf)D^{q) 



9a "'/_ 

'Ti-T2)i{ai-q){a2-q)-qH^l-ff2))oJ^A^iq). (3.82) 



1927r2F4 

.2 1,2 



9\K 



• Double A-excitation in the box graphs: 

^I:box-2 = -^^^l-4{AmfLHq) + E[HHq) + {E + 8{Arnf)D'Hq)]} (3.83) 

(ri • T2) <! {12E - uJ^)L^{q) + 3E[H^{q) + (8(Am)2 - E)D^{q)] 



77767r2 Fi 



' (<?! • q){a2 ■ q) - g'(<?i • <?2) ) <! 6L^(g) + (12(Am)2 - LO^)D''{q) 



34567r2F4 
207367r2Fi 



with 



<?i • g )(a2 • q) - q\ai ■ <T2)J |2L^(g) + (4(Am)^ + co^jO'^iq) 

L^{q) - L^(2V(Am)2 - M^) . (3.84) 



u;2 - 4(Am)2 
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Figure 21: Leading contributions to the 3NF due to explicit A's. For notation see Figs. 171 and I2UI 



Notice that some of the above contributions, in particular, those ones corresponding to the first and 
the last lines in Eq. H3.82() . have precisely the same structure as the corresponding N^LO terms in EFT 
without A's, see Eq. (|3.38|) provided one chooses 



- C3 = 2C4 



9Am 



(3.85) 



Using the large-A'^c v alue Ha = 3gA/V^ and comparing the LECs in Eq. (|3.85|) with the ones in EFT 
without explicit A's |l54| |. one concludes that the A provides the dominant (significant) c ontri bution 
to C3 (C4). Resonance saturation for these and other ttN LECs is discussed in detail in '203*]. The 
remaining contributions to the NN potential due to intermediate A~excitation in Eqs. (|3.8nj) - (|3.83j) show 
a non-trivial dependence on the A^A mass splitting. We further emphasize that s ome of the subleading 
corrections to the 2PE potential with intermediate A- excit ations are discussed in |ll9l | while the leading 
relativistic corrections of order u = A can be found in 120511. 



Let us now switch to the leading 3NF contributions due to intermediate A excitations, which also aris e 
at NLO and are shown in Fig. 1211 In the context of chiral EFT, they were first discussed in Ref. jl92l |. 
We remind the reader that in the EFT without explicit A's, the first nonvanishing contribution to the 



3NF in the energy independent formulation only arises at N LO, 



3, see section 13.3.11 The 2PE 



contribution from the first graph in Fig. |^has the same form as the corresponding terms in the theory 
without delta, see Eq. 1)3. 73() , provided one chooses ci = and 03^4 according to Eq. (|3.85() . Interestingly, 
the c ontributions of the two other diagrams in Fig. 1211 vanish as a consequence of the Pauli principle, see 
I194II for a related discussion. 



To summarize, the effects due to intermediate A excitations in the EFT with explicit A's first show up at 
NLO (NLO-A). At this order, one finds important contributions to the 2PE 2N and 3N potentials. In the 
EFT without explicit A's, the major part of these terms is shifted to N^LO (N^LO-^^) and represented 
by the contributions proportional to LECs 03^4 which are saturated by the A. In particular, the leading 
isoscalar central and isovector spin-spin and tensor components in Eq. (|3.82|) are exactly reproduced 
at N^LO-^^ For the 3NF, the functional form of the NLO-A contribution is completely recovered at 
N^LO-.^. We stress, however, that the strength of the 2PE 3NF at N^LO-^^is overestimated by a factor 
of 4/3 cor npar ed to the one at NLO-A if the values of the LECs 03^4 are fixed from vrA^ scattering at 
threshold \2Q(y^ . Clearly, higher-order counter terms remove this discrepancy. 

Despite the fact that many of the NLO-A contributions to the nuclear force are reproduced at N^LO-^ 
it might be advantageous to treat the delta isobar as an explicit degree of freedom. This might help to 
increase the range of applicability of the EFT. Further, one expects that such an approach would lead 
to the contributions to the nuclear force of a more natural size. For example, the large portion of the 
2PE potential oc Cj at N^LO-^ which is known to be very strong, is shifted to NLO in the EFT with 



explicit A's. Similar effects are also observed for isospin-breaking terms, see We emphasize, 



* ^^Antisymmetric (with respect to an interchange of nucleons i and j) nature of few-nucleon states is not fuUy incorporated 
in |l92|| . which resuhs in the presence of redundant terms. 
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however, that such an approach with explicit A's is much more compHcated since one has to deal with 
more structures and also needs e.g. to reanalyze pion-nucleon scattering (see Ref. Jio^ for an at tempt). 
For some recent related work in the single-nucleon sector the reader is referred to j208l . l209l . bid ] . 



4 Applications to few— nucleon systems 

4.1 The two— nucleon system 

4.1.1 Peripheral nucleon nucleon scattering 

Nucleon-nucleon scattering in peripheral partial waves withi n the chiral E FT fr amework has attracted a 



IN ucieon-nucleon scattermg m peripn eral partial waves witnm tne cmral tht 1. tr amework nas attracted a 
special interest in the past few years see also Ref. |213| for a related work. This 

has several reasons. First, the contribution of the short-range contact interactions is suppressed due to 
the centrifugal barrier in partial waves with large values of the orbital angular momentum. For example, 
no contact terms contribute to D~ (F-) and higher partial waves up to N^LO (N^LO). Consequently, 
peripheral phase shifts are entirely determined by the long-range part of the nuclear force and thus 
are expected to provide a sensitive test of the chiral 2PE potential. Secondly, the smallness of the 
corresponding phase shifts suggests that they should be describable in perturbation theory. This allows 
one to avoid the additional complication related to the nonperturbative treatment of the Schrodinger 
equation. The validity of perturbation theory is confirmed by the conventional scenario of nucl ear fo rces 
based on the existing boson-exchange models and various phenomenological potentials, see |214l |. It 
should be understood that the weakness of the chiral NN force in peripheral channels, which justifies the 
use of perturbation theory, may only be verified assuming a particular regularization scheme for the chiral 
potential, which shows a singular behavior at the origin. Using, for example, dimensional r egula rization, 
the iterated IPE potential was found to be numerically small in most peripheral channels |142l |. Notice, 
however, that the weakness of the NN interaction for high values of / is not related to the chiral expansion 
and does not follow from the chiral power counting. We further stress that the role of peripheral NN 
scattering and the resulting constraints on the values of the LECs should not be overestimated since the 
observables at low and moderate energies are, in general, only weakly affected by peripheral phase shifts. 

We will n ow consider D- and higher partial waves up to N^LO in chiral EFT following the lines of 



Ref. 143]. Using Born approximation for the scattering amplitude, the phase shifts and mixing angles 



in the convention of Stapp et al. |21Rl | are determined by the 2N potential as: 

St' = -^{lsj\V2^\lsj) , e, = -^(J - l,sj\V2N\j + l,sj) , (4.1) 

where s and j refer to the total spin and angular momentum, respectively, and q is the nucleon CMS 
momentum. Clearly, such an approximation violates unitarity. This violation is small provided that 
the corresponding phase s hifts are small. Alternatively, one can use the X-matrix approach in order to 



restore unitarity, see e.g. 21l|. Here and in what follows, we adopt the same notation for the matrix 
elements in the \lsj) basis as in Ref. |l62l |. The formulae for the partial wave decomposition can be 
found in appendix B of this reference. Up to N^LO, the potential V2N consists of the Ivr-exchange and 
the leading and subleading 27r-exchange terms given in Eqs. (|3.33|) . H3.36() and ()3.38|1 . respectively. The 
contact interactions do not contribute to D-waves at this order. For the sake of simplicity, we neglect 
isopi n-violating effects in this section. In what follows, we adopt the same values of the LECs as in 
Hsi]: gA = 1-26, = 92 A MeV and = 138.03 MeV. For the Goldberger-Treiman discrepancy in 
Eq. H3.33() . we use 6 = 0.03, which leads to g^rN — 13.2. In addition, we have to specify the values of 
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Ref. 


Cl 


C2 


C3 


C4 


^ AT 


12171 


-0.64(14) 


1.78(20) 


-3.90(9) 


2.25(9) 


7tN, Q3 


[203] 


-0.93(10) 


3.34(20) 


-5.29(25) 


3.63(10) 


ttN, Q3 (fit 1) 


u 


-1.23(16) 


3.28(23) 


-5.94(9) 


3.47(5) 


ttN, Q3 (fit 1) 


[154] 


-0.81(15) 


9.35(66.7) 


-4.69(1.34) 


3.40(4) 


NN 


[21JJ 


-o.sr 


3.28'^ 


-3.40 


3.40* 


pp 


[219] 


-0.76^ 




-5.08(28) 


4.70(70) 


NN 


[2201 


-0.76* 




-4.78(10) 


3.96(22) 



Table 3: Recent determinations of the LECs q. All values are in GeV . The values of the LECs used 
as input are marked by the star. 



the LEG ci,3.4 which determine the strength of the subleadin g 2P E potential. Table IHl summarizes some 
of the modern determinations of these LECs, see also Ref. |216l ] where a different notation was used. 
We also include the LEG C2 in Table IHl which first contributes to the potential at order u = 4. For a 
precise meaning of the indicated errors in the various determinations, the reader is referred to the original 
publications. Presumably, the most r eliab le determination of the Cj's from the vrA^ system in the third 
order HBCHPT is performed in Ref. |l54l ]. In this work, the vrA^ scattering amplitude is reconstructed 
in the unphysical region using dispersion relations, where the chiral expansion is expected to converge 
rapidly. Unfortunately, this method does not allow for a reliable determination of C2. Recently, the 
Nijmegen grou p inco rporated the chiral 2PE potential up to order i/ = 3 in an energy-dependent PWA 
of the pp data |219| | and pp + np data |22fll ] and was able to extract the values of the LECs 03^4, see 
Table 01 The resulting LECs 03^4 are consistent with the Q^-determinations in the ttN system (with the 
value of C4 being on the upper side). We stress that the errors indicated in these references are statistical. 
It is not fully clear how the complet e the oretical uncertainty can be estimated in this approach, which is 
not entirely based on EFT, see ;iso M for a related discussion. A somewhat smaller value of the LEC 
C3, C3 = —3.4 GeV~^, was found in Ref. |21l| to be consistent with empirical NN phase shifts as well as 
the results from dispersion and conventional meson theories. In that work, F- and higher partial waves 
were studied in chiral EFT at order = 4. Notice that the phase shifts depend most sensitively on C 3 
and are less sensitive on variations of ci^2,4- Similar smaller values for C3 were also found by Higa |l8fi| |. 
who looked at peripheral NN waves in the relativistic version of chiral EFT at order v = A. Interestingly, 
similar values for the LEC C3 were also extracted recently from matching the chiral expansion of the 
nucleon mass to lattice gauge theory results at pion masses between 500 and 800 MeV 152]. Given the 
present uncertainty in the v alue of C3, we will show the results corresponding to the following two choices: 
the central value from |l54l ]. C3 = —4.7 GeV~^, and the value from Ref. ^2 111] . C 3 = — 3.4 GeV~^. For the 
LECs ci_4 we adopt the central values from the Q^-analysis of the ttN system |l54| : ci = —0.81 GeV~^, 
C4 = 3.40 GeV-^ 

Let us begin with the D-waves which are shown in Fig. 1221 The LO result represented by the pure IPE 
potential already provides a good approximation to the phase shifts in the ^Di and partial waves. 
It is too weak in the ^D^ channel and does not describe properly the ^1)3 phase shift. The latter one 
appears to be quite small {\5\ ~ 4.6° at E'lab = 300 MeV) compared to the other D-wave phase shifts 
(|(5| ~ 9.7° . . . 25.5°). The reason is that the partial-wave projected IPE, taken on the energy shell, is 
strongly suppressed in this channel. Consequently, the ^D^ phase shift is sensitive to 2PE but also to 
the iteration of the potential which is neglected here. The NLO predictions show a visible improvement 
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Lab. Energy [MeV] Lab. Energy [MeV] 



Figure 22: D-wave NN phase shifts. The dotted curve is the LO prediction (i.e. based on the pure IPE 
potential). Dashed, dash-dotted and solid curves are the NLO, N^LO with C3 = —3.4 GeV~^ and N^LO 
with C3 = —4.7 GeV~^ results based on the DR potential. The dashed, light shaded and dark shaded 
bands refer to the NLO, N^LO with C3 = -3.4 GeV"^ and N^LO with C3 = -4.7 GeV~^ results using the 
SFR with A = 500 . . . 800 Me V. The filled circles (open triangles) dep ict the results from the Nijmegen 
multi-energy PWA [1^, E3 

(Virginia Tech single-energy PWA [223]). 



for the phase shift but go into the wrong direction in the and "^D-^ channels. The N^LO results 
based on DR are depicted by the dash-dotted and solid lines for the two choices of the LEG C3 specified 
above. The reasonable agreement with the data observed at LO and NLO is destroyed in all partial 
waves for energies E'lab > 50 MeV and the chiral expansion does not seem to converge. N otice tha t 
similar deviations from the data are also observed for mixing angle €2 which is not sho wn here 142l . 147 1 . 
This disagreemei it w ith the data for the DR potential was first pointed out in |l42| |^^^ and then also 
reported in Ref. 147l |. We stress that the results presented here are parameter-free. As pointed out by 



Kaiser et al. 142l |. the origin of the strong disagreement with the data can be traced back to the central 
part of the subleading 2PE potential, see Eq. (|3.38)) . which is proportional to the LEG C3 and shows an 
unphysically strong attraction at intermediate and short distances when DR is used to regularize the 
divergent integrals. This is demonstrated explicitly in Fig. El where the chiral potential is compared to 
phenomenological a {a + cu + p) contributions. Notice that another consequence of the strong attraction 
of the subleading 2PE potential is given by the unphysical bound states in low partial waves which could 

There are several minor difFerences between the results presented here, which are based on Ref. |l47ll . and the ones 
shown in |l42ll . In particular, this latter work uses relativistic kinematics when calculating phase shifts, slightly different 
valu es of the L ECs and also incorporates the contribution from the once iterated IPE potential. The numerical results of 
|l47|| and |l42| are, however, rather similar. 
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not be avoided in the N^LO analysis of Ref. [l39|. Although such spurious bound states do not influence 
NN observables at low energies, they are at leas t incon venie nt for the application to few-nucleon systems 



and might influence processes like e.g. Nd 214i | or -Kd |223l | scattering 



The strong attraction of the subleading 2PE potential can partially be explained by the large values of 
the LECs Cj. For example, the absolute value of the dimensionless coupling corresponding to t he L EG C3 
is of the order ~ 10 which has to be compare d wit h its expected natural size of the order ~ 1 |203| |. The 
origin of this enhancement is well understood |203l |: the LECs 03^4 are, to a large extent, saturated by the 
A-excitation. This implies that a new and smaller scale, namely — m ~ 293 MeV, enters these LECs 
in EFT without explicit A, see also the discussion in section ITU As shown in Fig. 1^1 the large numerical 
value of C3 appears to be consistent with the results based on phenomenological potential models at large 
distances, where the chiral potential is unaffected by the regularization procedure. On the other hand, 
the behavior of the 2PE potential at intermediate and short distances depends, to a large extent, on the 
way one regularizes the corresponding loop integrals, see section l5.2.11 The unphysically strong attraction 
of the 2PE potential in the region r < l/Mjr ~ 1.4 fm arises from high-momentum components of the 
exchanged pions, which cannot be properly treated in an EFT but whose contribution is, nevertheless, 
included in the potential obtained using DR. In the SFR approach, these high-momentum components in 
the mass spectrum are explicitly removed. As depicted in Fig. El using SFR with the reasonably chosen 
cut-off A = 500 . . . 800 MeV greatly reduces this unphysical attraction and the resulting potential is of the 
same order in magnitude as the one obtained in phenomenological boson-exchange models. The results 
for D-waves are strongly improved using the spectral function regularization instead of dimensional one, 
as documented in Fig. 1221 It should be understood that dimensional regularization is by no means 
ruled out by such consideration. In general, for quickly converging expansions, it can and should be the 
method of choice. If, however, the convergence for some well understood physical reason is slow and 
(some) observables become sensitive to spurious short-distance physics kept in DR, it is preferable to 
use SFR. Choosing DR, one picks up a portion of spurious short-distance physics which, under normal 
circumstances, can be compensated by the counter terms included in the potential at a given order. In 
the case at hand, however, unnaturally large counter terms (contact interactions) are required due to the 
large values of the LECs 03^4. The neglect of higher-order contact terms is, therefore, not justifled and 
results in a large disagreement with the data. It is known that a counter term at order = 4 is able t o 
remove the bulk of the disagreement observed for D-waves in the DR based approach |l39l . llG.'j . l224l . 
The unnaturally large size of the leading counter term in the partial wave was also found in Ref. |212l | 
based on the distorted wave method. On the other hand, applying SFR with the cut-off A of the order 
of the separation scale allows to remove the spurious short-range physics from the potential, which 
results in a more natural size of the counter terms. As demonstrated in Fig. I22| no N^LO counter terms 
need then to be taken into account at N'^LO. In other words, the convergence of the chiral expansion 
is substantially improved usin g SFR with the reasonably chosen cut-off A instead of DR. We refer the 
reader to Refs. |l48l . Il49l . Il50j | where a similar method has been applied to improve the convergence of 
the SU(3) baryon CHPT. 

The results for F-waves are presented in Fig. 1231 Although the situation with the DR subleading 2PE 
potential is much less dramatic compared to D-waves, a too strong attraction is clearly visible in the 
^i^2) ^-^3 and ^i<4 partial waves, especially if one uses C3 = —4.7 GeV~^, see also 1142]. Removing the 
short-distance components of the 2PE potential with the cut-off regularization improves the results in 
the ^-Fs channel, while additional repulsion is still missing in the partial wave. A much smaller 
sensitivity of F-waves to the choice of regularization and thus to short-range physics is, of course, the 
consequence of the stronger centrifugal barrier. The leading contact terms in F-waves are suppressed by 
(Q/A)^ compared to D-waves and appear flrst at order v = Q. Therefore, even if the corresponding LECs 
are large in the approach based on DR, neglecting those terms leads to much smaller effects in F-waves. 
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Lab. Energy [MeV] Lab. Energy [MeV] 

Figure 23: F-wave NN phase shifts. For notation see Fig. 1221 



The dependence of the phase shifts on the SFR cut-off being varied in a certain reasonable range reflects 
the influence of the higher-order contact terms and thus provides a natural estimation of the theoretical 
uncertainty at a given order. The typical uncertainty o f ^ 30% observed in the N^LO predictions at 



E'lab = 300 MeV is consistent with the power counting, see 147l | for more details. In view of this theoretical 
uncertainty, it appears to be impossible to make any preference in favor of the choice C3 = —4.7 GeV~^ 
or C3 = —3.4 GeV~^ based on the D- and F-waves at N^LO. 

Let us now briefly summarize the main results for D- and F-waves at N^LO: 

• The subleading 2PE potential obtained using DR leads to strong disagreement with the data in 
D-waves (F~waves) for E^iab > 50 MeV (^lab > 150 MeV) and for both choices C3 = -4.7 GeV"^ 
and C3 = -3.4 GeV-^ (for the choice C3 = -4.7 GeY'^). 

• Using SFR instead of DR with the cut-off A = 500 . . . 800 MeV strongly (sizably) improves the 
description of the phase shifts in D-waves (F-waves). 

• The theoretical uncertainty in the description of D- and F-waves is sizable at large energy. 

As already mentioned before , F- and higher NN phase shifts were also studied at order v = 4 in both 



the standard heavy-baryon j21ll | and relativistic framework 18u\. The authors of Ref. j21ll | use the 
perturbative approach similar to the one described above but also take into account the contribution of 
the once iterated IPE potential and adopt a different counting rule for the nucleon mass (m ~ A). F- 
wave phase shifts obtained in this work and based on the choice C3 = —3.4 GeV~^, are shown in Fig. 1241 
The results at LO, NLO and N^LO are similar to the ones shown in Fig. ESI The N^LO corrections are 
found to be small in the ^^3 and channels and noticeable in the '^^2 and '^^3 partial waves, where 
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Figure 24: F-wave neutron-pro ton phase shifts from Ref. |21ll ]. The sohd dots and open circles are 
the results from Nijmegen PWA [lOSl . Il09| and Virginia Tech analysis SM99 |222l |. Figure courtesy of 
Ruprecht Machleidt. 



they slightly increase the disagreement with the data. We will discuss NN phase shifts at order u = 4 in 
more detail in section 14.1.31 

Last but not least, we would like to comment on an early attempt of Ref. |214l | to avoid the above 
mentioned difficulties caused by the strong attraction of the subleading 2PE potential. In this work 
based on the DR version of the 2PE potential, smaller in magnitude values of the LECs 03^4, namely 
C3 = —1.15 GeV~^, C4 = 1.20 GeV~^, were adopted in order to reduce the unphysical attraction at N^LO. 
This allowed for a fairly good description of t he N N data with no spurious bound states and enabled 
applications to few-nucleon systems, see e.g. jl94| . Clearly, such an approach is still not satisfactory 
due to its incompatibility with vrN scattering. The framework based on SFR allows to avoid the above 
mentioned difficulties and, at the same time, stay consistent with ttN scattering. It sh ould, therefore, be 
the method of choice. We will show some results based on the reduced 03,4 from H in section^ 



4.1.2 Regularization and renormalization of the Schrodinger equation 

In the previous section, we considered peripheral NN scattering in perturbation theory. To properly 
describe low partial waves, where the interaction is strong, the Lippmann-Schwinger (LS) equation for 
the scattering amplitude has to be solved nonperturbatively. The nuclear potential derived in chiral EFT 
is only valid for small momenta and becomes meaningless in the large-momentum region. The results 
presented in section l3?2l show that it grows with increasing momenta. Consequently, the LS equation is 
ultraviolet divergent and needs to be regularized (and renormalized). The problem of renormalization in 
the non perturbative regime in the context of ch i ral EFT attracted a lot of interest in the past years, see 
e.g. 0, Eli, H Hliii, H Eii, H Eii, H H S Eii, H. Here, one difficulty is due to 
the fact that the two-nucleon scattering amplitude can only be obtained numerically if pion-exchange 
contributions are treated nonperturbatively. In addition, since the nuclear potential is nonrenormalizable 
(in the traditional sense), an infinite number of counter terms is needed to remove all ultraviolet diver- 
gences generated by the iteration of the potential in the LS equation. Notice that this feature is in strong 
contrast to the pion and single-nucleon sectors as well as to the KS W approach discussed in section 12.31 
where all divergences at a given order in the chiral expansion can be removed explicitly by a redefinition 
of a finite number of parameters in the Lagrangian. At first sight, this seems to be in contradiction with 
the power counting discussed in section 1231 which tells us that a finite number of counter terms have to 
be included in the potential at any finite order in the chiral expansion. Moreover, the infinite number of 
counter terms needed to remove all divergences in the LS equation raises the question, whether the EFT 
in the 2N sector still has predictive power. Fortunately, this indeed appears to be the case. We recall that 
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the power counting outlined in section 031 is formulated for renormalized contributions to the amplitude. 
Provided that the renormalized LECs are of natural size as explained in section [231 the contribution of 
higher-order counter terms is suppressed by powers of the low momentum scale Q and thus does not 
need to be taken into account. 



The standard procedure to renormalize the LS equation is based on Wilson's method I23q|] and implies 
the following two steps. First, one solves the LS equation regularized with the finite momentum (or 
coordinate space) cut-off and with the kernel represented by the potential truncated at a given order in 
the chiral expansion. Secondly, the LECs, which accompany the short-range contact interactions entering 
the potential, are determined by matching the resulting phase shifts to experimental data. Based on the 
"naive dimensional analysis", see e.g. 1^3] j one expects the resulting LECs to be of natural size.*^^ The 
contributions of the neglected higher-order terms in the potential are expected to be small and irrelevant 
at the considered order. Notice that the exact form of the regulator used in the LS equation is 
irrelevant (provided it does not affect the long-distance behavior of the potential). There is a certain 
freedom in choosing the value of the cut-off. For the sake of definiteness, let us consider the momentum 
space cut-off A. It is clear that taking a too small A will remove the truly long-dist ance physics and 
reduce the predictive power of the EFT. On the other hand, as argued in 

QiUliii, taking too large 



values of A leads to a highly nonlinear behavior and should be avoided as well. An example of such a 
nonlinear behavior is given by the appearance and accumulation of bound states for regularized 1/r" 
attractive singular potentials. In this case, the renormalization-group behavior of the LECs appears to 
be of a limit cycle type (this feature depends on the regulari zatio n employed js^), so that the "naive 
dimensional analysis" is not applicable anymore, '^^^ see Ref. for a recent related discussion. This 

example is also of relevance for the IPE potential, whose tensor component behaves at short distances as 
Taking the limit A — > oo in the LS equation will, therefore, generate bound states in partial waves 
with arbitrarily high values of the orbital angular momentum / and lead to strong deviations from the 
data. As a consequence, an infinite number of higher-order counter terms will be needed in this limit 
in order to restore the agreement with the data. This problem is, clearly, an artifact of the unphysical 
behavior of the point-like static IPE potential at large momenta (or short-distances). In this region, the 
chiral expansion is not applicable, and the IPE potential does not properly describe the true interaction 
between the nucleons. Using a finite cut-off A of the order of the characteristic hard scale of the theory 
allows to avoid this problem in a natural way. It should also be understood that no improvement in the 
description of the data can, in general, be expected increasing A beyond the pertinent hard scale unless 
new physics corresponding to this scale is properly incorporated in the theory. For a much inore exte nsive 
disc ussion of t hese issues including many explicit examples, the reader is referred to [o^. l237l . l238l |. see 
also llfil . 112 1 for a recent related work. 



Let us further mention that the consistency of Weinberg's approach has been questioned in Refs. [t^I.!?^. 
where it was pointed out, based on perturbative arguments, that the order u = 2 contact term oc 
needs to be taken into account in order to absorb the ultraviolet divergence generated by the iteration of 
the order v = Q potential. This problem with the formal inconsistency of the Weinberg scheme was also 



addressed in Ref. 114 1 based on a nonperturbative approach and found to be present in the ^So-channel 



*^^The determined LECs are, strictly speaking, the bare ones (with respect to the ultraviolet cut-oflt in the LS equation) 
and do not necessarily need to be of natural size. 

*'^'^This, however, does not mean that adding an additional higher-order term to the potential and keeping the values of 
the LECs unchanged will resu lt in a small correction to the scattering amplitude. Due to mixing between operators of 
different chiral dimension |237|I . the LECs have to be readjusted. Only then the effect of adding the higher-order term may 
be expected to be small. 

*'^^This should not be misunderstood in terms of impossibility to renormalize the amplitude in the presence of the l/r" 
potenti als w ith the regulator being completely removed (i.e. in the limit A oo). In fact, the opposite has been demonstrated 
in Ref. I23I. 
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while absent in the coupled ^S'l-^Di-channel. In that work, a new p ower counting scheme was invented 
based on the expansion of the potential about the chiral limit, see |ll4l | for more details. It remains, 
however, unc lear whether the problem with the formal inconsistency of Weinberg's approach raised in 
Ill4l | and based on the requirement that the amplitude is renormalizable with the regulator 
being removed, is relevant for calculations within the finite cut-off regularization framework as outlined 
above. Stated differently, it has not been shown that the renormalized LECs corresponding to contact 
terms oc take unnaturally large values within this renormalization scheme. On the con trary, the 
phenomenological success of Weinberg's approach and the qualitative arguments of Ref. jlisj ] based on 
perturbation theory suggest that these LECs are of natural size. 



Finally, we would like to emphasize that various alternative regularization schemes were also considered 
in the context of NN scattering. In particular, subtracti ve procedures are d iscussed in 0,0, IH^, 226l . 228l |. 
boundary condit ion r egularization scheme is applied in 124,^ 1 and dimensional regularization 

is considered in 227 [. In addition, a higher-derivative regularization scheme was applied in Ref. |244l |. 
see also |245| | for an early application of this method in the context of the nonlinear sigma model. 



4.1.3 Two nucleons up to N'^LO 

The first quantitative ap plicat ion of chiral EFT in Weinberg's formulation in the 2N sector was per- 
formed by Ordonez et al. This early study was based on the energy-dependent NN potential and 
incorporated the effects of the A-resonance. The calculations were performed in configuration space up 
to N^LO. Exponential cut-offs were used in order to regularize the divergent loop integrals entering the 
NN potential and the ones arising from iterations in the LS equation. A local form of the employed 
regulator induced effects of contact interactions in all partial waves. Global fits to Nijmegen PWA with 
parameters, many of them being redundant since Fierz reordering was not used, were performed 
for several choices of the cut-off and first numerical results for phase shifts and deuteron properties were 
obtained. The determined values for the the axial pion coupling and the pion decay constants are in 
good agreement with the experimental numbers. For more details on this w ork b v Ordonez et al., which 
is an important milestone, the reader is referred to the original publication 

Park et al. |24fi| | considered a series of interesting applications, mostly related to the deuteron proper- 
ties. They restricted themselves to one-pion exchange and contact interactions for the relevant phases 
^50,3 sJd, and the mixing parameter I 

The first complete N^LO analysis of neutron-proton scattering and the deuteron properties in the EFT 
without explicit A's a nd b ased on the energy-inde pend ent potential, obtained using the method of 
unitary transformation |l3lj |. was performed in Ref. |l39l |. see also 132] for more details. In this work, 
the DR expressions for the chiral 2PE potential in Eqs. 1)3. 36() . ()3.38|) were adopted. For M-,^, F-,^, qa 
and the constant 5 in Eq. (|3.33j) the following values were used : M^^ = 138.08 MeV, F^^ = 93 MeV, 
QA = 1.26, (5 = 0. For 01^3^4, the central values from Ref. |l54l | were adopted. No isospin-breaking 
interactions were included. Performing anti-symmetrization of the short-range part of the potential, the 
total number of independent contact terms was reduced to 9, see Eq. 1)3. 41(1 . The potential was multiplied 
by an exponential regulator function, which did not introduce any angular dependence, so that the 
contact interactions only contributed to S- and P-waves and to the mixing angle ei. The corresponding 
LECs were determined fitting each low partial wave separately to the Nijmegen PWA. As explained in 
sections 13.2.11 and 14 . 1 . H the subleading 2PE potential calculated using dimensional regularization shows 

*^*This number includes also the LECs which are usually taken from the pion and single-nucleon sectors, such as qa and 
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unphysically strong attraction at intermediate distances r 1 — 2 fm. The perturbative description of D- 
and F-waves based on DR potential leads to strong deviations froi n the data at rather low energies, see 
section 14.1.11 for more details. The nonperturbative analysis in Ref. ^ demonstrated, that phase shifts 
in S~, P- and D-waves could only be described simultaneously taking the momentum space cut-off A in 
the LS equation at least of the order of 1 GeV (unless one includes higher-order counter terms), which 
has to be compared with A = 500 MeV used at NLO. With such large cut-offs, the isoscalar central 2PE 
potential is already so strongly attractive that unphysical bound states in D- and in the lower partial 
waves are generated. In addition, phase shifts in D-waves, where the interaction is strong enough to 
produce bound states and no counter term s ap pear at N^LO, are strongly cut-off dependent. In spite of 
these difficulties, the N^LO potential from jl39l | allowed for a good description of the NN data, which was 
also visibly improved compared to the NLO results. In addition, NLO analysis with explicit A degrees 
of freedom was also presented in Ref. |l39l |. Using the large- A'^c value for the vrA^A coupling constant, 
the results were found to be similar to the ones at N^LO in the theory without isobars (including the 
appearance of the spurious bound states and strong cut-off dependence in D-waves). 

In order to get rid of the spurious bound states and enable few-nucleon calculations, a new version of 
the N^LO potential was introduced in i214j based on the numerically reduced values of the LECs 03^4: 
C3 = —1.15 GeV~^ and C4 = 1.20 GeV~^. This choice allowed for a fairly good description of the np 
data without unphysical bound states and using A ~ 500 MeV. The subleading 2PE potential then yields 
small corrections to the amplitude in most channels. The values of LECs 03^4 quoted above are, however, 
incompatible with vrN scattering. 

Entem and Machleidt studied the 2N system based on the chiral potential at N^LO (in the energy- 
independent formulation) obtained using DR and including contact terms up to N^LO UM- The values 
of the LECs c, were taken consistent with vrA^ scattering. The inclusion of higher-order counter terms 
at N^LO appears to be unavoidable in order to compensate for the unphysically strong attraction of the 
2PE potential obtained using DR. This, therefore, raises the question a bout the convergence of the chiral 
expansion for NN scattering in that framework. In their later work |l6,'j |. Entem and Machleidt also 
incorporated the 2P E potent i al at N^LO. Independently, the 2N system was studied up to N'^LO based 
on the SFR scheme Em, Ei3, 

[24^. In the following, we will show the results for various NN observables 



based on the most recent analyses of Refs. [162, 163]. 

As detailed in section ESI the chiral potential at N'^LO includes the IPE, 2PE and 3PE contributions, 
24 isospin-invariant contact terms as well as isospin-violating corrections. The leading 3PE potential 
turn s out to be rather weak (especially in the SFR framework), see section [3.2.21 and was neglected in 
[l63 . Il6^ . Both analyses use = 92.4 MeV and ^^(1 + S) = 1.29. In Table H we summarize the 
values of the LEC s c; and di adopted in these studies incpmparison with the ones extracted from vrA'^ 
scattering in Ref. |l54j | in th e case of 01^3^4 and in Ref. |218i] in the case of C2 and di. The LECs 02,3,4 
were fine tuned in Ref. |l63| |. wh ich r esulted, in particular, in the large value for C4 incompatible with 
ttN scattering. In the analysis of |l62l | , the central values fro m th e irN system are used for all LECs with 
the only exception of C3, for which the value dete rmined in j21ll | is adopted. This choice is on the lower 
side but still consistent with the result from |l54l ]. Usi ng th is value for C3 turns out to be important at 
N^LO in order to properly describe the ^Pq phase shift [243]. 



Isospin-breaking corrections are treated diffe rently in the analyses of |163l | and |162l ]. In |163l ]. they are 
incorporated following the lines of Ref. |24fll ]. Specifically, the authors of |l63l ] include the pion mass 
difference in IPE, see Eq. ()3.53() . the Coulomb potential in pp scattering, pion mass difference in the 
order u = 2 2PE potential as defined in Eq. (|3.59)) . the 7r7-exchange potential in Eq. (|3.52j) and th e tw o 
lowest-order isospin-breaking contact terms proportional to /?^^q and /9"^g in Eq. 1)3. 63() . In Ref. 



the long-range isospin-breaking corrections are treated in the same way as in the Nijmegen PWA |10SI |. 
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Cl 


C2 


C3 


C4 


di + d2 


d3 


(is 


dii — di5 




-0.81 


2.80 


-3.20 


5.40 


3.06 


-3.27 


0.45 


-5.65 




-0.81 


3.28 


-3.40 


3.40 


3.06 


-3.27 


0.45 


-5.65 


ttN 


-0.81(15) 


3.28(23) 


-4.69(1.34) 


3.40(4) 


3.06(21) 


-3.27(73) 


0.45(42) 


-5.65(41) 



Table 4: LECs used in the N^LO analyses NN" (Ref. [163]) and NN^ (Ref. [16; 



obtained from vrA^ scattering 1541 . 218| |. The Cj (di) are in units of GeV ^ (GeV 



compared to the values 



i.e. are based on the pion mass difference in IPE and the electromagnetic interactions in Eq. (|3.51|) . 
In addition, the two leading contact interactions proportional to /3^^q and /3"^q were included. Since 
phase shifts from the Nijmegen PWA were used in [162 ] to fix the values of the unknown LECs, treating 
isospin-breaking corrections differently from Ref. jl08| would lead to inconsistency.*^^ 



Also the treatment of the relativistic effects is rather different in [1631 ] and [1621 ]. The work of [1621 ] is 
based on the relativistic Schrodin ger /L ippmann-Schwinger equation as explained in section rj.2.41 *^^ On 
the contrary, the analysis of Ref. [l63i ] is based on the nonrelativistic Schrodinger equation a nd uses th e 
static IPE potential and the 1/m- and 1/m^-corrections to the 2PE potential from Refs. [l42l . [l65l ]. 
where no particular dynamical equa tion was specified. It is, therefore, not clear whether the relativistic 
corrections to the potential used in |l63| are consistent with the nonrelativistic Schrodinger equation. 

Further differences between the two analyses can b e attributed to different regularization schemes, fitting 
procedures and error estimations. In Ref. [l62j ]. SFR with the cut-off A was employed in ord er to 
regularize divergen t loop integrals entering the potential. Following the standard procedure, see e.g. 
139l . 147 . 2l4 . 248j ]. the resulting potential V{p, p') is multiplied with a regulator function 



V{p, p') - f\p) V{p, p') f^^ip') , (4.2) 

in order to remove the divergences in the LS equation, where the exponential regulator function 

/^(p) = expl-Z/A^] (4.3) 

was used. The following cut-off combinations (all values in MeV) were adopted in [l62| ]: 

{A, A} = {450, 500}, {600, 600}, {450, 700}, {600, 700} . (4.4) 

We remind the reader that while one, in principle, could further decrease the value of A (at the cost 
of the reduced accuracy), increasing A beyond ~ 650 MeV leads to the appearance of spurious bound 
states and one enters the regime, where the Weinberg power counting and "naive dimensional analysis" 
are not applicable, see the discussion in section 14.1.21 For A = 500 MeV, the value A = 600 MeV was 
already found to be close to its critical value corresponding to the appearance of a bound state. A typical 
dependence of the LECs on t he cu t-off A is exemplified in Fig. 1251 where we show the "running" of the 
LEC Cspi at N^LO from Ref. |248[ ]. The discontinuity in values of C3P1 for A ~ 700 MeV corresponds to 
the accumulation of the first spurious bound state in this channel. Further increasing the cut-off A leads 



*^^The major problem is that the np isovector phase shifts (except in the channel) in the Nijmegen PWA [lOSjl are not 
obtained independently from np data but rather extracted from the pp phase shifts by an appropriate change in the IPE 
potential and switching off the electromagnetic interaction. 

*^®As outlined in section l3.2.4l the relativistic Schrodinger equ ation can be cast in to equivalent nonrelativistic forms. The 
corresponding phase-equivalent potentials are also discussed in |l62j| . 
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Figure 25: "Running" of the LEG C^pi witli tlie cut-off A at N^LO. The cut-off A in the spectral function 
representation is fixed at the central value A = 600 MeV. 



to discontinuities in the values of C^pi corresponding to the accumulation of the second, third, etc. bound 
states. The limit cycle behavior of C^pi is typical for singular l/r" potentials'^^^ and is similar to the 
cut-off dependence of the strength of the contact 3N force observed in Refs. 



The analysis of Ref. [I63l | is based on the DR potential. The LS equation is regularized multiplying the 
potential with a regulator function 

fHp) = exp[ V7A'"] , (4.5) 

where the exponent 2n was chosen to be sufficiently large, so that the regulator generates powers which 
are beyond the order u = 4. Notice that n takes diffe rent values for various terms in the potential |25Cl| | 
(the actual values of n are not specified in Ref. |l6,l| |). The choice A = 500 MeV was adopted in this 
work. 

The 24 isospin-invariant + 2 isospin violating LECs corresponding to contact interactions (cf. Eqs. 
(|3.41l) and ()3.63|1 ') were fixed from a fit to the phase shifts with the orbital angular momenta I < 2. In 
Ref. '163], the fit was subsequently refined by minimizing the obtained from a dire ct co mparison with 
the data. The process of determination of the LECs is explained in detail in Ref. |l62| | and deserves 
some comments. First, we note that, in gen eral, one has to expect multiple solutions for the LECs. 
This problem has already been discussed in |139l | at NLO and N^LO. It is difficult to select out the 
true solution in the ^Sq channel, where five LECs /3fgQ, Ciso, Ciso, DI^q, -^150 ii^ed to be fixed from 
a fit to the Nijmegen pp and np phase shifts*^^ as well as for the ^S'i~'^-Di-channels, where eight LECs 
have to be determined simultaneously. In the simpler case of P-waves, where only two LECs need to 
be determined simultaneously (except in the coupled ^P2~^-^2 channels where one has 3 LECs), the 
choice of the solution is usually unambiguous if one makes use of the naturalness assumption. For 
example, the following two solutions for the LECs C^pi and D^pi were found for the cut-off combination 
{A, A} = {450, 500}: C3P1 = -0.6334x10^ GeV"'^, D3P1 = 4.2359 xlO"^ GeV'^ and C3P1 = 5.9620 xlO'^ 
GeV~^, -D3P1 = — 20. 6154 X 10^ GeV~^ |l62| |. Both sets of LECs lead to a similarly accurate description of 
the data. The first solution, however, fulfills the naturalness assumption while the second does not, i.e. the 
corresponding dimensionless coefficients are of the order ~ 10. Notice further that C3P1 = —0.6334 x 10^ 
GeV~'^ is close to the NLO and N^LO values (for the same cut-off combination) C3P1 = —0.4932 x 10^ 
GeV"^ and C3P1 = -0.7234 x 10^ GeV"^, respectively. The results for other P-waves are similar. No 



*^''The N^LO 2PE potential in the SFR scheme behaves at short distances as see Il47ll . 

*^*The LEG Piso was then obtained from the requirement to reproduce the nn S-wave scattering length. 
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multiple solutions arise in D-waves, where a single LEG needs to be determined in each channel. 



Let us now comment on the size of the obtained LECs. In general, the natural size for the LECs can be 
(roughly) estimated as follows: 

Att 47r Att 

^*~F2' ^ Z72A2 ' ^ F2A4 ' (4-6) 

where Alec is the scale entering the values of the LECs and the fac tor An results from the angular 
integration in the partial wave decomposition. All LECs determined in |l62l | except D\gQ and were 
found to be of natural size for all cut-off combinations. For example, the Z^j's expressed in the natural 
units defined in Eq. 1)4. 6() take the values in the range —1.1 . . . 2.0 if one adopts Alec = 500 MeV. For the 
LECs D\gQ and this range is, however, —5.2 ... — 11.2. Still, the higher-order contact interactions 
are suppressed compared to the lower-order operators at low momenta. For example, for the cut-off 
combination {450, 500} and p = p' = M^^, the contributions of the contact operators at various orders 
are given by: 

{'So\VZt{p, p)\'So) 



p,p'=M^ 



C^so + C,soip'+p") + [DlsoP'p" + Dlso {p'+p") 

p,p'=M„ 

0.091 + 0.057 + (-0.010 + 0.003)1 x lO'' GeV'^ . (4.7) 



For mo re de tails on the determination of various LECs and for their explicit valu es th e reader is referred 
to Ref. II62II. We further stress that the LECs of natural size were also found in 111911. 



We now turn to the discussion of phase shifts. Before showing the results, we would like to make a simple 
estimate for the expected theoretical uncertainty. Following the rules of the "naive dimensional analysis" , 
we expect for the uncertainty of a scattering observable at CMS momentum k at N^LO to be of the order 
~ (max[/c, MttJ/Alec)^- To provide a fair estimate, we identify the hard scale with the smallest value 
of the ultraviolet cut-off, i.e. we adopt Alec ~ 450 MeV. This value is consistent with the natural size 
of the determined LECs and yields the following estimations for the (maximal) theoretical uncertainty: 
- 0.5% for ^lab ^ 50 MeV and below, ~ 7% for Ei^b ^ 150 MeV and ~ 25% at E^iab ^ 250 MeV. One 
should keep in mind that the above estimations are fairly rough. For a more detailed discussion on the 
theoretical uncertainty the reader is referred to |248| |. 



The results at NLO, N^LO and N'^LO for S-, P- and D-waves and mixi ng angles ei and €2 from Ref. |1( 
are shown in Fig. [JHl in comparison with the PWA results of Refs. |loi.fl09. 222]. The bands correspond 
to variation of the cut-offs as specified in Eq. (|4.4|) . Low-energy observables should not depend on the 
cut-off(s) provided that all terms in the EFT expansion are included. In practice, however, calculations 
are performed at a finite order, so that some (small) residual dependence of observables on the cut-off(s) 
remains. One, in general, expects that this cut-off dependence gets weaker when higher order terms 
are included. While this is indeed the case for N^LO, the bands at NLO and N^LO are of comparable 
width. This does, however, not indicate any inconsistency as the following arguments show. The cut-off 
dependence of the scattering amplitude at both NLO and N^LO has to be compensated by inclusion of 
the counter terms (contact interactions) at order = 4 and higher. The contact interactions appear 
only at even orders i' = 2n, n > in the low-momentum expansion while pion exchanges contribute, in 
general, at both even and odd orders. Since the same contact terms enter the expression for the effective 
potential at NLO and N^LO, a similar cut-off dependence for the observables at these orders should be 
expected. The results shown in Fig. [2ni confirm these expectation. 

The uncertainty due to the cut-off variation in S-, P- and D-waves at N'^LO agrees well with the 
estimation given above. For example, at i?iab = 250 MeV, the P-wave phase shifts at N^LO deviate from 
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Figure 26: S-, and D-w aves np phase shifts. The dashed, hght shaded and dark shaded bands show 
the NLO, N^LO and N^LO [m2 results, respectively. The dashed line is the N^LO result of Ref. Il6i 
The filled circles (open triangles) depict the results from the Nijmegen multi-energy PWA |l08l . Il09l | 
(Virginia Tech single-energy PWA 222]). 

the data by an amount of up to ~ 8°, which has to be compared with the typical size of the P-wave phase 
shifts at this energy of the order ^ 25°. We further emphasize that all phase shifts at N^LO are visibly 
improved compared to N^LO (and NLO). The results at N^LO provide an accurate description of the 
data up to E'lab ~ 200 MeV. It is comforting to see that the bands at N^LO and N^LO overlap (except 
in some channels at higher energies). As pointed out before, one cannot expect the same for the NLO 
bands, which un deres timate the theoretical uncertainty at this order. Finally, we also show the N^LO 
results from Ref. |l63l | based on A = 500 MeV. The two N^LO analyses agree with each other in all cases 
except some minor deviations in the partial waves ^Pi and at intermediate and higher energies. 
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F- and selected higher partial waves, which are parameter-free at the considered orders, are shown in 
Fig. I27L In most channels, the predictions at N^LO are in agreement with the data. Contrary to the 
previously considered case of low partial waves, the bands do not get thinner at N^LO. This has to be 
expected due to the lack of the short-range contact terms in these channels. Such terms start to contribute 
to F-waves at N^LO {u = 6) and to G-waves at N'^LO {v = 8). Consequently, a similar uncertainty due 
to the cut-off variation should be expected up to these high orders in the chiral expansion. The largest 
deviation from the data is observed for the ^Gs partial wave. In this channel, the IPE and the leading 
2PE potentials are not sufficient to reproduce the PWA result at energies beyond ~ 100 MeV. The 2PE 
corrections at N^LO and N^LO improve the description of the data, but the effects are not big enough. 
This disagreement should, however, not be taken too seriously because of the exceptionally small size of 
the corresponding phase shift (m ore t han 10 times smaller in magnitude compared to other G-waves). 
The N^L O predictions from Ref. |l6,'j | are also plotted in Fig. 1271 They lie inside the N'^LO bands from 
Ref. jl62| | in all cases with the exception of the ^F^ and ^G^ partial waves. The relatively large deviations 
in th ese two channels might be caused by the different values of q (especially of the LEG C4) adopted in 

Once the NN phase shifts are calculated, all two-nucleon scattering observables can be obtained in a 
straightforward way using e.g. the formulae collected in (180i] . In Fig. I^Hlthe np differ enti al cross section 
and vector analyzing power at Eiab =25 and 96 MeV at N^LO and N^LO from [l^ are shown in 
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Figure 28: differential cross section and vector analyzing power at E'lah = 25 MeV (left panel) and 
-£'iab = 96 MeV (right panel). The Nijmegen PWA result is taken from [l09|. The cut-offs A and A are 
varied as specified in Eq. (|4.4() . For data see |l62l |. 



comparison with the data and the Nijmegen PWA results. This calculation is based on np partial waves 
up to J < 8 and does not include the magnetic moment interaction. It is comforting to see that N^LO 
and N'^LO results overlap and are in agreement with the Nijmegen PWA. The deviations in the analyzing 
power at forward direction are due to the neglected magnetic moment interaction. The small but visible 
deviations from the Nijmegen PWA for the differential cross section at forward and backward angles at 
N^LO are probably due to the neglect of partial waves with j > 8. 

We now regard the wp S~wave scattering length and effective range parameters. The results at various 
orders in the chiral expansion are summarized in Table |SI All observables are improved when going from 
NLO to N^LO and from N^LO to N^LO. The N^LO result for the scattering length fills a small gap 
between the N^LO prediction and the value of the Nijmegen PWA. For the pp ^Sq scattering length and 
effec tive range, the values Opp = —7.79 5 . . . — 7.812 fm and r^p = 2.73 . . . 2.76 fm were obtained at N^LO 
|Tfi3 |. which agree nicely with the data H: Upf = -7.8149±0.0029 fm and r^_= 2.769±0.014 fm. The 
"standard value" for the nn scattering length, a^^ = —18.9 fm, was used in \162\ as an input parameter 
in order to pin down the LEG Pi^q- The resulting prediction for the effective range, r^n = 2.76 . . . 2.80 
fm, agrees with the experimental value, r^rf = 2.75 it 0.11 fm. 

Finally, we woul d like to discuss various deuteron properties. The results at NLO, N^LO from |248l | and 
N^LO from |l62| | are collected in Table together with the experimental values. Notice that none of the 
deuteron properties were used in fits to determine the LECs. The predicted binding energy at N^LO is 
within 0.4% of the experimental value, which has to be compared with 1%-1.5% (~2%-2.5%) deviation 
at N^LO (NLO). The asymptotic S-wave normalization strength As is also visibly improved at N'^LO 
and deviates from the experimental (central) value by 0.3% as compared to ~1.1% (~1.9%) at N^LO 
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NLO [248] 


N^LO [248] 


N^LO [162] 


Nijmegen PWA 


a 150 [fm] 
nso [fm] 

0351 [fm] 
rssi [fm] 


-23.447... -23.522 
2.60... 2.62 
5.429... 5.433 
1.710... 1.722 


-23.497... - 23.689 
2. 62... 2. 67 
5.424... 5.427 
1.727... 1.735 


-23.585... -23.736 
2.64... 2.68 
5.414... 5.420 
1.743... 1.746 


-23.739 
2.68 
5.420 
1.753 



Table 5: np scattering length and effective ran ge for the ^5*0 and ^Si partial waves at NLO, N^LO and 



N'^LO compared to the Nijmegen PWA results 252l . 2 





NLO [248] 


N^LO [248] 


N^LO \162] 


Exp 


Ed [MeV] 


-2.171... -2.186 


-2.189... - 2.202 


-2.216... -2.223 


-2.224575(9) 


Qd [fm^] 


0.273... 0.275 


0.271... 0.275 


0.264... 0.268 


0.2859(3) 


m 


0.0256... 0.0257 


0.0255... 0.0256 


0.0254... 0.0255 


0.0256(4) 


VirYn. [fm] 


1.973... 1.974 


1.970. . . 1.972 


1.973... 1.985 


1.9753(11) 


As [fm-V2] 


0.868... 0.873 


0.874... 0.879 


0.882... 0.883 


0.8846(9) 


Pd [%] 


3.46... 4.29 


3.53... 4.93 


2.73... 3.63 





Table 6: Deuteron properties at NLO, N^LO and N'^LO compared to the data. Here, is the binding 
energy, the quadrupole moment, t]^ the asymptotic D/S ratio, \/{r^j^ the root-mean-square matter 
radius. As the strength of the asymptoti c S-wave normalization and the D-state pr obab ility. The 
data for E^ are from [255], for from |2 



2571 ]. for rjd from j256 



and f or A s from 12571 1. For the 

rms-radius we actually show the deuteron "point-nucleon" rms-radius from l2F 



(NLO). The results for the asymptotic D/S~ratio at all orders are in agreement with the data within the 
experimental uncertainty. The N^LO result for r]d agrees well with the one of the Nijmegen PWA ^252], 
rjd = 0.0253(2). No improvement is observed for the quadrupole momentum at N^LO, which shows an 
even larger deviation from the data compared to N^LO (6%-8% versus 4%-5%). One should, however, 
keep in mind that the effects of the two-nucleon currents and relativistic corrections to this observable 
were not accounted for in the calculation of Ref. 162]. The situation with the deuteron rms-radius is 
similar to the one with the quadrupole moment: one observes a larger deviation from the data at N^LO 
as compared to the NLO and N^LO re sults . Notice however that the N^LO result for this quantity is 
still within 0.5% of the value from Ref. 125311. 



Last but not least, we would like to mention that the N^LO analysis of Ref. '163] based on the cut- 
off in the LS equation A = 500 MeV was recently extended to A = 600 MeV |259| . In most partial 
waves, very close results for both choices of A with no visible differences for energies up to 300 MeV 
were reported. This seems to be in contradiction with the expectations based on the argument s and 
estimations presented above: even if A is chosen to be 500 MeV instead of 450 MeV as in Ref. |l62l ]. 
one can still expect only a very slow convergence of the low-momentum expansion at Siab = 300 MeV, 
which cor responds to the CMS momentum k ~ 375 MeV. It would be interesting to see how sensitive the 
results of 163l . EHi] are to a particular form of the regulator in Eq. 1)4. 5|) with n being chosen differently 



for different contributions to the potential. 
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LEG 


2PE (NLO) 


2PE (N^LO) 


a (NLO) 


Ci (n2lo) 


Bonn B 


Nijm-93 


Ciso 
Ciso 
Cssi 
C3S1 
C,i 

Cipi 
C3P0 
C3P1 
C3P2 


-0 004+°-°°° 

-n 570+0-036 

0.638+°:- 
-0 190+°°^2 

-0 067+°-°°^ 

u.uu/ _o oo5 

-0 425+°-°25 
246+°-°°^ 
-0 022+°-°°° 


-0 004+°-°°° 

_Q 443+0.078 
"-'•^^'-'-0.057 

-0 004+°-°°° 
-0 443+°-°^^ 

^•^^■^-0.057 

-0 090+°-°^=^ 

OOQg +0.003 

2 47 +0-032 

-0.044 

r) 1 C--I +0.020 
•^•-^^^-0.028 


-0 117+2-271 
"•-^-^ ' -0.042 

1 294+2-873 
■^•^^^-0.322 

1 or, +0.025 

U.iJO_Q Q21 

231+°-^^2 
^•^"^^-0.007 

n qoc; +0.000 

U.OZO_Q Q3g 

146+°-°°^ 
923+°-^^2 

-0 260+°-°°^ 

U.ZDU_g ggg 

-0 262+°-°32 

U.ZUZ g g.^3 


n 1 c;q +0.178 
-I 910 +0.408 

l.ZiO g 

-0 137+0-017 

"^•l"^' -0.027 
n c:9o+0.197 

U.OZO g g3g 

-0 395+0-007 

U.OJO_Q Q72 

126+°-°23 

U.lZD^g 

0.920 +J:?^^ 

-0.108l2.364 

-0.421 ins 


-0.117 
1.276 

-0.101 
0.660 

-0.410 
0.454 
0.921 

-0.075 

-0.396 


-0.061 
1.426 

-0.014 
0.940 

-0.343 
0.119 
0.802 

-0.197 

-0.467 



Table 7: The LECs Ci at NLO and N2L0 compared -with the results from the Bonn B and Nijmegen 93 
OBE potential models. Also shown are contributions from chiral 2PE as explained in text. The are 
in 10^ GeV-2 and the d in 10^ GeV"^. 



4.1.4 Resonance saturation for NN contact interactions 

In this section -we would like to confront the LECs determined from chiral effective field theory with 
the highl y su c cessf ul phenomenological/meson exchange models of the nuclear force following the lines 
of Refs. |248l . l26(]ll and restricting ourselves to N2L0. Consider, for example, the Bonn-B |26ll | and 
Nijmegen 93 j262l | potentials which are genuine one-boson-exchange (OBE) models. The long-range 
part of the interaction in these models is given by IPE (including a pion-nucleon form factor) whereas 
shorter distance physics is expressed as a sum over heavier mesons exchange contributions. For nucleon 
momentum transfer below the masses of the exchanged mesons, one can interpret such exchange diagrams 
as a sum of local contact operators with an increasing number of derivatives (momentum insertions). The 
LECs accompanying the resulting contact interactions are given within each model in terms of the meson 
masses, meson-nucleon coupling constants and corresponding form-factors. We now wish to compare 
these LECs to the ones enterin g the chiral NN potential in order to see, whether they can be understood in 
terms of resonance saturation |26nl |. In order to allow for a meaningful comparison with the OBE models, 
one needs to properly account for the chiral 2PE potential, which is absent in the OBE models. To achieve 
that, the 2PE potential at NLO and N2L0 is power expanded and the corresponding contributions to the 
LECs are identified, see |248l | for more details on this procedure and for explicit analytical expressions. 
The second and third columns in TableElshow the corresponding numerical results at NLO and N2L0 for 
the SFR cut-off A = 600 MeV. The indicated uncertainty refers to the cut-off variation A = 500 . . . 700 
MeV. The fourth and fifth columns contain the values of the LECs Ci at NLO and N2L0, where the just 
discussed contributions from 2PE have already been added. The numbers are presented for the cut-off 
combination {A, A} = {550, 600} with the uncertainties referring to the variations: A = 500 . . . 700 
MeV and A = 450 . . . 600 MeV (A = 450 . . . 650 MeV) at NLO (N2L0). Notice that in certain cases, 
Cj's show a rather strong cut-off dependence. This happens if the cut-off A becomes too large (i.e. close 
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to the critical value, at which spurious bound states arise) and one leaves the plateau-region for the 
corresponding LEG Ci(A). This situation is exemplified in Fig. 1251 Clearly, it only makes sense to 
discuss resonance saturation of the Cj's in the plateau-region of the first branch, where they only change 
modestly and where the effective potential is, at least, not strongly non-phase-equivalent to the OBE 
models. The last two columns in Table show the LECs as predicted by resonance saturation based 
upon the Nijmegen 93 and Bonn B potential models. One observes a remarkable agreement between the 
LEG values obtained from fit to NN phase shifts in the EFT approach and the ones resulting from the 
OBE models. For a related study based on a toy-model, the reader is referred to Ref. (263.] . 



4.1.5 Quark mass dependence of the nuclear force 

Since the absolute values of t he running up and down quark masses at the scale 1 GeV, niu — 5 MeV, ~ 
9 MeV, are rather small |l84l |. one expects that hadronic properties at low energies do not change strongly 
in the chiral limit (CL) of 0. This feature is crucial for the chiral expansion to make sense and 

is certainly true for the pion and pion-nucleon systems, where the interaction becomes arbitrarily weak 
in the CL and for vanishing external momenta. The situation in the few-nucleon sector is significantly 
more complicated due to the nonperturbative nature of the problem and also due to the fact that the 
interaction between nucleons does not become weak in the chiral limit. The M^-dependence of the 
nuclear force can naturally be studied in the chiral EFT framework. It is not on ly of academic interest, 
but also relevant for interpolating the results from lattice gauge theory, see also |ll4l |. For exampl e, th e 
S-wave scattering lengths have been calculated on the lattice using the quenched approximation j26J|. 
Another interesting application is related to imposing bounds on the time-dependence of fundamental 
couplings from the two-nucleon sector, as discussed in |265|]. 

Th e firs t, pioneering study of the NN system for vanishing quark masses was performed by Bulgac et 
al. |26fij | based upon the explicit M7r-dependence of th e IPE po t ential. More advanced studies and 



extensive dis cussi on on this topic can be found in Refs. |114l . Il34l . 12671 . 126a . I269l |. Here, we follow the 
lines of Ref. |l34l |. where the MTr-dependence of the nuclear force on the pion mass was analyzed at 
NLO based on the DR potential^^^ in the limit of exact isospin symmetry. We remind the reader that 
the potential at this order is given by the IPE and 2PE contributions and contact interactions with up 
to two derivatives or one M^-insertion as detailed in section 13.2.21 In addition, one has to include the 
corrections to IPE and the leading contact terms at the one-loop level, which lead to renormalization 
of the corresponding LEC s and therefore induce implicit quark mass dependence. These corrections are 
discussed in detail in |l34l |. For the IPE potential 

with (/ttAt being the pion-nucleon coupling constant, one has to account not only for explicit M^- 
dependence in the denominator, but also for implicit depen dence of the ratio g-Kw/^N on the pion 
mass. For an arbitrary value of the pion mass one obtains |l34l |: 

where = 1-26, F.^ = 92.4 MeV and M^r = 138 MeV denote the physical values of the nucleon axial 
vector coupling, pion decay constant and pion mass, respectively. Further, 1^, dig and die are LECs 
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As pointed out in section Pl.l.ll both DR and SFR lead to similar results at NLO. 
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Figure 29: Deuteron BE versus the pion mass. The shaded areas show allowed values. The light (dark) 
shaded band depicts the uncertainty due to the unknown LECs Ds,t {die)- The heavy dot shows the BE 
for the physical case M^r = -^tt- 



related to pion a nd pion-nucleon interactions. Following Ref. |l34l |. we use I4 = 4.3 p, die = "l-^S^Q^g 
GeV-2 lis, l27lj | and dig = —0.97 GeV ^. The co nstan t dig is determined from the observed value of the 



Golberger-Treiman discrepancy with Qt^n = 13.2 |272l |. Notice further that for the LEG die, we use an 



average of three values given in 27ll |. which result from different fits. The shown uncertainty is defined 



in the way to cover the whole range of values from ^21\- 

The remaining M^-dependence of the nuclear force at NLO is given by 2PE, see Eq. (|3.36|) . and by the 
short-range terms of the form 



F.^'i = Ml Ds + Driai • 0^2) - iPs + /3r(<?i • ^2)) In ^ ) , (4.10) 



where the constants Ps,T are given in terms of qa, F-,^ and Ct |273| |. All o ther contact terms d o no t 
depend on the pion mass and the corresponding LECs were adopted in |l34l | from the analysis of |l39l |. 
performed for the physical value = M,r- The essential difficulty in extrapolating the nuclear forces in 
the pion mass is due to the fact that the LEGs Ds, Dt cannot be fixed from the NN data.'^^'^ In order 
to proceed further, natural values for these LEGs were assumed in |l34l | 

Ds,T = ^ where - 3.0 < a5,T < 3.0 , (4.11) 

based on Alec — 1 GeV. This estimation is consistent with the size of othe r contact terms with the 
corresponding dimensionless coefficients ctj lying in the range —2.1 . . . 3.2 |26Cl ]. 

The resulting deuteron binding energy (BE) as a function of the pion mass is shown in Fig. 1291 The 
deuteron becomes unbound for ^ 1 70 M eV and is found to be stronger bound in the chiral limit 
with the BE = 9.6 ± 1.9 MeV IHi]. Here, the first indicated error refers to the uncertainty 

*^°They can be determined in the processes including pions such as e.g. pion-deuteron scattering. Such an analysis is 
however not yet available. 
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in the value of Dsg-^ and diQ being set to the average value die = —1-23 GeV~^ while the second 
indicated error shows the additional uncertainty due to the uncertainty in t he de termination of die as 
described above. Other deuteron properties in the chiral limit are discussed in |273l |. The resulting values 
for the two S-wave scattering lengths in the chiral limit are smaller in magnit ude a nd more natural: 
achCSo) = -4.1 ± 1.6+0° ™^ acL(^5'i) = 1.5±0.4+°3 fm. As pointed out in [lii], one needs lattice 
data for pion masses below 200 MeV to be able to perform a stable e xtrap olation to its physical value. 
Last but not least, we emphasize that the analysis of 268], see also |274l | for a related work, assumes 
a significantly larger uncertainty in the values of the LECs Ds^t and diQ, which is the main reason of 
different results for the_^i channel reported in that work. For a n ext ended discussion on this issue, the 
reader is referred to |269| |. We also emphasize that the work of |275l | seems to indicate the importance 
of including A's as explicit degrees of freedom in order to properly describe the quark mass dependence 
of the nucleon axial- vector coupling constant. Notice further that effects due to radiative pions might 
become important close to the CL. Finally, the possibility of an infrared renormalization group limit 
cycle in the 3N sector is discused in |276l |. 



4.2 The three- and four— nucleon systems up to N^LO 



We now turn to systems with three and four nucleons, which serve as an excellent testing ground for 
chiral forces and allow for a highly nontrivial check of the consistency of this approach since most of the 
unknown LECs are already determined in the NN system. 

Th e firs t quantitative application of chiral forces to Nd scattering in the 3N and 4N sectors was performed 
in |277t | at NLO. To that aim, the Faddeev-Yakubovsky equations were solved rigorously for the 3N and 
4N systems and the corresponding binding energies as well as various 3N scattering observables were 
computed. Since no 3NF appears at this order in the chiral expansion, the study of 277l | was entirely 
based on the 2NF yielding parameter-free results for the 3N and 4N systems. This work used the DR 
version of the chiral 2N potential with the LS cut-off varied in the range A = 540 . . . 600 MeV and 
demonstrated a good description of the nd elastic scattering data at E'lab = 3 MeV and ii^iab = 10 MeV 
as well as of some break-up observables at -Eiab = 13 MeV. The predictions for the triton and a-particle 
BE were found to be in a similar range as the ones based on phenomenological NN potentials. 



The NLO analysis of |277l | was extended to incomplete N^LO in |214l | where the 2N potential at N^LO 
was used without taking into account the corresponding 3NF. The first complete analysis of nd scat tering 
at this order including the 3NF was presented in Ref. ^194]. Both calculations of 214] and used 
the DR 2N potential with the numerically reduced values of 03^4 in order to avoid the appearance of 
spurious bound states, see section [4.1.31 for more details. In addition, certain nd scattering observables 
were studied in |278l ] using the 2N potential from 247 1 with no 3NF included. 



Recently, the 3N and 4N systems were reanalyzed at NLO and N^LO within the SFR framework |279 
We remind the reader that the SFR 2N potential at N^LO is based on the LECs Cj which a re co nsistent 
with vrA^ scattering. This is an important advantage compared to the previous study in Ref. 1941 ]. Notice 
th at th ese LECs also determine the strength of the 2PE 3NF. In addition, larger cut-off variation adopted 
in |l62l ] is expected to provide a more realistic estimation of the theoretical uncertainty. In the following 
sections we will present some results for 3N a nd 4 N observables based on the SFR 2NF. Further details 
and more results will be published elsewhere 127911 . In selec ted cases, we will compare the results based 
on the SFR and DR 2N potentials from Refs. Il62| and 121^1. F urther applications of chiral nuclear forces 
to 3N scattering can be found in Refs. 
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4.2.1 The formalism 



To describe the properties of the 3N and 4N systems, a corresponding Schrodinger equation needs to 
be solved. As described in section 14.1.31 we use a specific form for the regularization of the potential, 
which allows us to determined the LECs of the NN force partial wave by partial wave. This, however, 
implies that the interactions are strongly non-local. Therefore, a formulation in momentum space is most 
natural. 



In the past, the techniques were developed for solving r eliably the Schrodinger equation in momentum 
space using Faddeev- or Yakubovsky equations 28,4 28fi| |. For the 3N bound state problem, the Faddeev 



equations have the form 287 1 



t/; = Go t P V + (1 + Go t) Go v}^l (1 + P) V • (4.12) 

(i) 

Here V^23 is that part of the 3N force which singles out the particle i and which is symmetrical under 
the exchange of the other two particles. The complete 3NF is decomposed as 

V123 = V}^1 + VSI+VSI. (4.13) 

Further, ^p denotes the corresponding Faddeev component, t is the two-body t-operator. Go = 1/{E—Hq) 
is the free propagator of three nucleons and P is a sum of a cyclical and anticyclical permutation of the 
three particles. In case of nd scattering we follow the formulation described in Refs. i28fTi . i288l] and first 
calculate a quantity T related to the 3N break-up process via the Faddeev-like equation: 

T = tPcp+{l + tGo) Vill (1 + P) ^ + t P GqT + {1 + tGo) (1 + P) Go T , (4.14) 

where the initial state (p is composed of a deuteron and a momentum eigenstate of the projectile nucleon. 
The elastic nd scattering operator is then obtained as 

U = PGo^ + PT + Vl^l (1 + P) (1 + Go T) , (4.15) 

and the break-up operator via 

C/o = (l + P)r. (4.16) 



The s e eq uations are accurately solved in momentum space using a partial wave decomposition, see [2^ 
I290I . I291I for details. The partial wave decomposition of the chiral 3NF is given in the appendix of 
Ref H. 

Similarly, we use Yakubovsky equations (YE) to solve the 4N bound state problem. We rewrite the 
Schrodinger equation into two YEs and t hus decompose the wave function \I' into two independent 
Yakubovsky components (YCs): Vi and ip2 j28ft . .292 . The wave function then reads 

^ = (1 _ (1 + p) P3,)(i + p)^^ + (1 + p) (1 + p)^2 (4.17) 

and is again expressed with the help of the permutations P and P = P13P24 where P^ denotes transpo- 
sitions of particles i and j. In the case of a bound state it is, in principle, possible to solve directly the 
Schrodinger equation. The usage of two YCs is, however, advantageous since it naturally introduces two 
kinds of Jacobi coordinates which accelerates the convergence of the partial wave decomposition 293]. 
The YEs reduce to two coupled integral equations 

V'l = Got P[{l-Pu)^Pi+H + i'^ + Got)GoVSl^ (4.18) 
^-2 = Go t P [(1 - P34) V'l + H , (4.19) 



66 



which can be solved by similar techniques as the 3N problem. The high dimensionality (up to 10® x 10®) 
of the discretized YEs, however, requires massively parallel computers. 

Because these techniques are formulated in momentum space, the application of long-range interactions, 
like the pp Coulomb force, is a major difficulty. For the scattering problem, we therefore neglect the 
Coulomb interaction completely. At higher energies (above 50 MeV nucleon lab energy), it is expected 
to contribute mainly in forward direction. Here, the comparison of our results to the data should be 
less affected by the missing electromagnetic forces. For lower energies, we will compare our results 
to modified data, which have been corrected for t he electromagnetic interactions using calculations of 
Alejandro Kievsky |294l | based on the AV18 |295l | interactions with and without the electromagnetic 
forces. For the bound states, we include the pp Coulomb interactions in the t-matrices. Since the bound 
nucleons are confined to a small space region, we can put the Coulomb interaction to zero outside a radius 
of 10 — 20 fm, effectively making it short-ranged. Then the Fourier transformation of this interaction is 
nonsingular. The results are cut-off independent and numerically stable. 



4.2.2 Elastic Nd scattering 

Consider now elastic nd scattering. At N^LO, one has to take into account the 3NF which is discussed 
in section IH.3.11 Similarly to the 2N potential, the chiral 3NF behaves unphys icall y at large momenta 
and leads to ultraviolet divergences in the Faddeev-Yakubovsky equations. In i279| |. regularization was 
performed in the way analogous to the one adopted in the analysis of the 2N system: 



y^^-(p,g1p^r)-r(p,g1^^^^^(p,g1p^r)r(P^r), r{p,q) = e v ; , (4.20) 

where p and q {p' and q') are Jacobi momenta of the two-body subsystem and spectator nucleon before 
(after) the interaction. The regulator function f^{p,q) is chosen so that for g = it coincides with 
the function f^{p) in Eq. ()4.3() .^'^^ The results presented here and in what follows are based on the 
variation of the LS cut-off A in the range A = 400 . . . 550 at NLO and A = 450 . . . 600 at N^LO. The 
SFR cut-off A is varied in both cases in the range A = 500 . . . 700 MeV.#32 Specifically, the following 
cut-off combinations {A, A} were used at NLO and N^LO: 

NLO : {A, A} = {400, 500}, {550, 500}, {550, 600}, {400, 700}, {550, 700}, 
N^LO : {A, A} = {450, 500}, {600, 500}, {550, 600}, {450, 700}, {600, 700}. (4.21) 

As explained in section 13.3.11 the chiral 3NF at N^LO is given by the two-pion exchange, one-pion 
exchange with the pion emitted (or absorbed) by 2N contact interactions and 3N contact interactions. 
While the 2PE contribution given in Eq. ()3.73() does not introduce any new parameters, the two other 
terms depend on two new LECs, D and E, cf. Eq. (|3.74j) . which are not dete rmined in the 2N system 
and thus need to be fixed e.g. from the 3N data. As demonstrated in Ref. |l94l |. they can be determined 
using the '^H BE and the nd doublet scattering length "^Und, which are bona fide low-energy observables. 
In the present analysis, we use the coherent nd scattering length bnd instead of the ^a^d. This quantity is 
defined in terms of the doublet and quartet nd scattering lengths '^a^d and '^and and the neutron/deuteron 
masses as 

inn + rud 



bnd 



rrid 



1 2 1^4 

2 O-nd + 2 CLnd 



(4.22) 



* In Refs. a slightly different form of the re gula tor functions was employed. 

*^^The somewhat smaller values of A as compared to l 24a | were chosen in order to avoid the appearance of unnaturally 
large LECs, cf. Table □ 
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/i-EFT, N^LO 


NLO 


N^LO 


Exp 


and 




0.61... 1.19 


0.61... 0.63 


0.65 ± 0.04 [296] 


4 

and 


6.33 ± 0.05 


6.360... 6.366 


6.353... 6.362 


6.35 ± 0.02 [296] 


bnd 




6.656... 6.959 


6.669"^ 


6.669 ± 0.003 [297] 



Table 8: nd scattering lengths (in fm) in EFT in comparison with the data. For NLO and N^LO in the 
EFT with explicit pions, the cut-offs A and A are varied in the range given in Eq. ()4.21|1 . The value of 
bnd used as input is marked by the star. 



where ^a„rf is the nd quartet scattering length and nin {rrid) refers to the neutron (deuteron) mass. The 
coherent nd scattering length is much better known experimentally than ^a„rf, which allows to reduce the 
uncertainty in the determination of the LECs D and E. The resulting LECs are found to be of natural 
size, i.e. the constants cd and ce defined as 

are of order one. The only exception is given by the cut-off combination {600, 500}, for which the 
magnitude of cd is rather large {cd = —10.0 for Alec = 700 MeV). Notice further that the values for 
both LECs depend strongly on the cut-offs and differ significantly from the ones found in an earlier 
study |l94l |. Finally, we would like to emphasize that the BE and the central experimental value of 
hnd could be reproduced simultaneously only for {Aj, Aj} in Eq. (|4.2H) with i = 1,3,4. For two other 
cut-off combinations, it was not possible to find a solution for co and c^, which would describe both 
observables at the same time. In these cases, the scattering length bnd (the triton BE) was required to 
be reproduced exactly (as good as possible). We will discuss the resulting BE in section I4.2.4L More 
details on fixing the values oi cd and ce will be given in |279l |. 

In Table |H1 we summarize the EFT results for various nd scattering lengths in comparison with the 
experimental numbers. In the case of pionless EFT, the value shown in this table is taken from Ref. [h^, 
where more discussion on this approach and further references to earlier determinations of '^and can be 
found. 

With the LECs cd and ce being determined as described above, we are now in the position to predict 
various nd elastic scattering observables. In Fig. EOlwe show the results for the differential cross section 
(in the left column) and vector analyzing power (in the right column) at 3, 10 and 65 MeV. One observes 
a good agreement with the data at two lowest energies at both NLO and N^LO. This also holds true for 
tensor analyzing powers shown in Fig.|^for E'lab = 10 MeV. The only exceptions are given by the minima 
of T20 and T21 at N^LO, where some deviations from the data are visible. We remind the read er th at the 
data at energies 3 and 10 MeV have been corrected to account for the missing Coulomb force |294l |. The 
uncertainty due to the cut-off variation in the differential cross section is significantly reduced at N^LO 
compared to NLO at all energies considered. At 65 MeV, the results at NLO show large deviations from 
the data in the minimum of the cross section as well as for Ay. Similarly to the 2N system, the bands 
at NLO seem to underestimate the true theoretical uncertainty at this order, see section [4.1.31 for more 
discussion. The N'^LO results at 65 MeV are in agreement with the data for both observables. Notice 
that the deviation from the data for da/dQ in forward directions is due to the Coulomb force which is 
missing in the calculations. At this energy, the uncertainty at N^LO appears to be rather large. For 
the tensor analyzing powers at the same energy, the situation is similar: while the NLO results deviate 
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Figure 30: Differential cross section (in mb/sr) and vector analyzing power for elastic nd scattering at 
3 MeV (upper panel) 10 MeV (middle panel) and 65 MeV (bottom panel) at NLO (light -shaded bands) 
and N^LO (dark-shaded bands) in the SFR framework. The bands correspond to the cut-off variation 
as specified in Eq. (|4.2H) . For data see |l94t |. 



significantly from the data, the N^LO predictions are in agreement with the d a,ta b ut the uncertainty 
due to the cut-off variation is large. A more detailed discussion will be given in |279| |. 

Let us now discuss the nucleon vector analyzing power Ay, which is the most problematic observable in 
nd elastic scattering at low energy. This particular observable is underpredicted in the maximum by an 
amount of ^ 25 . . . 30% by modern high-precision nuclear potentials, which is known in the literature as 
Ay puzzle Egi,^^^. Augmenting the NN potentials with 3NF models such as the TM99' 3NF |300] or the 
Urbana-IX 3NF ^SQi] , which are frequently used in modern few-body calculations, does not substantially 
improve the description of this o bserv able. The only exception is given by the phenomenological spin- 
orbit 3NF introduced by Kievsky |302l | . which allows to describe the data. Similar discrepancies (but less 
pronounced compared to Ay) are also observed for the tensor analyzing power iTn. As demonstrated in 
Fig. I30[ The NLO result for Ay is in agreement with the data at 3 MeV and even slightly overpredicts 
the data at 10 MeV. Very similar results based on the DR NN potential at NLO were obtained in 
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Figure 31: Tensor analyzing powers for elastic nd scattering at 10 MeV at NLO (light-shaded bands) 
and N^LO (dark-shaded bands). The bands correspond to the cut-off variation as specified in Eq. (|4.2ip . 
For data see 



214l . l277| While this looks encouraging, one cannot conclude that the Ay-puzzle has been solved. This 



observable is well know n to be very sensitive to the spin-orbit 2NF and, therefore, to the triplet P 



278l . 285l |. which need to be reproduced accurately in order to h ave c onclusive results. It 



waves, see e.g. 

is instructive to look at the spin-orbit phase shift combination A^s defined as [2 

^LS = ^{ 26sp^ - 36,p^ + 55sp^) . (4.24) 

In Table EJ we show the results for this quantity at NLO and N^LO compared to the ones from Nijmegen 
PWA. Clearly, the spin-orbit force at NLO is enhanced compared to Nijmegen PWA, which also explains 
the enhancement for nd Ay at this order. Notice that the overestimation of A^^ at NLO (and, to a less 
extend, also at N^LO) is largely due to the failure to properly describe the partial wave, cf. Fig. 1261 We 
emphasize, however, that the quantity A^s is more accurately reproduced at N^LO, where the calculated 
nd Ay is in a reasonable agreement with the data (although the uncertainty due to the cut-off variation 
appe ars to be quit e sizable). A more detailed discussion including the role of the 3NF will be given in 
j279|, see also |303| for a related earlier work. 

Finally, we emphasize that at low energy the results for nd elastic scattering observables at NLO and 
N^LO in both SFR and DR |2l3] frameworks are very similar. The strongest differences are observed 
for Ay and iTn, where the N^LO correction is larger in the DR approach. Very different values of 03^4 
adopted in these analyses, which determine the strength of the 2PE 3NF, have only a little impact on the 
considered nd elastic scattering observables. At higher energies such as 65 MeV, the differences between 
the two sets of calculations, howeve r, be come quite significant. F or fur ther results in nd elastic scattering 
based on the NN potential of Ref. |247i ] the reader is referred to \21 
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Slab (MeV) 


NLO 


N^LO 


Nijmegen PWA 


10 


0.201... 0.210 


0.183. . . 0.199 


0.202 


20 


0.694... 0.721 


0.622... 0.673 


0.641 


30 


1.46. . . 1.51 


1.28... 1.39 


1.25 


40 


2.43... 2. 54 


2. 10... 2. 28 


1.97 


50 


3.57... 3.76 


2.97... 3.30 


2.73 



Table 9: Spin-orbit combination A^^ for np phase shi fts ( in degrees) at various energies at NLO and 
N^LO in comparison with the result of Nijmegen PWA [l09|. The cut-offs are chosen as in Eq. 1)4. 21() . 



4.2.3 Nd break up 

Let us now switch to break-up observables. In the fohowing, we will show some results at two energies, 
13 and 65 MeV, where a lot of pd data exist. Since we do not have reliable Coulomb corrections in the 
case of break-up, we show the non-corrected pd data in comparison to our nd calculations. 

In Fig. 1211 we present the results for a few often investigated configurations, the space-star (SST), a final 
state interaction (FSI) peak configuration and a quasi-free scattering (QFS) configuration, respectively. 
For a general discussion on various break-up observables and configurations the reader is referred to 
j285l ] . As demonstrated in Fig. 1321 the results at 13 MeV are very robust and essentially the same at 
NLO and N^LO in both the SFR and DR approaches. The description of the configuration dominated 
by FSI peaks is rather accurate (for a more elaborated procedure the angular openings of the detectors 
have to be taken into account, see ^285i] ) . The present theory for the break-up configuration including 
a QFS geometry fails in the central maximum. This might be due to missing Coulomb force effects. 
The third configuration, the so called space-star, is one of the long standing puzzles of 3N scattering, 
see e.g. H. Similar to the phenomenological interactions, we even fail to describe the nd data. The 
strong deviation between the nd and pd data indicates the importance of the Coulomb effects for this 
particular configuration. It is still an open question whether the pd data can be reproduced by the theory 
if Coulomb force is taken into account. 

At 65 MeV we show in Fig.|2Hlthe results for the cross section and nucleon vector analyzing power for two 
selected config urat ions. These and other configurations were studied in t he co ntext of phenomenological 
nuclear forces |306l | as well in the chiral EFT based on the DR potential |l94l ]. The situation at 65 MeV 
seems, in general, to be promising. In most cases, a clear improvement in the description of the data is 
observed when going from NLO to N-^LO. We stress, however, that the uncertainty due to the cut-off 
va riatio n is rathe r larg e at this energy. More results for break-up observables at this energy will be given 
in H, see also HiJ 

In addition, we would like to mention some further recent studies, where the results based on chiral 
EFT were shown. Differential cross section of the ^H(p, pp)n reaction at the proton energy 16 MeV in 
three kinematical configurations, the np FSI, the co-planar star, and an intermediate-star geometry, was 
analyzed based on the conventional nuclear forces as well as on chiral EFT at NLO and (incomplete) 
N^LO and N^LO. Various proton-to-proton and proton-to-deuteron polarizatio n tra nsfer coefficients in 
d{p, p)d and d{p, d)p reactions at the proton energy 22.7 MeV are considered in |284'] . For these observ- 
ables, the restriction to the forces in NLO is shown to be insufficient. At N-^LO a satisfactory description 
of the data is observed, similar to the one obtained with the (semi) phenomenological interactions. In 
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Figure 32: nd break-up cross section in [mb MeV~^ sr^^] along the kinematical locus S at 13 MeV in 
comparison to predictions at NLO (light shaded band) and N^LO (dark shaded band) in the chiral EFT. 
In the left (right) panel, the results based on the SFR (DR) scheme are shown. In the upper row a final 
state interaction configuration is depicted, in the middle one a quasi-free scattering configuration (both 
in comparison to pd data) and in the lower one a space star configuration (upper d ata nd, lower data 
pd). The precise kinematical description and references to data can be found in Ref. |285l |. 



addition, differential cross section of th e deuteron-proton break-up reaction at the deuteron energy 130 
MeV was studied recently in Ref. |283l | . In this work, which is mainly focused on the 3NF effects, the 
results based on modern NN potentials combined with 3NF models and on chiral EFT at N^LO are com- 
pared with the new high-precision data for 72 kinematically complete configurations. The description of 
the data at N^LO is found to be of a similar quality as the one based on realistic high-precision nuclear 
force models. 
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Figure 33: nd break-up cross section in [mb MeV~^ sr~^] and nucleon analyzing power along the 
kinematical locus S (in MeV) at 65 MeV in comparison to predictions at NLO (light shaded band) and 
N^LO (dark shaded band) in chiral effective field theory. The cut-offs A and A are chosen as in Eq. (|4.2H1 . 
The left and right panels sho w th e symmetric space star and symmetric forward star configurations, 
respectively, pd data are from |3r)5j | . 



4.2.4 Bound states 



The results for the triton and a-particle BEs are summarized in Table 1101 These observables are, in 
general, very sensitive to small changes of the interaction, as they come out as the difference of the large 

kinetic and potential energies^ A s a c onsequence, we found a rather large dependence of the BEs on 

the cut-off at NLO, see also Hii, |277| for similar results obtained in the DR scheme. At N^LO the 



BE, which was used as input in order to fix the LECs in the corresponding 3NF as explained in section 
14. 2. 21 is exactly reproduced for three cut-off combinations. In the two other cases, the deviation from 
the experime ntal value is less than 1%. Because of the strong correlation of the 3N and 4N BEs known 



as Tjon-line |3f)7l | . one can expect a good description of the a-particle BE. We, however, emphasize that 
3NFs break this correlation, see e.g. |30g]. Indeed, we observe a significant dependence of the a-particle 
BE on the value of the LEG cd with ce being fixed to reproduce the triton BE. For a recent work on 
the correlation between the 3N and 4N BEs in the context of pionless EFT the reader is referred to 
Is^. The predicted values for the a-particle BE are within 5% of the experimental value for all cut-off 
combinations considered. 

We would like to emphasize that the large numerical value of cd for the cut-off combination {600, 500} 
is reflected in the large expectation value of the 3NF in the triton and a-particle. The situation is 
similar for the cut-off combination {600, 700}, although the expectation values of the 3NF are smaller 
in magnitude. This might indicate that the cut-off A = 600 MeV is already too close to its critical value. 
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NLO DE [1941 


N^LO DE [1941 


NLO SFE 


N^LO SFE 


Exp 




-7.53... -8.54 
-23.87... - 29.57 


-8.68* 
-29.51... - 29.98 


-7.71 ... - 8.46 
-24.38... - 28.77 


-8.48... - 8.56* 
-27.77... - 29.61 


-8.482 
-28.30 



Table 10: "^H and ^He BEs (in MeV) at NLO and N^LO of the chiral expansion (for the cut-off range 
considered throughout) compared to experimental BEs. The values of Ei^ used as input is marked by 
the star. 



where spurious bound stat es appear and the naturalness assumption is violated. Further details on this 
topic will be given in |27 



4.3 More nucleons 



Due to the fast increase of the comput ational ab ilities, one is now able to solve the Schrodinger equation 
for light p-shell nuclei including 3NFs 3091 . 3ic| | . 



Using phenomenological forces, it has been established that the binding energies and spectra do depend 
on th e structure of the 3NF, even if models describe the triton and a-particle binding energy equally well 
|31lj |. Moreover, they also depend on the isospin T = 3/2 components of the 3NF, to which the previously 
discussed 3N and 4N observables are not sensitive. This makes the application of chiral interactions to 
light nuclei even more interesting. 

Chiral interactions are low momentum interactions. As discussed throughout this review, the unknown 
short distance part of the force is absorbed in a tower of contact terms. It turns out that one obtains 
a decent description of the NN data using rather small cutoffs. The experience with traditional models 
indicated, however, a need for hard cores not only in the NN interaction, but also in the 3NF in order 
to prevent strong overbinding in systems beyond the s-shell. This seems to be in contradiction with 
the basic EFT philosophy. In the following, we will discuss some results for systems with six and seven 
nucleons, which do not show any indication of unphysically increasing densities or binding energies. 

Let us begin with ^ Li. I n Fig. |^we show the parameter-free results for the ^Li binding and excitation 
energies from Ref. |312| obtained within the no-core shell model framework and based upon the DR 
NLO and N^LO chiral forces from Ref. |214l |. Going from NLO to N^LO one observes a reduction of 
the cut-off dependence. At N^LO, the uncertainty due to the cut-off variation is 1.8 MeV or 5.7% (170 
keV or 7.6%) of the binding (excitation) energy. Notice that for the BE, the estimated correction of the 
experimental value due to the missing CIB is also shown in Fig. 1341 Calculation based on the N^LO 
Idaho NN potential 1 163'] and N^LO 3NF led to similar results |312l |. In this study, the LECs cd and ce 
were fixed from the and ^He BEs which yielded two different sets: cd = —1.11, ce = —0.66 (3NF-A) 
and CD = 8.14, ce = -2.02 (3NF-B). 

We now turn to ^Li. Here, the only calculations performed so far are based on the N^LO Idaho NN po- 
tential augmented with the N^LO 3NF. The results fo r the binding energy are summarized and compared 
to other calculations and the experiment in Table ITTI |313l |. As one can expect, ''Li is underbound for the 
NN interaction only. Both sets of 3NFs provide more binding. However, in both cases the final binding 
energy result is still short of the experiment by 1.2 and 2.5 MeV, respectively. This slight underbinding 
is encouraging in view of the general expectation that strong repulsion at short distances is required to 
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Figure 34: ^Li bin ding (left panel) and excitation (midle pane l) energies at NLO and N^LO based upon 
the DR 2NF from |214l | and the corresponding 3NF from |l94l | in comparison with experimental values. 
The bands correspond to the A variation between 500 and 600 MeV. Right panel shows the results for 
^Li based upon the N^LO Idaho 2NF combined with the N^LO 3NF (sets A and B) in comparison with 
experimental values. Figure courtesy of Andreas Nogga. 



avoid a collapse of nuclei. Clearly, we do not observe any sign of overbinding, though neither the NN force 
nor the 3NF have a strong repulsive core. Both interactions are very soft, yet the binding energies are 
reasonable. So far, it wa s believed that only the addition of a repulsive core, like in the Urbana-IX and 



Illinois models j30ll . l31l| . can cure this overbinding problem. Here, it is demonstrated that the additional 



structures of chiral 3NFs also prevent overbinding. 

Table ITT] also shows the radii obtained in this way in comparison with the experimental values and other 
calculations. The results for the chiral interactions are comparable to the ones based on phenomenological 
models Urbana-IX and Illinois. No indication is observed that soft, chiral interactions fail to saturate 
nuclear systems with a realistic binding energy and density. 

Finally, the predictions for the excitation energies are summarized in Fig. 1341 All combinations of the 
interactions, Idaho N'^LO NN force alone, with 3NF-A or 3NF-B, do predict the right ordering for these 
states. The splitting of the 3/2~ and l/2~ states is small. The agreement to experiment for 3NF-A 
seems to be superior. Because the splitting itself is very small, this might be accidental. More significant 
deviations are observed for the 7/2~ and 5/2~ states. Both, the position of this multiplet and the splitting 
are strongly affected by the 3NFs and the agreement with the experimental results is clearly best for 
3NF-B, which is opposite to the case of BE. Clearly, further studies are needed in order to clarify the 
situation. As discussed before, the two N'^LO NN interactions available at present use different values 
for the LECs Cj wh ich d etermine the strength of the 2PE 3NF. The results presented here are based on 
the choice of Ref. and it is conceivable that the detailed description of the binding energies and 





NN 


NN + 3NF-A 


NN + 3NF-B 


AV18+Urbana IX 


AV18+IL2 


Exp 


Egs [MeV] 


34.6 


38.0 


36.7 


37.5 


38.9 


39.2 


r [fm] 


2.40 


2.19 


2.31 


2.33 


2.25 


2.27 



Table 11: Comparison of the ground state binding energy results Egs and the point proton rms radius 
r for ''Li for chiral interactions and several phenomenological combinations to the experiment. 
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spectrum is affected by the choice of Cj's. For more details the reader is referred to |313| . 

To summarize, the first results for the 6N and 7N systems based on chiral forces look promising. No hard 
repulsive core is found to be necessary to provide realistic binding energies and densities of p-shell nuclei. 



5 Miscellaneous omissions 



In this section we give a list of related topics which were not discussed in this work. We stress that this 
list and the given references are not meant to be complete and should merely provide the reader some 
guidance to further studies. 



Electroweak and pionic probes in the nuclear environment. The exchange vector a nd axial- vector 

using 



NN currents were considered in chiral EFT by Park, Min and Rho at one-loop level 3141 . 



the Feynman graph technique and dimensiona l reg ularization. Most of the practical a pplication s 



including Co mpt o n sc attering on the deuteron |31fil | , radiative up cap ture (at t hresho ld) 31RI . 1311 
solar fusion |3isl . l31ij | and other solar neutrino-deute ron reactions 13201 . l32lj |. hen l322 1 and hen 
|323l | processes, pi qn-d euteron scattering at threshold |ll3l . l223l . |32 4|, pion photo- |325l . l326l | and 
electroproduction 327 1 . pion production in NN collisions, see [iis] and references therein, and 



others were performed in the so-called "hybrid" approach. In such a scheme, the interaction kernel 
is derived within chiral EFT while the wave functions for few-nucleon initial and final states are 
calculated using phenomenological potentials. Certain reactions were already studied in a consistent 
way base d on the few-nucleon wave functions o btain ed in chiral EFT, see e. g. Refs. liii,^,! 



33 ll . l332l | , as well as in the pionless framework |335 



• Nuclear parity violation. A systematic study of nuclear pa ritv- violation within the framework of 
effective field theory was recently carried out by Zhu et al. |33J|, who derived the parity-violating 
IPE and 2PE potentials and discussed the ways of fixing the corresponding unknown LECs. 

• Nuclear forces and the large Nc limit. The role of the large A^c~liiiiit of Q CD for the nucleon-nucleon 
interaction wa s original lv addressed in the seminal work by Witten j335l | and then more recently by 
Kaplan et al. |336l . 133'^ . The consi steiicy of the large A'^c'-limit with the meson-exchange picture 



of the nuclear force is discussed in 33Fjl . 339l . 3401 . l34l]. 



Chiral effective field theory and nuclear matter. Applications ofchiral EFT to th e nuclear many- 
body problem were considered by several groups, see |342l . l343l . l344l . 134,4 for some recent 
references. This is presently an active research field. 

Nucleons in (partially) quenched QCD. The NN potential in quenched and partially quenched QCD 
is discussed in [267l |. Application of a partially-quenched extension of an effective field theory 



347l |. These studies will help to relate lattice 



to nucleon-nucleon scattering was performed in 
simulations in the 2N sector to experimental data. 

Low-momentum NN interaction Viow k ■ An effective interaction Vi n,,, k acting in the Hilbert space 
of low-momentum modes was first explicitly constructed in Refs. |263l . 1348.] using the method of 
unitary transformation and based on a potential of the Malfiiet-Tjon type. A different method 



was employed in Refs. 349l . 



to construct Viow k for various high-precision NN potentials. The 
universality of the resulting low-momentum interactions and the implications for the nuclear many- 
body problem ar e discusse d in |35ll |. Some recent applications of VJo^ ^ in the few-nucleon sector 
are presented in 352l . 353| | . 
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6 Outlook 



In this work we outlined in detail the structure of the nuclear forces in the framework of chiral effective 
field theory and discussed recent applications in the few-nucleon sector. We also focused on some related 
topics including isospin violating effects and chiral extrapolations of the two-nucleon observables. 

In the future, these studies should be generalized in various ways. First, it is important to extend the 
N^LO analyses jl62l . Il63l ] of the 2N system to few-nucleon systems which requires the derivation and 
numerical implementation of the corresponding N'^LO contributions to the 3N and 4N forces. It remains 
to be seen whether the presently observed difficulties in the theoretical description of certain low-energy 
3N scattering observables like e. g. the nucleon vector analyzing power Ay can be overcome at this 
order in the chiral expansion. Second, the electroweak nuclear current operators should be constructed 
to the same accuracy as the nuclear forces and applied to a rich variety of electroweak reactions in 
the nuclear environment without using the "hybrid" approximation. Furthermore, given the recent 
theoretical progress in understanding isospin and parity-violating corrections to the nuclear force within 
the chiral EFT framework, a systematic study of these effects in few-nucleon systems should be pursued. 
In addition to these fairly straightforward extentions, a further effort is called for to achieve a better 
understanding of the nonperturbative renormalization in the context of the few-nucleon problem. This 
issue is being currently investigated by several groups. Increasing the range of applicability of the chiral 
EFT approach forms another challenging direction of future research. This might require the inclusion of 
A-isobar as an explicit degree of freedom. One particular difficulty in such an approach is related to our 
lack of knowledge of the values of the corresponding LECs, since applications of the EFT with the A's 
in the single-baryon sector were not yet performed to the same detail as in the A-less theory. Finally, it 
would be highly desirable to incorporate in the EFT more constraints from QCD using for example the 
large- Ai'c expansion or lattice gauge theory. 
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A 2PE potential from the third order ttN amplitude 

In this appendix, we give the expressions for the 2PE 2NF resulting from the d iagra ms which contain 
the third order pion-nucleon amplitude, see Fig. d and were obtained by Kaiser jl4f)l |. It appears to be 
convenient to use the (subtracted) spectral function representation: 
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Wc,s(q) = [ 



2a^ f'°° 

J2AU 



(A.l) 



For the spectral functions Pi{fi) {r]i{fi)) one finds |l45l |: 
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2M^ + 1 — — In ^ 



2fi 
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where we have introduced the abbreviations 

1 



r=-V/x2-4M^ 



and 



1 + 



1^ 2 2 



2 \ 3/2 



In 



(A.2) 

(A.3) 
(A.4) 



We use the scale-independent LECs di, ^2, ^3, ^5, di4 and dis defined in |218l |. In the limit A ^ 00 
some of the integrations in E qs. (jA.ll) with the spectral functions given in Eqs. (|A.2|1 have been carried 
out analytically in Ref. 211 1. 
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